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iX  OlSTRISUTION/AVAILASILITV  OF  ABSTRACT 
UNCLASSIFIEO/UNLIMITEO  is  same  as  rft.  C^^c  users 


ESa  name  of  responsible  individual  ,  ^ 

Li. Cot. 


TjEzaas 


ABSTRACT 


This  report  presents  a  brief  summary  of  the  principal  results  obtained  in 
a  research  program  on  damage  development  in  fibrous  composite  laminates.  The 
following  technical  subjects  are  described:  (1)  Effect  of  transverse  cracks 
and  fiber  breaks  on  stiffness  changes  in  unidirectional  and  laminated  plates. 
(11)  Effect  of  ply  thickness  on  initial  failure  and  on  progressive  cracking  in 
brittle  matrix  laminates,  and  (111)  Analysis  of  progressive  cracking  in  metal 
and  polymer  composite  laminates.  Work  in  progress  is  discussed  as  well. 
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INTRODUCTION 


This  research  project  was  conducted  as  a  cooperative  effort  of  two  Inves¬ 
tigators.  Dr.  George  J.  Dvorak  was  the  principal  Investigator  of  the  program 
at  Rensselaer  Polytechnic  Institute,  which  was  the  primary  contractor.  Dr. 
Norman  Laws  was  the  principal  Investigator  of  the  part  of  the  program  sub¬ 
contracted  from  RPI  to  the  University  of  Pittsburgh. 

The  program  addressed  several  basic  problems  In  damage  mechanics  of  both 
brittle  and  ductile  fibrous  con^ioslte  materials,  and  produced  many  new  results. 
The  principal  accon^llshments  can  be  summarized  as  follows: 

1.  Formulation  of  self-consistent  and  bounding  techniques  for  a  deterministic 
evaluation  of  stiffness  changes  caused  by  transverse  cracks  or  by  fiber 
breaks  in  plies  of  any  orientation  In  fibrous  laminates.  The  accuracy  of 
stiffness  predictions  obtained  from  the  self-consistent  method  was  verified 
by  comparison  with  available  and  our  own  experimental  data. 

2.  Evaluation  of  the  effect  of  ply  thickness  on  Initial  failure  and  on  pro¬ 
gressive  cracking  in  polymer  matrix  laminates  under  incrementally  Increas¬ 
ing  load.  The  significant  conclusion  obtained  in  the  analysis,  and  in  the 
related  experimental  confirmation,  was  that  damage  development  In  such 
brittle  matrix  laminates  could  be  retarted  or  completely  eliminated  by 
keeping  the  thickness  of  each  ply  to  a  minimum. 

3.  Analysis  of  progressive  cracking  In  metal  and  polymer  composite  laminates. 
In  this  effort  we  were  able  to  predict  stiffness  changes  which  were  ob¬ 
served  experimentally  under  monotonlcally  Increasing  loads  In  glass  and 
graphlte/epoxy  laminates.  Also,  we  analyzed  the  pi.  cess  of  fatigue  damage 
development  In  metal  matrix  B/A£  laminates,  which  Is  caused  by  cyclic 

plastic  straining  of  the  matrix.  We  found  that  damage  In  this  system 
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served  as  a  shakedown  mechanism  which  allowed  the  composite  to  resume  an 
elastic  deformation  response  in  the  saturation  damage  state.  An  incremen¬ 
tal  procedure  was  formulated  to  find  the  stiffness  changes  associated  with 
saturation  damage  states  reached  under  any  steady-state  loading  program. 
The  predicted  stiffness  changes  were  found  to  be  in  excellent  agreement 
with  experimentally  measured  magnitudes. 


These  principal  findings,  together  with  other  related  results  are  de¬ 
scribed  in  the  sequel. 


2 


SIGNIFICANT  ACHIEVEMENTS 


2. 1  Loss  of  Stiffness 

The  first  important  technical  results  obtained  in  the  course  of  this  re¬ 
search  relate  to  the  loss  of  stiffness  of  unidirectional  composites,  or  indi¬ 
vidual  plies  in  a  laminate,  due  to  transverse  matrix  cracks  and  due  to  fiber 
breaks  and  consequent  aligned  penny-shaped  cracks.  The  general  analysis  of 
Laws,  Dvorak  and  Hejazi  [1]  was  simplified  and,  for  transverse  matrix  (slit) 
cracks  the  self-consistent  results  were  extensively  discussed  by  Dvorak,  Laws 
and  Hejazi  [2].  Subsequently  these  results  were  favorably  compared  by  Laws 
and  Dvorak  [3]  with  some  exact  results  of  Delameter,  Herrmann  and  Barnett  [4], 
see  Table  1.  Further,  Dvorak  and  Laws  [5]  showed  that  the  predictions  of  the 
self-consistent  model  were  in  good  agreement  with  the  experimental  data  of 
Highsmlth  and  Reifsnider  [6]  for  (0,903)3  E-glass  epoxy  laminates,  see  Fig.  1. 
In  addition  we  note  that  the  model  may  be  used  to  give  results  for  transverse 
cracking  of  (90j,0)g  laminates.  However,  Highsmlth  [7]  advises  that  the 
reported  (6]  data  for  (903,0)g  E-glass  epoxy  laminates  is  not  reliable. 

It  is  worth  noting  here  that  the  model  developed  by  the  present  authors 
(3,5]  is  entirely  consistent  [8]  with  bounds  obtained  by  Hashin  [9]  for  the 
loss  of  stiffness  of  cracked  laminates. 

At  the  time  of  writing,  a  popular  alternative  to  the  self-consistent 
model  is  the  differential  scheme.  We  remark  that  the  first  application  of  the 
differential  scheme  to  cracked  solids  was  given  by  Laws  and  Dvorak  (10].  As 
far  as  cracked  laminates  are  concerned,  it  is  straightforward  to  use  the  dif¬ 
ferential  scheme  rather  than  the  self-consistent  scheme  in  the  Dvorak-Laws 
model  for  cracked  laminates.  Indeed,  in  many  cases,  the  results  are  almost 


Identical. 

We  emphasize  that  the  Dvorak-Law  model  for  the  loss  of  stiffness  of 
cracked  laminates  Is  completely  deterministic  In  the  sense  that  the  model  does 
not  contain  any  adjustable  parameters. 

Next  we  turn  to  the  analysis  of  fiber  breaks  In  a  unidirectional  com¬ 
posite,  or  In  an  Individual  ply  of  a  laminate.  An  essential  prerequisite  to 
this  study  Is  the  analysis  of  crack  opening  displacements  by  Laws  [10].  It 
Is,  perhaps,  relevant  to  point  out  that  the  analysis  of  Laws  [10]  corrects  a 
rather  serious  error  In  the  literature.  The  paper  by  Laws  and  Dvorak  [11] 
breaks  new  ground  In  many  Important  areas.  In  particular  we  gave 

(1)  the  self-consistent  model  for  the  loss  of  stiffness  of  a  unidirec¬ 
tional  composite  due  to  fiber  breaks  accompanied  by  penny-shaped 
cracks  at  the  ends  of  the  broken  fibers, 

(11)  the  differential  scheme  for  the  same  problem, 

(111)  the  appropriate  Hashln-Shtrlkman  bound  for  cracked  solids, 

(Iv)  a  succinct  derivation  of  the  Mori-Tanaka  [12]  model  which  showed 
that  the  model  coincided  with  the  Hashln-Shtrlkman  bound, 

(v)  a  complete  numerical  comparison  of  the  various  models,  see  e.g. 
Figs.  2,  3, 

(vi)  the  calculation  of  various  energy  release  rates,  thus  Indicating 
the  likelihood  of  multiple  cracking  rather  than  the  catastrophic 
propagation  of  a  single  crack. 

We  note  that  the  results  described  briefly  above  have  been  used  exten¬ 
sively  by  Laws  and  Brockenbrough  [13,14,15],  and  others,  in  the  study  of 
mlcrocracked  polycrystalline  ceramics. 
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I 

I  2.2  First  Ply  Failure  of  Composite  Laminates 

I 

The  mechanics  of  crack  initiation  in  an  elastic  fibrous  ply  were  explained 
In  [16]  and  [17].  Cracks  were  assumed  to  Initiate  from  a  nucleus  created  by 
localized  fiber  debonding  and  matrix  cracking.  Conditions  for  the  onset  of 
unstable  cracking  from  such  nuclei  were  evaluated  with  due  regard  to  the 
Interaction  of  cracks  with  adjacent  plies  of  different  elastic  properties.  It 
was  found  that  crack  propagation  In  the  direction  of  the  fiber  axis  controls 
the  strength  of  thin  piles,  whereas  cracking  In  the  direction  perpendicular  to 
the  fiber  axis  determines  the  strength  of  thick  plies.  The  theory  relates  ply 
thickness,  crack  geometry  and  ply  toughness  to  ply  strength. 

In  order  to  apply  the  theory  presented  by  Dvorak  and  Laws  [16,17]  to  a 
specific  laminate,  It  Is  necessary  to  evaluate  the  reduction  In  the  stress 
Intensity  factor,  at  the  crack  tip  of  a  transverse  crack,  due  to  the  inter¬ 
action  with  the  adjacent  plies  in  the  laminate.  This  calculation  Is  especially 
difficult  and  Involves  the  solution  of  a  complicated  singular  Integral 
equation  -  the  details  are  to  be  reported  in  a  forthcoming  Ph.D.  thesis  [18]. 
When  use  Is  made  of  these  stress  Intensity  reduction  coefficients  together 
with  experimental  data  reported  by  Wang  [19],  one  obtains  the  comparisons 
shown  In  Fig.  4. 

2.3  Progressive  Cracking  of  Cross  Ply  Laminates 

A  further  significant  achievement  has  been  the  formulation  of  a  simple 
shear-lag  model  to  predict  progressive  transverse  cracking  of  cross  ply 
laminates  under  monotonlc  loading.  The  model  Is  based  upon  fracture  mechanics 
but  accounts  for  the  statistical  nature  of  the  transverse  cracking  process 


observed  in  experimental  tests.  The  theory  has  been  developed  in  references 


[8]  and  [20].  In  particular  we  suggested  a  definite  choice  of  probability 
density  function  to  account  for  the  statistical  variations  in  the  locations  of 
the  respective  cracks.  Thus  the  model  does  not  involve  any  adjustable  param¬ 
eters.  The  model  allows  us  to  predict  the  loss  of  stiffness,  progressive 
crack  density.  Initiation  of  H-cracking  at  Interfaces,  etc.  when  we  are  given 
basic  data  on  ply  thermo-elastic  properties,  ply  toughness  and  laminate 
geometry. 

Figure  6  shows  the  predicted  loss  of  stiffness  for  the  bench-mark  data  of 
Highsmlth  and  Relf snider  [6].  Also  shown  in  Figure  6  are  the  Hashln  [9] 
lower  bound  and  the  Highsmith-Reif snider  [6]  prediction.  We  observe  that  the 
Ovorak-Laws  self-consistent  result  is  not  shown  In  Fig.  6  since  it  is  almost 
Indistinguishable  from  the  Hashln  [9]  bound.  In  the  final  analysis  one  can 
conclude  that  predicting  the  loss  of  stiffness  is  quite  a  forgiving  process. 
Indeed,  when  one  considers  graphite-epoxy  systems,  the  loss  of  stiffness  due 
to  transverse  cracking  is  so  small  that  all  models  must  give  almost  Identical 
results.  However,  when  we  consider  the  prediction  of  progressive  crack 
density,  a  totally  different  picture  emerges. 

As  far  as  we  are  aware  the  only  other  work  which  addresses  the  problem  of 
determining  crack  density  under  monotonic  load  is  due  to  Wang,  Crossman  and 
co-workers  [19,21,22,23,24].  The  predictions  of  the  model  developed  by  the 
present  authors  [8,20]  are  extremely  close  to  the  observed  values,  see  Figs.  9 
and  10.  We  note  that  we  have  omitted  Wang's  [24]  numerical  predictions  from 
Fig.  9  since  it  is  impossible  to  do  justice  by  attempting  to  reproduce  the 
published  graphs.  Nevertheless,  it  is  significant  that  both  the  Wang-Crossman 
model  and  the  Laws-Dvorak  model  give  excellent  predictions. 

In  our  own  experimental  program  conducted  by  Dvorak  and  Martine  at  RPI, 
we  examined  glass-epoxy  tube  specimens  of  the  On/^®m/*^n  (n  •  1,  2;  m  •  1,  2,  4) 


layup  under  Incrementally  applied,  combined  tension  and  torsion  loads.  The 
crack  density  In  the  90“  piles  and  the  stiffness  losses  In  the  laminated  tubes 
were  measured  In  the  experiments.  In  agreement  wxth  the  predictions  of  our 
analysis  of  the  role  of  ply  thlckne<^s  In  damage  development,  we  found  no  pre- 
failure  damage  In  single-ply,  l.e.  ,  0/90/0  specimens.  The  single  ply  thickness 
In  our  specimens  was  0.254  mm.  Some  Incomplete  cracks  were  found  In  O/9O2/O 
tubes,  and  many  circumferential  cracks  were  observed  In  four-ply,  O/9O4/O 
laminates.  The  Implication  Is  that  the  amount  of  damage  In  laminates  can  be 
reduced  If  ply  thickness  Is  kept  to  a  minimum. 

The  above  experimental  Investigation  also  showed  that  transverse  cracking 
In  the  tested  system  was  caused  almost  exclusively  by  transverse  normal 
stresses  In  the  90“  layer.  Combined  loading  In  tension  and  shear  showed  that 
the  shear  stress  made  no  significant  contribution  to  transverse  cracking,  but 
that  It  was  responsible  for  delamination  between  plies  prior  to  final  failure. 

The  experiments  also  confirmed  predictions  of  stiffness  losses  by  the 
self-consistent  method.  The  agreement  was  very  good  in  the  case  of  axial 
elastic  modulus  of  the  tube,  but  less  satisfactory  in  the  case  of  the  shear 
modulus.  Viscous  effects  were  apparently  responsible. 


2.4  Fatigue  of  Metal  Matrix  Laminates 

The  experimental  results  of  Dvorak  and  Johnson  [26]  on  B/A£  laminates,  as 
well  as  subsequent  studies  by  Johnson  on  SiC/Ti  plates  suggest  that  cyclic 
plastic  straining  of  the  matrix  Is  the  principal  cause  of  fatigue  damage 
growth.  No  damage  is  typically  observed  if  the  laminate  Is  loaded  by  an 
elastic  load  cycle,  either  In  the  Initial  or  shakedown  state.  This  argument 
can  be  extended  to  damaged  laminates.  In  particular,  one  can  assume  that  all 
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damage  growth  will  terminate  If  the  laminate  reaches  an  elastic  deformation 
state.  As  Illustrated  by  the  example  In  Fig.  8,  elastic  deformation  can  be 
restored  under  an  Initially  Inelastic  cycle  of  loading*  If  the  amount  of 
damage  In  the  piles,  and  the  applied  strain  cycle  cause,  respectively,  expan¬ 
sion  and  translation  of  the  ply  relaxation  surfaces  such  that  the  prescribed 
load  cycle  or  the  corresponding  strain  cycle  can  be  accommodated  within  the 
new  elastic  range.  The  damage  evolution  process  can  then  be  regarded  as  a 
mechanism  that  the  con^)oslte  laminate  employs  to  reach  an  elastic  state.  In 
this  new  state,  the  originally  inelastic  part  of  the  total  strain  in  each  ply 
is  accommodated,  in  part,  by  the  strain  caused  by  opening  of  the  cracks. 

Viewed  from  a  different  perspective,  the  damage  process  can  be  thought  of 
as  a  part  of  a  shakedown  mechanism  in  the  composite  laminate.  According  to 
the  static  or  Melan  shakedown  theorem,  an  elastic-plastic  solid  or  structure 
will  shake  down  if  any  admissible  residual  stress  field  can  be  found  such  that 
its  superposition  with  the  stresses  caused  by  the  applied  loading  will  not 
violate  the  yield  condition  anywhere  In  the  solid.  In  other  words,  the  lamin¬ 
ate  will  shake  down  if  a  subsequent  yield  surface,  or  its  relaxation  surface 
counterpart  in  the  strain  space,  can  be  found  which  contains  the  applied  load 
or  strain  cycle.  Of  course,  shakedown  can  take  place  only  if  the  structure  is 
loaded  within  its  failure  envelope,  and  if  early  collapse  by  incremental  plas¬ 
tic  straining  can  be  prevented.  That  is  usually  the  case  in  laminated  plates 
where  the  elastic  fibers  support  a  major  part  of  the  load  so  that  the  total 
strains  are  small,  yet  substantially  larger  than  the  initial  yield  strains  of 
a  ply. 

Modeling  of  the  incremental  shakedown-damage  process  in  a  metal  matrix 
laminate  can  be  illustrated  with  the  help  of  relaxation  surfaces  shown  in 
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Fig.  8.  In  this  example,  the  yield  surfaces  in  strain  space,  so-called  relax¬ 
ation  surfaces,  were  plotted  for  a  (0/90)g  laminate  subjected  to  in-plane 

normal  strains  e  ,  e  .  The  top  bars  indicate  laminate  strains,  which  .are 
22  33 

equal  in  all  plies;  the  x^  direction  coincides  with  the  fiber  orientation  in 
the  0“  ply,  and  x^  is  the  fiber  direction  in  the  90®  ply.  Cracks  have  been 
added  after  cong)letlon  of  the  load  cycle,  and  the  relaxation  surfaces  are 
therefore  plotted  from  the  end  position  at  50  MPa.  The  dashed  line  indicates 
the  crack  opening  condition.  As  one  would  expect,  a  larger  overall  strain  is 
needed  to  cause  the  initial  yield  strain  in  the  matrix  of  a  cracked  ply. 
Therefore,  the  open  crack  branches  of  the  relaxation  surfaces  expand  with 
Increasing  crack  density  3. 

It  is  useful  to  point  out  that  cyclic  plastic  loading  of  the  laminate 
creates  cyclic  plastic  strains  in  Individual  plies  which,  as  illustrated  by 
Fig.  8,  tend  to  reach  a  steady  state  after  relatively  few  cycles.  On  the 
other  hand,  the  plastic  deformation  cycle  also  promotes  low  cycle  fatigue 
damage  growth  which,  in  comparison,  proceeds  very  slowly.  Typically,  several 
thousand  cycles  may  be  needed  to  cause  a  significant  change.  One  may  then 
expect  the  relaxation  surfaces  to  translate  much  more  rapidly  than  expand. 

The  direction  of  translation  should  be  such  as  to  minimize  the  magnitude  of 
plastic  work  per  cycle.  Under  such  circumstances,  the  relaxation  surfaces 
will  tend  to  translate  into  such  most  favorable  positions  which  will  assure 
that  the  amount  of  expansion  -  which  is  to  say  extent  of  ply  damage  -  will 
reach  only  the  minimum  amount  necessary  for  an  elastic  accommodation  of  the 
loading  program. 

Another  consequence  of  the  large  disparity  between  deformation  and  damage 
rates  is  that  in  each  damage  state  the  effect  of  past  deformation  history  will 


fade  very  quickly.  That  is  to  say  that  the  deformation  field  In  the  laminate 
at  a  particular  state  of  ply  damage  will  be  very  similar  to  that  which  one 
would  reach  If  this  amount  of  damage  was  Introduced  Into  an  elastic  laminate 
prior  to  the  application  of  the  corresponding  load  cycle.  The  Implication  Is 
that  the  final  deformation  state  In  a  damaged  laminate  subjected  to  a  constant 
load  cycle  must  be  Independent  of  previous  loading  and  damage  history.  There¬ 
fore,  In  arriving  at  a  final  damage  state,  one  may  follow  any  convenient  path. 
For  example,  the  damage  analysis  of  the  typical  case  of  constant  amplitude 
loading,  which  causes  large  excursions  Into  the  plastic  range  during  many 
Initial  deformation  cycles,  can  be  replaced  by  analysis  of  damage  caused  by  a 
load  cycle  which  expands  at  a  rate  comparable  to  that  of  damage  growth.  In 
this  particular  case  each  Increment  In  amplitude  Is  followed  after  few  cycles 
by  a  saturation  damage  Increment  which  restores  elastic  straining  within  a  new 
shakedown  sate. 

Apart  from  the  above  arguments,  available  results  obtained  In  tension- 
tension  fatigue  tests  of  B/Al  metal  matrix  composite  laminates  support  the  path 
Independence  concept.  For  example,  Johnson  found  that  a  saturation  damage 
state  at  a  certain  maximum  load  amplitude  could  be  reached  by  dissimilar  load¬ 
ing  sequences.  Another  support  for  this  concept  was  established  by  experiments 
which  showed  that  the  amount  of  damage  In  a  laminate  was  determined  primarily 
by  the  load  amplitude  and  not  by  variations  of  the  mean  stress. 

We  now  present  some  results  which  Illustrate  certain  aspects  of  damage 
development  in  B/Al  laminates  subjected  to  cyclic  tension  loading.  Our  objec¬ 
tive  Is  to  find,  for  several  different  load  amplitudes,  the  amount  of  damage 
In  each  ply  that  Is  needed  to  reach  a  shakedown  state  In  the  laminate.  Of 
course,  the  stiffness  loss  and  the  Internal  stress  distribution,  particularly 
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the  fiber  stress »  are  also  of  Interest.  One  specific  laminate  under  consider¬ 
ation,  of  0/90  layup,  was  already  discussed  in  connection  with  Fig.  8.  In  ad¬ 
dition,  a  similar  analysis  was  performed  for  a  0“  plate.  The  final  load  cycle 
we  wish  to  reach  In  both  laminates  is  from  S^in  “  50  MPa  to  Sjj^x  *  SOOmax  MPa. 
The  actual  path  we  follow  starts  with  cycling  of  the  laminate  to  a  steady 
state,  as  in  Fig.  8.  Then,  while  is  kept  constant,  S^ax  reduced  to 

bring  the  laminate  into  an  elastic  state.  Next,  Sy^ax  Increased  in  small 

increments.  After  each  Increment,  cracks  are  Introduced  in  the  plastically 
deforming  plies  to  the  extent  needed  to  accommodate  the  deformation  path 
within  the  expanded  relaxation  surfaces.  This  Involves  both  expansion  and 
translation  of  the  surfaces,  which  cause  a  new  plastic  strain  state  and  a 
stress  redistribution  through  the  laminate.  Details  of  the  procedure  have 
been  described  by  Wung  (27]  and  Dvorak  and  Wung  [28].  As  an  example.  Fig.  9 
shows  the  current  relaxation  surfaces  at  two  levels  of  Sq,3x*  Note  that  nigh 
values  of  crack  density  8  are  needed  here  to  accomplish  the  accommodation.  It 
would  be  unrealistic  to  expect  crack  densities  exceeding  8  »  1  in  each  ply. 
Therefore,  8  must  be  regarded  here  as  a  damage  parameter  which  accounts  for 
extension  of  the  cracks  onto  adjacent  plies,  for  crack  intersection  at  ply 
boundaries,  and  for  crack  extension  by  delamlnatlon  along  fiber-matrix  Inter¬ 
faces.  The  incremental  expansion  of  the  loading  range  continues  until  one 
reaches  the  desired  final  magnitude  of 

Fig.  10  shows  the  change  In  the  axial  elastic  modulus  caused  by  satura¬ 
tion  damage  in  the  two  laminates  as  a  function  of  the  applied  tension  stress 
range.  The  computed  results  are  plotted  together  with  experimental  data  of 
Dvorak  and  Johnson  [26].  In  the  experiments,  the  saturation  damage  state  was 
defined  as  the  damage  scate  after  2x10^  cycles  at  constant  stress  amplitude. 
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as  noted  In  Fig.  10,  but  actual  measurements  of  stiffness  loss  Indicated  that 
damage  usually  stabilized  after  5x10^  cycles.  Finally,  Fig.  11  Indicates  the 
computed  magnitudes  of  the  axial  stress  In  the  0°  layer  fibers.  In  the  satura¬ 
tion  state  at  different  levels  of  This  stress  change  has  been  plotted 

up  equal  to  the  experimentally  observed  endurance  limit.  Note  that  while 

the  endurance  are  quite  different  In  the  two  laminates,  the  terminal 

fiber  stresses  are  nearly  Identical.  The  Implication  Is  that  fatigue  failure 
occurs  In  these  composite  systems  by  overloading  of  the  0°  fibers.  Of  course, 
the  maximum  stress  Is  not  seen  by  the  fibers  until  the  laminate  reaches  the 
saturation  damage  state  at  the  endurance  In  the  Initial  stages  of 

damage  development,  a  part  of  the  load  Is  carried  by  the  undamaged  off-axis 
piles,  but  as  damage  grows  more  stress  Is  absorbed  by  the  0°  fibers. 
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3.  WORK  IN  PROGRESS 


Two  Ph.D.  dissertations  are  still  in  progress,  one  by  Mr.  A.  Kaveh 
Ahangar  [29]  at  RPI,  and  one  by  Mr.  J.B.  Wang  [18]  at  the  University  of 
Pittsburgh.  This  work  in  progress  can  be  described  as  follows. 

Mr.  Ahangar 's  dissertation  is  concerned  with  evaluation  of  overall  stiff¬ 
ness  changes,  and  the  local  fields  in  plies,  for  cracked  laminates  of  any  lay¬ 
up.  The  approach  is  based  on  variational  principles  of  elasticity,  which 
allow  one  to  select,  among  a  class  of  chosen  admissible  fields  in  the  plies, 
the  best  approximation  of  the  actual  stress  field  in  a  cracked  ply.  The 
chosen  fields  do  not  take  into  account  crack-tip  singularities,  this  seems  to 
be  reasonable  in  fibrous  systems  where  the  crack  tip  geometry  is  not  well  de¬ 
fined.  However,  interaction  between  cracks  in  a  ply,  and  of  the  cracks  with 
adjacent  plies  is  taken  into  consideration*  This  technique  leads  to  much 
better  estimates  of  ply  fields  than  the  self-consistent  method  which  gives 
only  ply  stress  averages.  We  have  already  completed  the  analysis  of  laminates 
with  specified  crack  densities.  Current  work  addresses  the  problem  of  crack 
growth  under  Incrementally  applied  load.  Results  of  this  kind  have  been 
available  so  far  only  for  0/90  laminates.  The  method  relies  on  a  numerical 
procedure,  but  the  conqputlng  effort  is  relatively  moderate. 

Work  is  also  continuing  on  the  exact  stress  analysis  of  cracked  laminates. 
To  a  certain  extent  this  will  be  reported  in  the  forthcoming  Ph.D.  thesis  by 
Wang  [18].  But  we  continue  to  strive  for  a  computer  code  which  is  sufficiently 
user-friendly  and  cost  effective  to  analyse  laminates  of  arbitrary  lay-up. 

This  is  not  easy  and  Indeed  is  impossible  without  use  of  the  supercomputer. 

As  remarked  earlier  in  this  report,  this  program  provides  the  stress  Intensity 
factors  which  are  essential  in  the  Dvorak-Laws  model  for  first  ply  failure. 
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The  worlc  on  progressive  cracking  Is  continuing  and  problems  associated 
with  other  lay-ups,  e.g. ,  (±25,90)g,  (  (0,90)ii)8»  ^re  being  addressed.  In 
addition  we  are  continuing  with  our  work  on  the  onset  of  additional  damage 
modes  together  with  the  progressive  nature  of  such  modes.  A  particular  prob¬ 
lem  of  concern  here  relates  to  the  development  of  damage  during  bending. 

But  the  major  part  of  the  work  In  progress  Is  now  related  to  ceramic 
matrix  cooqposltes  and  laminates.  Here  we  are  Investigating  the  change  of 
thermomechanical  properties  of  both  unidirectional  composites  and  laminates. 
Typical  problems  of  Interest  are  the  loss  of  stiffness  and  Improvement  of 
toughness  of  various  ceramic  matrix  ceramic  fiber  systems.  It  Is  clear  from 
the  literature  and  from  discussions  with  Or.  Nicholas  Pagano  and  with  Or.  Ted 
Nicholas  at  the  Air  Force  Materials  Laboratory  that  unidirectional  materials 
pose  significantly  different  Issues  to  those  which  must  be  addressed  for 
laminates.  Further,  It  Is  also  clear  that  the  damage  modes  in  ceramic  matrix 
ceramic  fiber  systems  (unidirectional  or  laminates)  are  different  from  the 


m 


modes  observed  In  graphite-epoxy  or  glass-epoxy  systems.  These  problems  are 
under  Intensive  study. 
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Experimental  data  from  reference  [6]. 

Figure  6:  Theory  versus  experiment  for  progressive  cracking  of  AS-3501-06 
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Figure  7:  Theory  versus  experiment  for  progressive  cracking  of  T300/934 

laminates.  Data  from  Wang  {24]. 
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elastic  modulus  of  two  B-A£  laminates.  Comparison  of  prediction 
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Figure  11:  Fiber  stresses  in  zero-degree  plies  after  damage-induced  shakedown 
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ABSTRACT 

Overall  stiffnesses  and  compliances,  thermal  expansion  coefficients,  and  stress  and 
strain  averages  are  evaluated  for  a  fibrous  composite  lamina  which  contains  a  given 
density  of  open  transverse  cracks  and  is  subjected  to  uniform  mechanical  loads  and 
thermal  changes.  The  evaluation  procedure  is  based  on  the  self-consistent  method  and 
is  similar,  in  principle,  to  that  used  in  finding  elastic  constants  of  unidirectional  com¬ 
posites. 

1.  INTRODUCTION 

IN  MANY  FIBROUS  COMPOSITE  SYSTEMS  THE  FAILURE  STRAIN  OF  THE  FIBER 
far  exceeds  that  of  the  matrix.  Under  load,  the  difference  is  usually  accom¬ 
modated  by  matrix  cracking.  This  is  frequently  observed  in  monotonically  or 
cyclically  loaded  laminated  plates,  where  each  layer  may  contain  a  system  of 
aligned  slit  cracks  which  grow  in  the  direction  of  the  fibers  and  across  the 
thickness  of  the  ply.  Such  cracks  are  often  called  ply  cracks  or  transverse 
cracks,  although  it  is  more  appropriate  to  reserve  the  latter  for  cracks  which 
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are  perpendicular  to  the  fiber  axis,  and  refer  to  cracks  which  grow  parallel  to 
the  fibff  direction  as  axial  cracks. 

In  polymer  matrix  composites  axial  matrix  cracking  typically  stans  at  low 
strain  levels  in  the  weakest  off-axis  ply  (1-31.  As  loading  continues,  cracks  ap¬ 
pear  in  other  off-axis  plies;  also,  their  density  increases  until  it  reaches  a  cer¬ 
tain  saturation  level.  For  example,  in  statically  loaded  (0/90),  graphite-epoxy 
laminates  the  minimum  crack  spacing  was  observed  to  be  equal  to  3.S-4.0  ply 
thicknesses  [4].  In  metal  matrix  composites,  matrix  cracking  appears  to  be 
caused  only  by  cyclic  loads  which  exceed  the  shakedown  limit  of  the  laminate 
[5].  Under  such  circumstances  the  matrix  experiences  cyclic  plastic  straining 
and,  consequently,  low-cycle  fatigue  failure.  Both  axial  and  transverse  crack¬ 
ing  is  present,  the  former  in  off-axis  plies,  the  latter  in  zero  degree  plies.  The 
crack  patterns  and  densities  are  generally  similar  to  those  found  in  polymer 
matrix  systems.  However,  saturation  density  increases  with  load  amplitude, 
and  it  is  not  unusual  to  find  cracks  as  close  as  one  ply  thickness. 

In  a  typical  part  of  a  laminated  composite  structures,  removed  from  con¬ 
centrated  loads  and  free  edges,  matrix  cracking  is  the  initial,  low-stress 
damage  mode  under  applied  load,  it  is  eventually  followed  by  other  types  of 
damage,  such  as  cracking  between  layers  and  fiber  breaks;  but  these  appear  at 
relatively  high  loads  which  may  exce<Kl  allowable  design  magnitudes.  In  con¬ 
trast,  matrix  cracks  grow  at  low  loads,  and  they  can  significantly  impair  stiff¬ 
ness  and  strength  of  composite  structures,  especially  those  containing  many 
off-axis  plies.  For  example,  fatigue  tests  on  toth  polymer  and  metal  matrix 
laminated  plates  indicate  that  stiffness  and  residual  static  strength  reductions 
caused  by  cracking  in  plies  may  be  as  high  as  10-50%  after  2  x  10*  cycles  of 
loading  [5,6].  It  is  therefore  desirable  to  consider  the  effect  of  matrix  cracking 
on  composite  properties  in  design. 

Sufficiently  general  theoretical  models  of  progressive  cracking  in  com¬ 
posite  laminates  are  apparently  not  to  be  found  in  the  literature.  Such  results 
as  are  available  for  angle-ply  laminates  have  been  obtained  from  finite  ele¬ 
ment  calculations  [4,7],  while  other  studies  have  focused  on  (0/90„), 
laminates  [2,3]. 

The  purpose  of  this  group  of  papers  is  to  introduce  a  modelling  procedure 
for  axial  cracking  which  can  be  applied  to  any  laminate  geometry  under 
general  loading  conditions. 

The  first  step  in  the  analysis,  undertaken  in  this  paper,  is  an  evaluation  of 
overall  thermoelastic  properties  of  a  fibrous  composite  which  contains  a  cer¬ 
tain  density  of  aligned  slit  cracks.  Figure  la.  In  subsequent  papers  we  intro¬ 
duce  appropriate  fracture  criteria  for  a  unidirection^  ply,  and  develop  a 
technique  for  evaluation  of  instantaneous  crack  densities  and  stiffnesses  of  a 
single  ply  which  is  strained  in  a  prescribed  way.  Finally,  these  results  will  be 
incorporated  into  an  analysis  of  laminated  plates  which  are  loaded  by  in¬ 
plane  stresses  and  contain  many  cracking  layers  of  different  orientation. 
Although  the  subsequent  parts  utilize  the  results  obtained  herein,  the  pro¬ 
cedures  developed  there  do  not  depend  on  the  particular  method  used  to 
evaluate  the  effect  of  cracks  on  thermoelastic  properties  of  a  lamina. 
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Figum  1. 1  A)  an  infinite  fibrous  medium  with  aligned  sHt  cracks:  (B)  a  fiber  lamina  with  parataf  slit 
cracks. 


In  this  paper  we  are  concerned  with  evaluation  of  overall  compliance, 
thermal  expansion  coefficients  (or  thermal  strain  and  stress  vectors),  and 
strain  and  stress  averages  in  the  phases  of  a  fibrous  lamina  which  contains 
aligned  slit  cracks,  and  is  subjected  to  uniform  mechanical  loads  and  thermal 
changes.  The  approach  to  the  problem  was  outlined  in  our  recent  study  [8], 
where  we  suggested  that  the  effect  of  matrix  crack  systems  on  properties  of 
fibrous  composites  could  be  analyzed,  in  principle,  by  the  techniques  which 
are  commonly  used  to  evaluate  the  elastic  constants  of  composite  materials 
and  fibrous  laminates,  e.g.,  by  the  self-consistent  method. 

The  essential  approximation  in  the  evaluation  of  stiffness  and  compliance 
changes  of  laminates  consists  in  the  replacement  of  a  cracked  layer.  Figure 
1  b,  by  an  effective  medium  which  contains  many  cracks.  Figure  la.  The  crack 
densities  can  be  exactly  matched  to  provide  identical  stiffnesses.  However,  the 
cracks  in  the  layer  are  not  entirely  surrounded  by  the  layer  material,  instead 
they  interact  with  neighboring  layers  which  have  different  elastic  properties. 
This  interaction  is  limited  to  the  vicinity  of  the  crack  tip  (9),  thus  it  may  be 
important  in  analysis  of  local  crack  growth  at  the  interface,  but  it  has  only  a 
minor  effect  on  lamina  stiffness.  We  note  that  a  similar  approximation  is 
commonly  accepted  in  evaluation  of  elastic  moduli  analysis  of  laminated 
composite  structures  reinforced  by  monolayers  of  large  diameter  fibers. 


2.  CRACK  DENSITY 

It  was  emphasized  in  our  earlier  work  (8)  that  interest  in  cracked  infinite 
unidirectional  composites  is  largely  determined  by  the  desire  to  study  com¬ 
posite  laminates.  Certainly  in  papers  of  this  series  we  assume  that  the  proper¬ 
ties  of  a  cracked  ply  are  adequately  predicted  by  the  infinite  model.  It  is, 
therefore,  expedient  to  interpret  the  analysis  in  this  paper  with  particular 
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reference  to  cracked  laminates.  Thus  the  terms  cracked  ply  and  cracked  in¬ 
finite  unidirectional  composite  are  here  synonymous. 

Consider  a  unidirectional  fibrous  medium  with  aligned  slit  cracks  of  length 
2a,  Figure  la.  By  way  of  contrast,  there  are  several  ways  of  viewing  the  com¬ 
posite  on  the  microscale.  In  the  first  place  one  could  distinguish  between  the 
three  coexistent  phases;  fiber,  matrix,  and  cracks.  On  the  other  hand,  one 
might  consider  the  uncracked  fibrous  composite  as  an  effective  homogeneous 
solid  into  which  cracks  have  been  introduced.  In  {8)  we  referred  to  the  two 
possibilities  as  the  three  phase  model  and  two  phase  model,  respectively.  The 
distinction  reflects  the  fact  that  in  certain  laminated  composite  systems  the 
crack  length  and  fiber  diameter  are  of  similar  magnitude,  while  in  other 
systems  the  fiber  diameter  is  much  smaller  than  the  crack  length.  A  signifi¬ 
cant  conclusion  which  may  be  drawn  from  [81  is  that  for  many,  purposes  the 
predictions  of  the  two  models  are  indistinguishable.  In  the  present  paper  we 
concentrate  on  the  simpler  two  phase  model. 

In  analysis  of  cracked  materials,  it  is  well  known  that  cracks  can  be 
modeled  either  as  limiting  ellipsoids  or  as  slits— the  final  results  are 
equivalent.  It  was  convenient  in  our  earlier  work  [8]  to  regard  cracks  as 
limiting  ellipsoids  and  we  adopt  the  same  approach  here.  On  the  other  hand 
we  note  that  in  a  study  of  crack  geometries  different  from  ours,  Gottesman, 
Hashin  and  Brull  (10]  use  slits  a  priori.  Thus  we  can  consider  the  cracks  in 
Figure  1  to  be  elliptic  cylindrical  cavities.  Let  a,  b  respectively  denote  the 
major  and  minor  semi-axes  of  a  cavity,  and  let  n  be  the  number  of  cavities  per 
unit  area  of  the  x,Xi-plane.  Then  the  volume  fraction  of  cavities  is  1/4  n^d 
where  6  <■  b/a  is  the  aspect  ratio  and  fi  »  4r)a*  is  the  crack  density  parameter. 
For  simplicity,  suppose  that  the  cracks  have  length  2a,  then  p  is  the  average 
number  of  cracks  in  a  square  of  side  2a.  Alternatively,  when  the  cracks  are 
located  in  a  ply.  Figure  lb,  p  measures  the  average  distance  between  cracks. 
Clearly,  if  the  cracks  are  equally  spaced  in  a  ply,  the  distance  between  suc¬ 
cessive  cracks  is  2a/|3.  Thus  when  (3  °  1  the  distance  between  cracks  is  equal 
to  2a.  As  p  decreases  the  distance  between  cracks  increases  and  when  p  *0 
there  are  no  cracks.  We  note  that  in  some  experiments  [4]  the  observed 
minimum  separation  of  cracks  (in  a  saturation  state)  in  Gr-Ep  laminates  was 
3.S  to  4.0  ply  thicknesses.  The  corresponding  values  of  p  are  0.28  and  0.25 
respectively.  These  values  may  be  contrasted  with  B-AJ  composites  where 
values  as  high  as  /J  =  1  have  been  observed  (5].  Therefore  we  present  all 
results  for  0  <  <  1 . 

As  cracks  are  generally  initiated  at  preexisting  flaws  which  are  randomly 
distributed  in  the  ply,  the  crack  distribution  at  any  given  moment  is  not 
regular.  However,  we  assume  that  it  is  statistically  homogeneous,  and  that  p 
is  a  measure  of  average  crack  density  in  the  in-plane  area  of  the  ply.  Also, 
since  cracks  sometimes  tend  to  propagate  rapidly  across  the  entire  loaded 
width  of  a  ply,  the  crack  density  does  not  change  continuously  but  may  in¬ 
crease  in  discrete  steps.  These  steps  can  be  measured  in  terms  of  a  change  LP 
caused  by  a  single  crack  in  a  utut  plane  area  of  a  ply. 

To  illustrate  the  relationship  between  the  crack  density  P,  and  the  number 
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Table  1.  Cracks  in  a  1/8  mm  thick  piy  12a  >  0.012S  cm). 


Crack  danaity  fi 

0.1 

0.25 

0.5 

0.75 

1.0 

Crack  spacing,  2a/ P  mm 

1.25 

0.50 

0.25 

0.17 

0.125 

Numbar  of  cracks/ Unit  square  plana 

area  of  pt^ 

No.  cracks/cm' 

8 

20 

40 

60 

80 

No.  cracks/in' 

20 

SO 

100 

ISO 

203 

No.  cracks/ 

244 

610 

1220 

1830 

2438 

of  cracks  in  a  plane  area,  we  focus  for  a  moment  at  the  XiXi-plane  in  Figure 
1  b.  The  trace  of  each  crack  in  this  plane  is  a  straight  line  parallel  to  Xj.  Table  1 
presents  some  typical  numbers  calculated  for  a  1/8  mm  thick  ply.  It  shows  the 
distance  between  cracks  2a//},  and  the  number  of  cracks  in  a  unit  area  (XjXj- 
plane  in  Figure  lb)  of  the  ply,  at  specified  magnitudes  of  /).  Note  that  even  at 
low  /)  values,  and  in  small  areas,  the  crack  counts  are  relatively  large. 

These  results  also  indicate  the  magnitude  of  an  increment  i/9  caused  by  a 
single  crack  in  a  unit  area  of  the  ply,  which  does  not  depend  on  /}.  For  exam¬ 
ple,  even  in  a  small  area  of  1  cm^  a  single  crack  changes  ft  by  0.012S.  This  in¬ 
crement  is  relatively  small,  so  that  p  may  be  regarded  as  continuous  in  most 
applications.  A  notable  exception  occurs  at  stress  concentrations  where 
stresses  may  change  significantly  within  a  small  area,  e.g.,  1  mm'.  A  single 
crack  will  cause  hfl  =  0.125/crack  in  a  1  mm'  area,  hence  p  may  not  be 
regarded  as  continuously  changing.  Moreover,  the  crack  count  becomes  small 
(~1  crack/mm  at  /)  =  0.1),  so  that  the  concept  of  average  crack  density  may 
no  longer  be  applied.  It  follows  that  the  averaging  approach  is  admissible  in 
evaluation  of  stiffness  changes  caused  by  cracks  in  large  in-plane  areas  of  a 
ply  or  laminate,  but  that  it  may  not  be  suitable  for  analysis  of  local  properties 
at  small  geometrical  imperfeaions  such  as  notches  or  free  edges. 

It  is  also  appropriate  to  discuss  the  definition  of  crack  density  in  exterior  or 
surface  plies  of  a  laminate.  As  in  Figure  I  b,  such  cracks  are  open  at  the  sur¬ 
face,  hence  their  effective  length  is  not  equal  to  the  ply  thickness.  However, 
the  surface  layer  may  be  regarded  as  half  of  a  layer  of  thickness  4a,  Figure  2, 
so  that  open  surface  cracks  may  be  approximated  by  interior  cracks  of  length 
4a.  The  effective  crack  density,  P‘,  in  exterior  layers  is  easily  obtained  in  terms 
of  the  crack  density,  p,  defined  for  interior  layers.  In  fact,  from  Figure  2 


2a/p  =  4a/P', 


P'  =  2P. 
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Rgura  2.  Evaluation  of  crack  density  parameter  ji'  in  exterior  layers  of  a  laminate. 


3.  EVALUATION  OF  OVERALL  COMPLIANCE 

On  the  macroscale,  the  cracked  unidirectional  composite  of  Figure  1  can  be 
regarded  as  an  orthotropic  homogeneous  solid.  The  elastic  properties  of  the 
"matrix"  are  identical  with  those  of  the  Tibrous  composite  and  can  be  easily 
evaluated.  When  cracks  are  introduced,  the  macroscopic  or  overall  elastic 
moduli  of  the  solid  change.  To  make  the  concept  of  overall  moduli  mean¬ 
ingful,  it  is  necessary  to  consider  overall  uniform  loading.  Thus,  we  introduce 
uniform  overall  average  stresses  o  and  strains  |,  with  components  arranged  in 
(6^1)  column  vectors  and  related  by  constitutive  equations* 


5  =  I|,  |  =  A/s.  (1) 

where  L ,  and  M  are  the  overall  stiffness,  and  compliance  (6  ^  6)  matrices  of 
the  cracked  composite,  respectively.  M  ^  L'  when  the  inverse  exists. 

Since  we  are  here  concerned  only  with  a  2-phase  model  (in  which  the  un¬ 
cracked  composite  is  the  matrix)  it  is  convenient  to  adopt  a  minor  variant  of 
the  notation  of  our  earlier  work  [8];  we  use  the  suffix  o  to  denote  the  proper¬ 
ties  of  the  uncracked  composite.  For  example,  L,  is  the  stiffness  of  the  un¬ 
cracked  composite. 

According  to  the  derivation  given  in  [8],  the  self-consistent  estimates  for 
the  overall  stiffness  and  compliance  matrices  are  given  by 


L  =  L.-PU  M 
A?  =  A/,  +  ^  A, 


(2) 


where 


•As  in  j8).  (6  "  6)  matrices  arc  denoted  capital  Latin  or  Greek  lcticr>.  L.  \t.  \  0.  and  tfr  «  II  matrices 

by  lower  case  letters  underlined  bv  a  tilda.  e.g..  o.  c.  m.L 


Wi^JWW^VWJWS^ 
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?  =  ^  n/J. 


(3) 


The  matrix  A  has  only  three  nonzero  components,  which  are  expressed  in 
terms  of  compliances  A/y  of  the  effective  medium  as: 

+  (4) 

Mii 

A44  =  fA/44A/5V>’"  (5) 

Mii 

where  a,  and  aj  are  the  roots  of 

fMjjA/jj  -  M\i)a^  -  —  A/,jA/u>}a  +  M,,A/j3  -  Af?3  =  0. 

(7) 

These  results  imply  that  only  three  compliance  coefficients  A/j,,  A/«4,  and  A/^ 
are  affected  by  the  introduction  of  cracks,  the  remaining  six  terms  in  M  are 
unchanged,  i.e.,  they  remain  equal  to  those  of  the  uncracked  fiber  composite. 
In  particular 

A/|i  ~  A^i  ,  A/i2  =  A/?2  ,  A/|j  —  , 

Miy  —  Mli  ,  Mjs  ~  A/Sj  >  A/js  =  A/j5  , 

a/,2  =  a/?,  +  |3fA/22A/„  -  A^2/('a,‘'^  +  <ri'V/A/„  .  (8) 

A/«  =  An.  +  ,  (9) 

A/«  =  A/S»  +  P(MllMyy  -  -  M]y)''^(aV^  +  ai'^J/Myy;  (10) 

and  from  (7); 

®1®,  ~  (A/nA/33  ~  ~  M\y) 

a.  +  o,  =  [MiiMtt  +  ~  Mt^diy)]/lMxiM3i  ~  M\yJ. 

The  unknown  shear  compliance  A/m  can  be  obtained  directly  from  (9): 

A/m  =  An.  +  ^  A/,5  +  (P  Mli  +  /A/?4  A/,s> -'7  (11) 


The  remaining  unknowns  A/22  and  A/„  are  found  from  (8)  and  ( 10).  The  solu- 
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tion  of  these  equations  is  perhaps  best  explained  with  the  help  of  additional 
notation.  Thus  let 

X  =  -  Mii  ,  Jt.  =  MliMu  ~  Mii  (12) 

y  =  MJvl„  =  (13) 

p  =  -M\,  .  O'  =  2(MnM»  -  (14) 

Now,  it  is  easy  to  show  that 

+  (p/x)''^(x  -xj 

and  that  x  is  the  solution  of 

F(x)  sx-x,-  Pl3p''^  X  ■>r  (y,->r  q)  X-  p'’^  x,  x  =  0.  (15) 

At  this  stage  it  is  not  dear  that  the  self-consistent  estimate  for  x  is  unique. 
This  apparent  deficiency  is  easily  remedied  by  noting  that  the  square 
bracketed  terms  in  (15)  is  always  positive.  Thus  x'^  x,.  Next  one  could  for 
example  rewrite  (15)  as  a  quartic  in  x''^  and  use  the  Descartes’  rule  of  signs  to 
show  that  Equation  (15)  has  one  and  only  one  root  x>  x^. 

Results  obtained  for  several  composite  systems  indicate  that  an  approx¬ 
imate  value  of  the  positive  root  x  =  x,  can  be  found  as: 


~  FVxJ 


(16) 


This  is  illustrated  in  Figure  3,  where  FM  plotted  against  x  is  almost  a  straight 
line  between  x,  and  Jt, . 

Once  X,  has  been  found,  Afn  and  A/m  are  obtained  from  (12)  and  (13), 
respectively. 

The  preceding  argument  may  suggest  that  it  is  better  to  solve  equation  (2i) 
for  the  compliance  matrix  then  equation  (2,)  for  the  stiffness  matrix.  This  is 
not  the  case.  The  analysis  was  presented  in  terms  of  the  compliance  to  furnish 
a  simple  proof  of  the  uniqueness  of  the  self-consistent  solution.  From  a  com¬ 
putational  standpoint,  it  is  just  as  easy  to  solve  (2,)  for  the  stiffness  as  it  is  to 
solve  (2i)  for  compliance.  Once  either  is  known  the  other  is  found  by 
inversion. 

If  the  A/4,  are  known,  then  one  can  find  the  three  nonzero  components  of 

in  (6)  from  (4).  The  overall  stiffness  L  can  then  be  written  in  terms  of 
the  expressions  appear  in  (A-6)  in  Appendix. 

We  emphasize  that  the  results  obtained  above  are  valid  only  for  the  case  of 
open  cracks  (c.f.,  (32)  below).  If  cracks  are  closed,  the  overall  response  may 
be  approximated  by  that  of  the  uncracked  composite. 
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x(10'  «Pa)‘2 

Figure  3.  Solution  of  Equation  1 15). 


4,  RESPONSE  TO  THERMAL  CHANGE 

Let  0  denote  a  uniform  temperature  change  from  a  reference  temperature 
0,,  applied  to  the  composite.  The  constitutive  Equations  (1)  must  be  replaced 
by 

a- LI- 91  (17) 

where  m  is  the  thermal  expansion  vector  and  I  the  thermal  stress  vector  (11]. 
From  (17)  one  obtains 

l^Lm.  (18) 

It  is  probably  obvious  that  the  presence  of  cracks  does  not  affect  free  ther¬ 
mal  expansion  of  the  composite.  Thus 


and  from  (2)  and  (18): 


m  =  m„ , 


I  =  (/  -  '(iL\)L 


Explicit  forms  of  m  and  i  are: 


m 

io 


P 

P 


m 


m. 


V' 


Analysis  of  Progressive  Matrix  Cracking  in  Composite  Laminates 


225 


Ot 

fLu  +  LtiJOf  +  Ltj  ot^ 

ar 

(Lxi  Lxi)Qj  +  Lix 

Oa 

;  1  = 

0 

0 

0 

0 

_0_ 

_  0  _ 

(21) 


where  or  and  are  linear  thermal  expansion  coefficients  of  the  uncracked 
composite  in  the  transverse  and  axial  directions,  respectively.  Since  the  coeffi¬ 
cients  of  change  with  /3,  f  is  also  a  function  of  ft,  while  m  remains  con¬ 
stant. 

Again,  the  derived  expressions  are  valid  only  if  the  cracks  are  open.  Closed 
cracks  do  not  affect  thermal  response  of  the  composite. 


5.  LOCAL  STRAIN  AND  STRESS  AVERAGES 

Certain  applications  of  the  above  results  require  information  about  local 
stress  and  strain  fields  in  the  phases.  Estimates  of  local  fields  are  also  needed 
to  distinguish  between  loading  conditions  which  lead  to  either  open  or  closed 
cracks  and  thus  delineate  limits  of  validity  of  the  theory. 

It  is  clear  that  when  a  fiber  reinforced  composite  containing  cavities  is  sub¬ 
jected  to  general  uniform  mechanical  loading,  part  of  the  applied  strain  (i.e., 
the  average  strain)  is  accommodated  by  cavity  expansion.  Whereas  the  strain 
of  a  typical  cavity  wall  becomes  unbounded  in  the  limit  of  vanishing  aspect 
ratio,  the  important  quantity  in  the  macroscopic  study  of  composites  is  the 
wall  strain  multiplied  by  the  area  of  the  cross-section  (or  aspect  ratio)  which 
tends  to  a  finite  limit.  Following  Laws  (12),  let  c'  be  the  strain  of  a  typical 
cavity  wall,  and  set 


^  =  Urn  6  t'  .  (22) 

•'  d  -*  o' 

Then  it  is  easy  to  show  that  (c.f.,  [12]) 

^  =  A  0  .  (23) 

Furthermore,  from  the  formulae  given  in  [8],  and  in  Section  4  above,  it  is 
not  difficult  to  show  that  in  the  limit  of  slit  cracks  (d  -*  0),  the  overall  strain, 
I  is  given  by 


I  ~f»  ■'■I'  ■ 


(24) 


where  ij,  is  the  average  strain  in  the  composite  “matrix": 


ww?toct*ctvwa: 
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^  =  f/-pAL}c  +  l)ALm6.  (25) 

In  addition,  ^  is  the  crack  accommodation  strain 

i,  =  p  A  Ui  -  mQ)  ;  (26) 

it  represents  that  part  of  the  applied  (average)  strain  which  is  accommodated 
by  cracks. 

In  practice  Equations  (25)  and  (26)  provide  us  with  estimates  of  the  local 
stress  and  strain  fields  within  the  composite. 

Viewed  from  another  standpoint.  Equations  (25)  and  (26)  immediately 
furnish  the  strain  concentration  factors.  As  far  as  the  composite  “matrix”  is 
concerned. 


^  =  A  I  -  da  (27) 

where 

A  =  I -p  AL 

(28) 

a  - (A -l)m--pALm 


It  is  perhaps  important  to  emphasize  the  different  physical  interpretations 
of  i,  and  i, .  Whereas  ^  is  the  average  strain  in  the  matrix,  e,  represents  that 
part  of  the  overall  strain  which  is  accommodated  by  all  the  cracks.  Therefore, 
one  can  write  (26)  as 


|r  =  -4,  £  -  , 

where 

A^^^PAL 
a^==p  A  L  m  . 


(29) 


(30) 


For  example,  it  is  immediately  obvious  that  the  average  thermal  strains  in  a 
fully  constrained  cracked  composite  (c  =  0)  are 

u  ~  ~  I,  =  P  A  L  m  B  . 

The  matrix  expands  at  the  expense  of  crack  closing. 

For  computational  purposes  it  is  advantageous  to  review  (28,)  with  the  help 
of  (2)  in  the  form 
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It  is  often  convenient  to  have  explicit  forms  of  the  above  results.  The 
nonzero  components  of  A  in  (28)  are; 

A,  I  —  Afj  —  Aff  =  1 

Alt  “  “  ^ii)  Lti 

All  ~  (Mil  ~  Mu)  Lii  +  1 

(32) 

An  =  (Mil  ~  Mu)  Lu 

An  —  Mn 

Afytt  —  M^tt^Mtttt . 


Again,  A/J  are  compliances  of  the  uncracked  composite;  A/y,  L„  are  listed  in 
(A-3)  to  (A-6)  in  the  Appendix. 

The  components  of  thermal  strain  concentration  factor  q  in  (28)  are  all 
zero,  except  for 


Qu  ~  p  ^u{(i-ll  Lu)Ot  OaJ  (33) 


The  average  strain  ^  (26).  (29)  in  the  cracks  has  the  following  nonvanishing 
components: 

Ell  ~  P  Ajj/I.ij(7ii  ~  ajP)  +  Lu(Tii  ~  ard)  +  LufCji  ~  o^d)J 

(34) 

ejj  ^  (1  ~  Mm/M^His  ell  ~  (I  ~  AiS*/A/66)£ii , 

We  note  that  the  distinction  between  open  and  closed  cracks  is  determined 
by  the  inequality 

eU>0,  (35) 

which  can  be  evaluated  with  the  help  of  (34). 

The  results  presented  here  are  valid  only  if  this  inequality  is  satisfied.  If  the 
sign  is  reversed,  or  if  the  two  sides  are  equal,  then  cracks  are  closed  and  do 
not  affect  the  mechanical  and  thermal  response  of  the  composite.  An  excep¬ 
tion  arises  when  the  closed  crack  faces  slide  in  shear.  Treatment  of  this  case  is 
beyond  the  scope  of  the  present  paper  [13]. 

Finally,  we  note  that  evaluation  of  stress  averages  in  the  phases  is  trivial. 
Since  open  cracks  do  not  support  any  stress, 

4 
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6.  STIFFNESS  AT  LARGE  j3 

When  a  fibrous  lamina  is  embedded  in  a  laminated  structure,  it  may  crack 
rather  extensively.  Under  such  circumstances,  certain  components  of  L„  in  (2) 
become  independent  of  /J,  or  equal  to  zero.  It  is  worth  evaluating  these 
limiting  stiffnesses  since  they  correspond  to  the  worst  possible  damage  to  a 
particular  ply  of  a  laminate.  A  good  approximation  to  the  actual  residual 
stiffness  can  be  obtained  as 

L:,  =  lim  L„  .  (37) 

P  -*  oo 

Now  as  /}  ~*  0“  it  is  easy  to  see  that 

A/«.  A/„  -*  oo 

whereas  the  remaining  components  of  the  compliance  matrix  remain  finite 
(and  equal  to  the  uncracked  value).  It  then  follows  from  simple  matrix  inver¬ 
sion  that 

LI.  = 


f.;,  =  0  1^2  =  0  L‘zs  =  0 


L'u  =  -Afls/y  =  L?,  -  LIzLWLh 

Lu  ~  A/?t/y  =  L%i  ~  (Lii/^/Lu 


where  y  *  Also, 


L\,  -  L%, ,  Z.«  =  /.«  =  0  .  (39) 

For  a  moderately  cracked  ply,  say  /)  =  1,  a  reasonable  approximation  to 
the  stiffnesses  L,.  (ij  =  1,2,3)  is  still  given  by  (38),  while  Lss  is  of  course  equal 
to  L%i .  However,  I*,  apd  must  be  calculated  either  direaly  from  (2)  and 
(A-6),  or  indirectly  from'(8),  (10),  and  (1 1). 

7.  SELECTED  RESULTS 

To  illustrate  the  evaluation  of  compliances  of  a  cracked  composite,  and  of 
the  “matrix”  strain  averages,  we  consider  a  T300  Gr-EP  system.  Table  2  lists 
the  constituent  properties  of  fiber  and  matrix,  and  of  the  uncracked  com¬ 
posite.  The  comtKJsite  compliances  were  obtained  from  self-consistent 
estimates  of  moduli  [14]. 

Figure  4  shows  changes  in  the  three  compliances  M22,  M44.  and  Mm, 
calculated  from  Equations  (8-11)  for  given  values  of  crack  density  p.  Of 


w 


■  L^ 'M I ».■! H.»  ^.■. »■  i."  k.^  ,i l/ '4.. v.TvlVTLrr 


Tabtt  Z  Coiut/tuant  pnparties  and  compllancaa  of  tha  T300  Gr-BP  systam. 


Unit 

E« 

G>. 

Vji 

E., 

G.> 

Symmetry 

eti - 

rtD0r 

10^ 

33.00 

3.36 

2.25 

0.78 

Transversely 

(T300) 

lO^MPa 

227.5 

23.2 

0.410 

15.5 

5.4 

Isotropic 

Matrix 

lOIwi 

0.50 

0.19 

0.50 

0.19 

Isotropic 

(Epoxy) 

lO^Pa 

3.46 

1.26 

0.350 

3.4 

1.26 

Compoaita  Compllancaa: 


Compllanca 

c,  «  0.2 

c,  -  0.4 

0 

1 

<r 

M?, 

0.2069 

0.1596 

0.1192 

Mf. 

-0.1037 

-0.0779 

-0.0661 

Mf, 

-0.0075 

-0.0041 

-0.0028 

M}, 

0.2069 

0.1566 

0.1192 

M!, 

-0.0075 

-0.0041 

-0.0028 

MS, 

0.0207 

0.0107 

0.0073 

M84 

0.5106 

0.2726 

0.1244 

MS. 

0.5106 

0.2726 

0.1244 

Mi. 

0.6211 

0.4746 

0.3506 
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Pgurm  4.  CompHanct  changes  in  the  T300  Gr  Ep  system  caused  by  cracks  of  density  (I. 
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Rgum  5.  Changes  in  strain  concentration  factor  components  in  composite  "matrix"  caused  by 
cracks  of  density  fl  17300  Gr-Epl. 


course,  all  these  components  increase  with  increasing  /3,  but  their  change  is 
quite  gradual,  especially  at  low  fi.  This  contrasts  but  is  not  in  conflict  with  the 
relatively  rapid  reduction  in  stiffnesses  which  we  found  in  18]. 

We  note  that  the  composite  without  cracks  is  transversely  isotropic,  and 
has  five  independent  compliances.  A/., .  When  cracks  are  introduced  the  coef¬ 
ficient  Mii,  M44,  and  change,  the  material  is  no  longer  transversely 
isotropic,  and  the  number  of  independent  elastic  compliances  increases  to 
eight. 

The  changes  caused  by  in  components  of  strain  concentration  factor  A  in 
(28)  are  shown  in  Figure  5.  As  required  by  (32),  five  of  the  eight  components 
of  A  change  with  ft.  Note  that  the  fiber  volume  fraction  appears  to  have  a 
small  effect  on  -4,, .  However,  crack  density  fi  can  have  a  significant  effect  on 
the  strain  concentration -factor  components,  especially  at  relatively  low  values 
of/). 
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SYMBOLS 


2a  =  thickness  of  a  lamina 

A,  a  =  mechanical  and  thermal  strain  concentration  factors  for 
unbroken  segments  of  a  cracked  lamina 


> 

V 


tpilQI  >  <3>  O  Q  ^3  >3  ^ 


Analysis  of  Progressive  Matrix  Cracking  in  Composite  Laminates 


231 


=  6x6  identity  matrix 
=  instantaneous  overall  stiffness  of  a  lamina 
j  “  overall  lamina  stiffness  at  large  crack  density 

=  instantaneous  overall  thermal  stress  vector  of  a  lamina 
=  thermal  stress  vector  of  the  unbroken  segments  of  a  cracked 
lamina 

=  instantaneous  overall  compliance  of  a  lamina 
=  overall  thermal  strain  vector  of  a  lamina 
=  thermal  strain  vector  of  the  unbroken  segments  of  a  cracked 
lamina 

=  axial  linear  thermal  expansion  coefficient  of  the  uncracked 
composite  lamina 

=  transverse  linear  thermal  expansion  coefficient  of  the  ' 
uncracked  composite  lamina 
=  crack  density  parameter  in  an  interior  ply 
=  2|3  =  crack  density  parameter  in  a  surface  ply 

=  ',4n/3 

=  overall  uniform  strain 

=  volume  average  of  strain  accommodated  by  cracks  (crack 
accommodation  strain) 

=  volume  average  of  strain  in  the  unbroken  segments  of  a 
cracked  lamina 

=  uniform  temperature  change 
=  crack  tensor 
=  overall  uniform  stress 

=  volume  average  of  stress  in  unbroken  segments  of  cracked 
lamina 
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APPENDIX 

Selected  Properties  of  Orthotropic  Materials 

The  constitutive  equation  for  orthotropic  materials  have  the  following 
form: 


a  =  L  t ,  t  =  M  o 


Stresses  and  strains  are; 
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Stiffness  and  compliance  matrices  are  given  in  the  form 
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Note  that  M  is  symmetric,  hence  vu/ £11  =  vu/Eu,  etc.  _ 

Stiffness  L  can  be  written  as  the  inverse  of  M.  or  as  an  expansion  of  (  ) 

terms  of  K ;  the  result  is 
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ABSTRACT 

The  mechanics  of  transverse  cracking  in  an  elastic  fibrous  composite  ply  is  explored  for 
the  case  of  lo<v  crack  density.  Cracks  are  assumed  to  initiate  from  a  nucleus  created  by 
localized  fiber  debonding  and  matrix  cracking.  Conditions  for  onset  of  unstable  cracking 
from  such  nuclei  ate  evaluated  with  regard  to  interaction  of  cracks  with  adjacent  plies  of 
different  elastic  properties.  It  is  found  that  cracks  may  propagate  in  two  directions  on 
planes  which  ate  parallel  to  the  fiber  axis  and  perpendicular  to  the  midplane  of  the  ply. 
In  general,  crack  propagation  in  the  direction  of  the  fiber  axis  controls  the  strength  of  thin 
plies,  while  cracking  in  the  direction  perpendicular  to  the  fiber  axis  determines  the 
strength  of  thick  plies.  The  theory  relates  ply  thickness,  crack  geometry,  and  ply  tough¬ 
ness  to  ply  strength.  It  predicts  a  significant  increase  in  strength  with  decreasing  ply  thick¬ 
ness  in  constrained  thin  plies.  The  strength  of  thick  plies  is  found  to  be  constant,  but  it 
may  be  reduced  by  preexisting  damage.  Strength  of  plies  of  intermediate  thickness,  and 
of  unconstrained  thick  plies  is  evaluated  as  well.  Results  ate  illustrated  by  comparison  with 
experimental  data. 

1.  INTRODUCTION 

IN  THE  FIRST  PAPER  OF  THIS  SERIES  Jll  WE  DESCRIBED  THE  EVALUATION  OF 
thermoelastic  properties  of  a  unidirectional  composite  containing  a  certain 
density  of  aligned  slit  cracks  which  grow  in  the  direction  of  the  fiber  and  across 
the  thickness  of  the  ply.  These  cracks  are  usually  referred  to  as  ply  cracks  or 
transverse  cracks.  Now  we  focus  our  attention  on  a  single  ply  in  a  typical  part 
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of  a  laminated  strucmre,  removed  from  free  edges  and  other  stress  concentra¬ 
tions.  and  examine  the  mechanics  of  transverse  cracking  in  the  ply.  In  particular, 
we  consider  the  initial  stage  of  the  cracking  process  which  consists  of  first  ply 
hulure,  and  perhaps  subsequent  hiilures  of  the  ply  at  locations  which  are  far  apart 
from  each  other,  so  that  the  cracks  do  not  interact.  Our  goal  is  to  identify  the 
failure  modes  which  cause  transverse  cracking,  to  evaluate  the  corresponding 
energy  release  rates,  and  to  present  failure  criteria  for  the  cracking  process. 

Of  course,  if  the  composite  laminate  is  loaded  beyond  the  first  ply  failure,  then 
many  transverse  cracks  may  develop  in  each  ply.  While  this  progressive  cra~king 
process  is  essentially  a  sequence  of  repetitive  ply  failures  of  the  type  discussed 
herein,  the  formation  of  each  individual  crack  is  influenced  by  interaction  with 
adjacent  cracks.  This  subject  will  be  discussed  in  a  separate  paper. 

2.  CRACK  GEOMETRY 

Consider  a  unidirectionally  reinforced  ply  of  thickness  2a.  which  is  embedded 
in  a  laminated  plate,  shell,  or  a  similar  composite  structure.  The  structure  is  sub¬ 
jected  to  a  certain  incremental  loading  program  which,  at  each  loading  step, 
causes  a  known  instantaneous  strain  state  6  in  the  ply.  A  uniform  thermal  change 
$  can  also  be  prescribed.  When  the  ply  is  free  of  cracks,  the  local  strains  r  are 
assumed  to  be  uniform,  or  nearly  uniform  in  the  sense  that  their  wavelength  is 
much  larger  than  the  ply  thickness.  The  local  strains  are  then  equal  to  their 
volume  averages,  or  overall  strains,  c.  The  given  strain  state  |  and  uniform  ther¬ 
mal  change  d  are  related  to  the  resulting  overall  stress  a  in  the  ply,  in  analogy 
with  (1-17)*.  as: 


5  =  Lo<  —  0  f  =  M.a  -t-  0  ,  (1) 

where  L,.  are  stiffness  and  compliance  matrices  of  the  uncracked  ply,  and  f„ 
m,  are  the  thermal  stress  and  strain  vectors.  We  emphasize  that  5.  |  are  ply 
averages,  not  laminate  averages. 

Suppose  that  at  a  certain  magnitude  of «,  which  is  much  smaller  than  the  failure 
strain  of  the  fiber,  a  single  crack  appears  in  the  ply.  The  crack  geometry  depends 
on  the  strength  anisotropy  of  the  fibrous  ply,  and  on  the  state  of  stress  g.  Clearly, 
cracks  may  propagate  most  easily  on  planes  parallel  to  fiber  direction.  In  a  typi¬ 
cal  part  of  a  laminated  plate  or  shell,  each  ply  usually  supports  a  plane  stress 
state.  In  the  coordinates  of  Figure  1  this  suggests  that  overall  stress  components 
O]],  djs,  and  033  are  large  in  comparison  to  d,,<  and  a,j.  It  is  easy  to  see 
that  stresses  Ou  and  du  favor  cracking  on  x,jri-planes.  The  stress  ajj  does  not 
affect  these  cracks;  it  may  cause  cracking  on  planes  perpendicular  to  the  fibers 
[2],  but  such  cracks  are  seldom  observed  in  polymer  matrix  composites.  The  re¬ 
maining  minor  shear  components  Ou  and  a,j  may  exist  in  bent  plates,  a  small 
5|j  contributes  to  crack  growth  on  x,xj-planes,  a„  does  not  have  an  effect.  Of 


*Fqittfiffnt  which  appear  in  Reference  It|  will  be  referred  to  by  ihetr  numbers  preceded  1.  c  g  .  (I '  H)  denotes 
Equation  (17)  in  |l). 
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FIgurt  1.  Slit  crack  in  a  uniformly  strained  composite  ply.  Fibers  and  crack  are  aligned  with 
X,  axis. 


course,  it  is  essential  that  a,,  be  low  or  zero  so  that  splitting  or  delamination  of 
the  lamina  or  laminate,  respectively,  do  not  take  place. 

Under  such  circumstances  the  ply  will  crack  as  indicated  in  Figure  1.  A  typical 
ply  crack  is  a  slit  crack  which  extends  in  the  fiber  direction  Xy  across  the  entire 
loaded  area  of  the  ply.  Therefore,  its  length  can  be  very  laige,  whereas  its  width 
is  essentially  equal  to  the  ply  thickness  2a.  The  laminate  layers  on  either  side 
of  the  ply  have  certain  thicknesses  h  and  b' ,  let  h  2  b' .  If  the  ply  is  an  exterior 
layer,  then  b '  =  0,  and  the  crack  becomes  a  surface  crack. 

3.  STAGES  OF  CRACK  FORMATION 

Analysis  of  ply  cracking  requires  a  reasonably  accurate  model  of  the  process. 
As  in  other  fracture  events,  one  may  anticipate  that  there  will  be  an  initiation 
stage  leading  to  formation  ^^f  a  ’ncalized  crack  nucleus,  and  a  propagation  stage 
in  which  the  crack  extends  in  the  fiber  direction  across  the  entire  loaded  area  of 
the  ply. 

In  polymer  matrix  systems,  cracks  are  often  initiated  at  locations  where  several 
fibers  debond  from  the  matrix  under  load  (31.  Such  debonded  regions  may  also 
be  present  as,  or  start  from,  initial  flaws  in  the  ply,  but  these  should  be  infrequent 
in  well-made  materials  (4).  There  is  no  need  to  assume  the  presence  of  initial 
flaws.  Instead,  fiber  debonding  can  be  regarded  as  a  consequence  of  a  statistical 
distribution  of  bond  strength,  within  a  certain  stress  interval,  along  the  length  of 
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0  CRACK  NUCLEUS 

(D  RUNNING  LONGITUDINAL  (TYPE  L)  SLIT  CRACK 

FIgun  2.  Schematic  of  crack  nucleus  and  type  L  slit  crack. 


the  fiber.  Debonding  is  also  influenced  by  stress  distribution  on  the  microscale 
which  depends  on  fiber  spacing.  Garrett  and  Bailey  [5]  point  to  numerous  studies 
which  suggest  that  debonding  is  more  likely  to  occur  in  regions  of  high  fiber  den¬ 
sity.  or  if  the  fibers  are  in  contact. 

The  net  effect  of  these  factors  on  the  total  number  of  debonded  regions  is 
difficult  to  quantify.  However,  since  the  number  of  cracks  increases  slowly, 
especially  in  the  initial  stages  of  cracking,  one  can  conclude  that  debonded 
regions  of  significant  size,  which  can  actually  initiate  cracks  at  a  given  level  of 
applied  strain,  are  not  particularly  numerous.  Clearly,  while  the  crack  density 
is  low,  the  distances  between  initiation  sites  are  many  times  larger  than  ply  thick¬ 
ness.  Therefore,  each  crack  propagates  from  a  single  initiation  site  at  which  a 
sufficiently  laige  crack  nucleus  has  formed  from  a  debonded  region. 

In  the  absence  of  experimental  information  about  the  actual  geometry  of  a 
debonded  region,  we  assume  that  when  one  or  more  adjacent  fibers  debond  from 
the  matrix,  and  these  debonds  join  up,  a  larger,  stable  microcrack  appears  in  the 
ply.  This  microcrack  lies  in  the  xiXs-plane,  where  the  ply  crack  eventually 
forms.  Figure  1.  The  in-plane  contour  of  each  fiber  debond  should  resemble  an 
elongated  ellipse,  or  an  oval,  and  it  seems  reasonable  to  adopt  the  same  contour 
for  the  initial  microcrack.  The  resulting  shape  appears  in  Figure  2.  Let  the  width 
of  the  microcrack  be  denoted  by  2d;  the  length  dimension  in  .tj,  which  may  be 
much  larger  than  26.  need  not  be  specified.  The  thickness  of  the  ply  is  2a.  and 
the  micFocrack  is  located  at  mid-thickness. 

Now,  if  the  ply  is  strained  in  an  incremental  way.  the  microcrack  will  grow 
slowly  until  it  becomes  a  crack  nucleus  of  certain  critical  width  26c ,  and  starts 
to  propagate  as  a  Griffith  crack.  The  slow  growth  process  is  probably  time- 
dependent  in  most  polymer  matrix  systems,  and  one  should  relate  the  growth  rate 
to  loading  history.  The  rapid  crack  propagation  can  be  regarded  as  inviscid,  at 
least  in  brittle  resin  matrices. 
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4.  ENERGY  RELEASE  RATES 

To  evaluate  the  critical  size  26  of  the  crack  nucleus  under  a  given  stress,  and 
to  determine  the  direction  of  crack  propagation,  it  is  necessary  to  derive  expres¬ 
sions  for  energy  release  rates  of  the  crack  nucleus  in  a  ply  at  constant  overall 
strain.  Clearly,  the  crack  may  extend  either  in  the  x,  direction,  in  the  Xj  direction, 
or  simultaneously  in  both  directions.  To  distinguish  between  the  direction  related 
quantities,  we  indicate  the  relevant  direction  by  the  letter  L,  or  T,  written  in 
parentheses  next  to  the  appropriate  symbol.  For  example,  the  energy  release  rate 
G(T)  refers  to  crack  extension  in  the  transverse  direction  x,,  and  G(L )  to  crack 
extension  in  the  longitudinal  direction  Xy.  The  corresponding  toughnesses  are 
GAT),  GAL),  critical  crack  widths  dAT),  6AL),  etc. 

In  general,  we  shall  refer  to  a  crack  that  propagates  in  the  x,  direction  as  type 
T  crack,  and  to  that  which  propagates  in  the  jc,  direction  as  type  L  crack. 

It  is  probably  obvious  that  due  to  an,  both  crack  types  open  in  the  Xi  direction, 
or  in  Mode  I.  In  addition,  there  are  two  shear  modes,  one  along  the  fiber  direc¬ 
tion.  and  one  transverse  to  the  fiber  direction.  These  will  be  designated  for  both 
crack  types  as  follows;  Mode  II  represents  longitudinal  shear  displacement  in  x,, 
caused  by  an.  and  Mode  III  is  the  transverse  shear  mode  in  Xj  caused  by  a,2. 
This  notation  is  derived  from  the  type  L  cracking  and  adopted  for  type  T  cracks. 
Subscripts  I,  II,  and  III  will  denote  the  modes.  While  the  crack  front  may  change 
direction,  the  mode  designation  remains  related  to  stress  components  in  a  fixed 
coordinate  system. 

Consider  first  the  case  of  a  thick  ply  with  a  small  crack  nucleus  d  <<  a. 
Figure  2.  The  nucleus  is  not  affected  by  the  presence  of  adjacent  plies,  and  can 
be  regarded  as  a  slit  crack  in  an  infinite  orthotropic  medium.  In  this  case,  the 
energy  release  rates  for  both  type  T  and  type  L  can  be  obtained  from  the  interac¬ 
tion  energy  H'i,  which  is  equal  to  the  total  energy  released  by  a  slit  crack  of  unit 
length  X3  and  width  6  under  constant  e.  Expressions  for  Wi  in  anisotropic  solids 
were  derived  by  Stroh  [6]  and  in  the  orthotropic  case  by  Laws  (Tj.  In  the  notation 
used  in  [1];  for  6  <  <  a, 

—  2  tt  5^[Ali  all  +  A54  oJj  -I-  A«6  Uijj  (2) 

Here.  are  nonvanishing  components  of  the  crack  tensor  derived  in  [1],  taken 
for  a  dilute  concentration  of  cracks,  i.e.  for  crack  density  parameter  0  —  0.  With 
the  help  of  (1-4)  to  (1-6)  and  (1-A5)  one  can  write 


where  £i  =  En,  Et  =  Eu  =  £ji  are  the  longitudinal  and  transverse  Young's 
moduli,  Gi  =  L*,  =  Lss  in  [1]  is  the  longitudinal  shear  modulus,  and 
vi  =  v^l  =  vii  is  the  corresponding  Poisson’s  ratio. 
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The  energy  release  rate  of  a  crack  of  length  25,  for  an  elementary  extension 
in  the  Jt,  direction  in  an  unbounded  orthotropic  solid  is  equal  to  '/i  dM^t/dS,  or 

G(T)  =  J  ^  §[^22  ff22  +  Ai4  oil  +  AS»  ai2l,  for  6  <  <  a.  (4) 

One  can  easily  verify  that  (2)  and  (4)  reduce  to  their  familiar  plane-strain 
forms  if  the  solid  becomes  isotropic. 

We  now  turn  our  attention  to  extension  of  an  oval  crack  of  width  26  in  the  fiber 
direction's.  Figure  2.  The  energy  release  rate  G(L)  of  this  (type  L)  crack  is  equal 
to  the  interaction  energy  1^,  released  by  a  crack  of  width  25  as  it  extends  a 
unit  distance  Xj,  and  thus  creates  a  new  surface  area  254,1.  This  is  the  amount 
of  energy  change  caused  by  unit  extension  of  the  prismatic  part  CC  of  the  stage 
2  crack  in  Figure  2.  Hence, 

G(L)  =  Ws/26  (5) 

or, 

G(Z.)  =  i  j  CHT)  db  (6) 

From  (2)  and  (4)  one  obtains 

G{L)  =  ^  f  6(Al2  oil  +  A44  ffii  +  AXe  Sfsl  for  5  <  <  a  (7) 

Note  that  G(L)  does  not  depend  on  the  actual  length  of  the  crack  nucleus  in 
x,  direction,  i.e.,  of  the  segment  CC  in  Figure  2.  The  crack  width  6  is  the  only 
parameter  reflecting  crack  geometry,  and  if  its  value  is  fixed,  then  CiL)  is  con¬ 
stant. 

Consider  now  the  case  of  a  crack  nucleus  in  a  ply  of  intermediate  thickness 
0  <  8/a  <  I.  As  the  crack  approaches  the  interface  between  its  own  ply  and  an 
adjacent  ply  of  different  elastic  properties,  W,  is  influenced  by  the  proximity  of 
the  adjacent  ply.  In  the  direction  x,,  perpendicular  to  the  crack  plane,  the 
adjacent  plies  are  usually  much  suffer  than  the  ply  under  consideration,  and  this 
causes  a  reduction  in  W,  as  6/a  —  1.  To  reflect  this,  (4)  can  be  written  in  the 
form  valid  for  type  T  crack  of  any  length 

G(T)  =  2  ^  8[rffAU^i  -h  rffiAuals  +■  Vm  oh)  for  (0  <  6/0  <  1)(8) 


where 
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are  stress  intensity  reduction  coefficients  which  depend  on  crack  geometry  and 
elastic  properties  of  layers.  The  stress  a  =  for  t  =  /,  etc. 

For  a  unit  extension  in  the  Xj  direction  of  a  crack  nucleus  of  width  26,  Figure 
2,  when  0<  6/a <  1,  the  energy  release  rate  is  given  again  by  (5)  and  (6).  How¬ 
ever,  since  the  interaction  energy  is  now  influenced  by  the  proximity  of  the 
next  ply,  (7)  now  becomes,  for  0  <  6/a  <  1: 


CHL)  ~  All  (til  +  ill  Ojj  -t-  iiii  Ait  oil] 


where,  from  (6)  and  (8); 


-  -  i  =  i.ii,in 


Again,  the  effect  of  crack  geometry  is  reflected  only  through  6,  and  if  that  is 
fixed,  say  at  6  =  a,  then  GfZ,)  is  constant. 


0  0.2  0.4  0.6  0.0  I.O 

a 


Figufw  3.  Stress  intensity  reduction  coefficients  7,  =  K,/(a  Vi ;  tor  the  three  fracture  modes. 
The  longitudinal  shear  mode  is  denoted  by  i  =  II. 
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Coefficients  jj,  can  be  ft)und  for  many  crack  configurations  in  nonhomogeneous 
media.  For  transverse  cracks  in  layered  media  one  can  refer  to  the  results  of  Hil¬ 
ton  and  Sih  [8),  and  Cook  and  Erdogan  {9}  in  the  case  of  isotropic  layers.  Similar 
results  for  laminates  made  of  orthotropic  layers  can  be  found  in  [10],  and  a  finite 
element  evaluation  of  tj,  for  a  graphite-epoxy  laminate  appeared  in  {Uj. 

These  and  related  results  indicate  that  n.  are  influenced  by  the  differences  in 
elastic  constants  between  the  layers,  and  particularly  so  by  the  ratio  Eli ’/£]]>  of 
the  Young’s  moduli  of  the  cracked  layer  (1)  to  that  of  adjacent  layer  (2),  in  direc¬ 
tion  X,,  perpendicular  to  crack  plane  x^x,.  Additional  effects  of  elastic  constants 
exist  and  are  discussed  in  [10).  A  significant  feature  of  these  results  is  that  17,  are 
not  very  sensitive  to  values  of  the  moduli  ratio  for  EiJ’/Eli’  <  0.1  or  so,  the 
outer  layer  is  perceived  by  the  crack  as  very  stiff  or  rigid.  The  relative  thick¬ 
nesses  of  the  layers  can  have  an  impc>rtant  effect,  particularly  in  thin  laminates 
where  >?,  increase  if  the  crack  interacts  with  a  surface  ply  (8.9,12).  Thermal  stress 
effects  in  cracked  laminates  have  apparently  not  received  much  attention,  except 
for  [il),  and  at  this  time  can  be  examined  only  in  terms  of  their  contribution  to 
overall  applied  stress  d,j,  via  Equation  (1). 

For  illustration.  Figure  3  shows  the  coefficients  tf,  for  a  crack  in  a  laminate  con¬ 
sisting  of  isotropic  layers,  for  =  0.1.  v  =  0.3  (8.9).  These  coefficients 

can  be  used  as  an  approximation  for  laminates  consisting  of  transversely  isotropic 


0 


Figun  4.  Interaction  energy  reduction  coefficients  £.  for  tfis  tftrae  fracture  modes. 
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layers,  providing  that  one  takes  =  £F’  and  Ei}'  =  The  selected  ratio 
of  moduli  corresponds  to  that  of  (0/90)  graphite-epoxy. 

Figure  4  shows  values  of  (5/a)£  calculated  from  the  rj,  in  Figure  3  and  (11). 
These  refer  to  a  crack  which  has  a  certain  width  6.  and  are  plotted  as  functions 
of  8/a. 

We  note  that  17,  are  generally  available  only  up  to  81a  =  0  90.  Also,  the  order 
of  stress  singularity  changes  when  the  crack  reaches  the  interface. 

As  in  the  case  of  77, ,  the  estimates  of  £  in  Figure  4  are  rather  conservative, 
taken  for  adjacent  plies  of  very  different  properties  such  as  one  encounters  in 
(0/90)  graphite-epoxy.  Somewhat  higher  values  can  be  expected  in  glass-epoxy 
laminates,  and  also  in  thin  laminates  where  the  crack  interacts  with  one  or  two 
surface  plies.  For  example.  Chou  et  al..  (11)  give  values  for  (±  25/90j),  and 
(0j/90j).  in  AS-3501-06  graphite-epoxy  system  which  are  about  25%  higher  than 
those  shown  in  Figure  4. 

Finally,  consider  the  case  of  a  thin  ply  in  which  the  crack  nucleus  extends 
across  the  entire  ply  width,  i.e..  6/a  =  1.  Only  CHL)  is  of  interest  in  this  case, 
and  it  can  be  found  from  (10)  provided  that  £  are  known.  Unstable  type  T crack¬ 
ing  cannot  occur,  however;  one  should  keep  in  mind  that  the  crack  tips  at 
.r,  =  ±a  reside  in  a  soft  matrix  which  may  not  be  able  to  support  the  high  crack 
tip  stresses.  Therefore,  local  crack  extension,  across  or  along  the  ply  interface, 
may  take  place  to  relieve  the  stress  singularity.  Such  localized  crack  extensions 
have  been  observed  experimentally  [5).  More  recent  studies  suggest  that 
delaminations  as  well  as  numerous  fiber  breaks  can  be  found  in  the  axial  plies 
next  to  transverse  cracks  (13).  These  observations  are  not  surprising  in  view  of 
analytical  studies  of  cracks  terminating  at  interffices,  which  suggest  that  the  stress 
intensity  ffictor  in  the  next  ply  and  the  magnitude  of  tension  stresses  perpendicu¬ 
lar  to  the  interface  both  increase  at  higher  ratios  of  £ii'’/£jj*  [9,12).  Therefore, 
one  should  expect  that  in  the  presence  of  additional  cracking  at  ply  interface  the 
coefficients  are  actually  larger  than  those  predicted  theoretically  for  cracks  ter¬ 
minating  at  the  interface.  Until  the  effect  of  interface  damage  is  clarified,  one  can 
use  a  simple  estimate,  e  g.,  =  1.  as  in  [1). 

For  a  selected  £  one  can  find  GiL)  of  a  thin  ply  from  (10).  Note  that  (8)  and 
(9)  are  not  affected  by  interface  damage  but  (6)  and  (11)  no  longer  apply. 


5.  FAILURE  CRITERIA 


If  a  composite  ply  is  loaded  only  by  a  simple  stress  state,  such  as  On  or  6,, 
acting  alone,  then  unstable  propagation  of  type  T  or  type  L  cracks  takes  place 
when  the  appropriate  value  of  G,  or  G„  equals  or  exceeds  the  corresponding  ply 
toughness  G,c  or  Gm,  respectively. 

It  is  useful  to  point  out  that  fracture  surface  topography  depends  on  fracture 
mode.  Modes  I  and  II  in  a  composite  ply  always  create  asperities  on  the  fracture 
surface,  as  the  crack  seeks  its  way  between  densely  packed  fibers.  These 
asperities  can  be  interlocked  and  thus  impede  mode  111  cracking.  Also  in  a 
typical  part  of  laminated  plate  or  shell,  the  corresponding  stress  6,2  is  usually 


• .  •  ^  *1  • 


318 


George  J.  Dvorak  and  Norman  Laws 


small  or  equal  to  zero.  Therefore,  mode  HI  need  not  be  considered  in  many  ap¬ 
plications.  For  simplicity,  we  neglect  mode  III  in  the  sequel. 

Apart  from  the  creation  of  asperities,  fracture  surface  topography  is  quite  dif¬ 
ferent  in  modes  I  and  II.  Mode  II  favors  formation  of  hackles,  and  higher  energy 
absorption.  For  example,  a  recent  review  of  experimental  results  [14]  quotes 
toughness  values  for  graphite-epoxy  and  glass-epoxy  for  which  K„jKu  =  5. 

In  mixed  mode  fracture,  the  toughness  values  are  interrelated.  Specific  criteria 
which  describe  this  effect  need  to  be  derived  from  experimental  data.  Hahn  [14] 
offers  a  form  based  on  fracture  mechanics  considerations  that  fits  data  obtained 
for  certain  glass-epoxy  and  graphite-epoxy  systems; 


where  g  =  Gu  'G,u. 

If  G,  are  independent  of  loading  mode,  then  g  —  I,  and  (12)  as.sumes  a  familiar 
form.  In  actual  composite  plies  with  epoxy  matrices  g  =  0.1  [14). 

More  complex  criteria  can  be  derived  as  needed.  In  what  follows  we  shall  use 
(12)  as  a  representative  mixed  mode  criterion,  both  for  type  T and  type  L  cracks. 
However,  we  allow  for  possible  differences  between  toughness  values  in  the  two 
modes  by  writing; 


G,c(L)  =  7,  GiAT) 

GiiAL)  =  7//  GiiAT) 


(13) 


Results  reported  in  (4|  and  (15]  indicate  that  y,  ranges  from  0.6  to  0.8  for  certain 
graphite-epoxy  plies.  Data  of  this  kind  are  apparently  very  rare,  therefore  no  gen¬ 
eral  conclusions  about  y,  and  y,i  magnitudes  can  be  drawn  at  this  time. 

6.  FAILURE  OF  THICK  PLIES 

Consider  a  ply  in  a  composite  laminate  which  is  subjected  to  a  certain 
prescribed  strain,  such  that  within  the  ply  the  stress  state  is  Bn  =  q  and  is 
applied  incrementally  until  first  ply  failure.  The  ply  fails  by  unstable  propagation 
of  a  crack  nucleus.  Figure  2,  when  the  width  5  grows  to  a  certain  critical  magni¬ 
tude  5c  ■  We  recall  that  5  may  increase  slowly  under  load,  possibly  as  a  function 
of  time,  especially  during  cyclic  loading 

5  =  6[f(t)|  for  5  <  5,  (14) 

but  once  6,  has  been  reached,  the  propagation  process  is  rapid  and  regarded  as 
time-independent. 

The  thickness  2a  of  the  ply  is  large,  so  that  5rla  <  <  1. 

Our  objective  is  to  determine  the  relationship  between  6c  and  the  measured  ply 
strength,  and  the  directional  sequence  of  the  cracking  process. 
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Gi{T)  =  2  5 

Gii(T)  =  ',  It  8  A«<  oil 


(15) 


are  the  energy  release  rates  for  type  T  cracking  in  jc,  direction.  Also,  from  (7), 


Gi(L)  —  ^  IT  6  Ajj  (Tij 
GtiiL)  =  ^  6  A44  oh 


(16) 


are  the  rates  for  type  L  cracking  in  j:j  direction. 
At  any  given  stress  state. 


GAL)  =  j  GAT) 

for  6  <<  a 

G,AL)  =  ^G,AT) 


(17) 


Suppose  that  for  the  special  case  of  loading  ou  ^  0.  Ou  -  0,  the  ply  fails  at 
ou  =  (5ii)cr.  At  failure  load,  at  least  one  of  the  following  conditions  must  be 
satisfied: 


GAT)  >  G,AT).  S  >  5,.(r) 

and/or,  (18) 

GAL)  >  G,AL).  8  2  8,AL) 

where  8, AT),  81AT)  are  critical  widths  of  the  crack  nucleus  in  Figure  2,  at  onset 
of  crack  propagation  in  directions  Xi,  or  x,,  respectively. 

When  (18)  are  taken  as  equalities,  it  follows  from  (13)  and  (15)  to  (17)  that 

8, AT)  =  2G,AT)lTrMi(dii)h  (19) 


^iAL)  —  4Gfc(i')/TrAii(<7ji): 


(20) 
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We  remark  that  a  slightly  different  form  of  (19)  has  been  obtained  by  Hahn  and 
Johnnesson  in  [14).  Since  the  available  data  indicate  that  2y,  >  1,  then 
5i.(l)  <  5,c(3).  Therefore,  first  ply  failure  in  thick  plies  occurs  as  a  result  of  type 
Tcmctdng  followed  by  type  L  cracking.  The  crack  starts  to  propagate  as  the  type 
r  crack  in  the  jc,  direction.  After  the  crack  has  propagated  under  constant  stress 
ou  =  (au)cr  to  width  5  =  6u(L),  it  also  may  start  to  propagate  in  the  .v,  direc¬ 
tion,  as  the  type  L  crack.  This  bidirectional  propagation  continues  until  26  has 
reached  a  certain  maximum  value  which  may  not  exceed  the  ply  thickness  2a. 
but  may  be  smaller  than  2a  if  the  type  T  crack  becomes  arrested  because  of  its 
interaction  with  an  adjacent  ply.  In  any  case,  the  type  L  cracking  continues  in 
the  direction  at  any  5  >  8„(L)  across  the  entire  area  of  the  ply  where 
>  (oiihr-  We  recall  that  G(L)  does  not  depend  on  crack  length,  hence  the 
type  L  crack  may  be  arrested  quite  easily  if  the  crack  enters  a  region  of  low  stress 
or  locally  higher  toughness. 

Next,  suppose  that  the  ply  is  loaded  in  pure  shear,  on  =  0.  an  ^  0. 
Analogous  to  the  tension  case,  one  finds  at  failure  load  Ojj  =  (ojj)er  that 

G,AT)  >  G„.(7).  62  6„AT)  (22) 

and/or 

G,X£.)  2  G„AL),  52  SnAL)  (23) 

where 

5„,(r)  =  2G„AT)lir  \Uan)l 
h„AL)  =  4G„.(L)/t 
and 


b„AL)  =  2y„  b„AT)  (24) 

Now,  if  one  defines  g  in  (12)  for  the  two  crack  types  as 

g(T)  =  G,AT)IG„AT)  (25) 

g(L)  =  G,AL)IG„AL)  (26) 

then  it  follows  that 

6, AT)  =  8(T)  S„AT)  (27) 

d,AL)  =  glL)  8„AL)  (28) 


Recall  that  g  =  O.l  according  to  [14]  for  certain  graphite  epoxy  systems.  This 
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value  is  actually  g{L).  no  information  seems  to  be  available  about  giT).  However, 
it  appears  that  the  critical  width  of  the  crack  nucleus  in  thick  plies  is  generally 
smaller  for  Mode  I  cracking. 

Again,  if  the  ply  is  loaded  by  simple  shear,  the  crack  nucleus  starts  to  propa¬ 
gate  in  the  .Ci  direction  if  lyn  >  1  and  6,iAT)  <  buAL).  Bidirectional  propaga¬ 
tion  may  start  when  6  =  buAL).  As  in  the  previous  case.  GiLi  is  constant  at  fixed 
6.  Type  L  cracking  continues  across  the  loaded  area  of  the  ply  as  long  as 

Oiy  > 

Finally,  we  consider  the  general  case  of  loading  oiy  =  qdn  in  the  ply,  which 
is  applied  until  failure  at  =  ^ilonh,.  Critical  width  of  the  crack  nucleus 

must  be  found  from  (12)  for  each  q.  For  type  T  cracking  one  can  write 

[1  -  y( r )((/?,( T) I""  -f-  g(TlR,(T)  ■+  RiAT)  2:  I  (29) 


where 


R,(T)  =  GAT  MG, AT) 
R,AT)  =  G,AT)IG„AT) 


(30) 


Equation  (29)  can  be  solved  to  yield  6.(T), 

Similarly,  for  type  L  cracking,  the  equation  that  can  be  used  to  find  6AL)  is 

(1  -  g(L)U/?,(E)r''  -t-  g(L)  RAL)  +  R,AL)  >  1  (31) 

where  RAL).  R,AL)  are  defined  in  analogy  with  (30). 

Type  T  cracking  precedes  type  L  if 

SAT)  <  sad  (32) 

and  vice  versa.  In  any  event,  bidirectional  crack  propagation  starts  when 

6  =  max!dc(T).  SAD\  (33) 

The  results  derived  so  far  suggest  that  if  many  thick  plies  of  different  thickness 
are  tested  in  such  a  way  that  the  loading  history  |(f )  is  kept  constant  in  each  test, 
then  6  in  (14)  and  (33)  is  also  constant,  and  the  measured  strength  must  be  con¬ 
stant.  regardless  of  actual  ply  thickness. 

For  the  special  case  of  loading  by  Ot,  =  0.  ojj  ^  0.  the  strength  is  given  by 
the  smaller  of  the  two  expressions 


(5„)„  =  (2C,.(7')/ir  \h  5, .(71)'" 


(34) 


=  \AG,AL)ltt  5,.(L))'" 


‘-’T  1  •  V.  'X'? 


P 


I 


I 


jg 

•S 

f 


322 


George  J.  Dvorak  and  Norman  Laws 


For  5,1  =  0,  5,3  ^  0,  the  strength  is  given  by  the  smaller  of  the  two  expres¬ 
sions 


(^23)0,  =  l2Gnc(T)/r  Ai,  5„c(r)]''' 
(®J3)«r  =  [4G//e(L)/')r  A44  5//c(£-)]’^* 


(35) 


Constant  strength  is  seen  in  experimental  results  discussed  in  the  sequel.  Simi¬ 
lar  relations  can  be  found  from  (29)  and  (31)  for  combined  loading. 

The  strength  of  unidirectional  laminates,  or  unconstrained  thick  plies,  can  be 
estimated  along  similar  lines  by  assuming  that  the  crack  nucleus  may  most 
probably  form  at  the  surfece.  If  the  rate  of  slow  growth  of  6  depends  only  on  ap¬ 
plied  strain  history  as  in  (14),  one  should  anticipate  that  crack  nuclei  of  whidth 
26c  may  be  found  at  the  surface  of  unconstrained  unidirectional  laminate,  as  well 
as  in  the  interior.  This  implies  that  the  estimate  of  transverse  strength  of  un¬ 
constrained  plies  is 

(a„);,  =  (5„)„/(1.12Ny?)  (36) 

where  (5,,)„  follows  from  (34)  and  vice  versa.  The  multiplier  1.12  accounts  for 
the  stress  intensity  magnification  of  a  surface  crack. 

7,  STRENGTH  OF  PLIES  OF  INTERMEDIATE  THICKNESS 

Here  we  consider  plies  of  any  thickness  2a,  subjected  to  a  certain  strain 
history,  such  that  at  onset  of  unstable  crack  propagation,  the  critical  width  of  the 
crack  nucleus  26,:  is  comparable  to  2a,  i.e. 

0  <  6cla  <  1. 

To  evaluate  this  condition  in  an  approximate  manner,  one  may  utilize  5,  found 
for  a  thick  ply  of  the  same  system. 

We  recall  that  according  to  (8) 


G^T)  =  I  ir  5  jjfAl,  5i, 


Gii{T)  —  2  ^  ’t?/AS4  ui. 


(37) 


are  the  rei.ase  rates  for  type  T  cracking  in  x,  direction,  and  that  from  (10),  for 
'ype  L  cracks. 


Gi{L)  —  ^  IT  6  All  Oiz 


Gii(L)  —  ^  T  6  {/,  AJ4  01} 


(38) 
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At  any  given  stress  state. 

G,(L)  =  ^  (i,/r,})  G.m 

for  0  <  5/ci  <  1  (39) 

Gii{L)  =  ^  G//(7*) 

Consider  first  the  special  case  of  loading  Oi,  *■  0,  Sjj  =  0.  At  failure  load 
On  =  (o2i),,.  the  conditions  (18)  must  again  be  met.  When  (18)  are  taken  as 
equalities,  it  follows  from  (13)  and  (37)  to  (39)  that 


5,J  T )  =  2G;.(  T )/ir  rji  Alt  (ffji)^ 


5/c(  f- )  —  ^Gui  L)l  Tt  A22  ((t2l)c. 


(40) 


and 


&iAL)  =  2y,{^, ini)  d,AT)  (41) 

Hence,  it  2yA^iirii)  >  1.  then  6/J  T)  <  &,AL).  and  the  crack  starts  to  propa¬ 
gate  as  a  type  T  crack  in  the  .t,  direction,  and  vice  versa.  In  any  case,  bidirec¬ 
tional  propagation  starts  when 

5  =  max  b,AT).b,AL),  (42) 


Analogous  formulae  can  be  found  for  the  special  case  Sjj  =  0.  Sjj  0: 


bnAT)  -  2G,u(TUir  rii,  Ait 
biiAL)  ~  AGiuAL)/ It  A34 


(43) 


and 


b„AL)  =  2y„{^„/rih)  6„AT)  (44) 

Under  general  loading,  ajj  =  q  bn.  which  continues  until  a  critical  stres.s 
level  (bi,),.  =  q<bii)„  has  been  reached,  it  is  neces.sary  to  find  bAT)  and  bAL) 
from  (12).  The  equation  for  bAT)  is  formally  identical  with  (29)  and  (30)  but 
GaT).  GiAT)  must  now  be  taken  from  (37).  Also,  for  bAL)  one  can  use  (31)  with 
GAL)  and  G,AL)  from  (38).  Once  bAT)  and  bAL)  have  been  found  from  the 
measured  ply  strength  and  toughness,  one  can  find  that  type  T cracking  precedes 
type  L  if 


bAT)  <  bAL) 


(45) 
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and  vice  versa.  Bidirectional  propagation  starts  when 

5  =  max  16.(7).  5.(L)|  (46) 

In  theory,  one  can  use  inversion  of  (40)  and  (43)  in  conjunction  with  Figures 
3  and  4  to  calculate  ply  strength.  Current  experimental  data  is  insufficient  to 
justify  a  detailed  analysis  of  (40)  and  (43),  as  well  as  (41)  and  (44). 

The  propagation  process  of  both  types  of  cracks  is  affected  by  crack  interaction 
with  the  adjacent  ply.  For  type  7  cracking  this  is  evident  from  Figure  3,  where 
(rji  61a)  reaches  a  maximum  at  a  certain  value  of  (6/n).„  and  then  decreases, 
i^cording  to  (37),  G(7)  also  decreases.  Therefore,  one  may  find  that 
G(7)  <  Gc(7)  along  the  crack  path,  and  the  crack  will  stop.  In  fact,  if  6c(T)/a 
coincides  with  or  is  larger  than  {6/a)„^,  then  unstable  type  7  cracking  may  be 
prevented  altogether,  and  only  slow  growth  may  take  place  in  jc,  direction  until 
5  =  6.(L). 

On  the  other  hand,  type  L  cracks  may  propagate  across  the  entire  loaded  area 
of  the  ply  as  soon  as  6  =  6.(L).  It  is  not  certain  that  the  crack  retains  its  initial 
width  during  propagation.  The  width  will  probably  increase  to  26  =  la  because 
this  will  enhance  G,{L)  and  GtAL).  The  width  adjustment  may  occur  even  during 
slow  growth  of  the  crack  nucleus  in  direction.  The  crack  will  probably  try  to 
fan  out  across  ply  width,  especially  when  y,  and  yi,  are  small.  This  is  a  three- 
dimensional  crack  problem  that  cannot  be  resolved  on  the  basis  of  two- 
dimensional  solutions  which  have  been  used  therein. 

&  STRENGTH  OF  THIN  PLIES 

Here  we  consider  plies  of  any  thickness  2a,  subjected  to  a  certain  strain 
history,  such  that  at  onset  of  unstable  crack  propagation 

26.(1)  =  2a  (47) 

It  follows  from  (40j),  (432),  and  their  mixed  mode  equivalents  that  the  strength 
of  thin  plies  is 

(52i)«r  —  (4G/t(£.)/ir  A22  aj*'^*  ,  at  0^  —  0 

(48) 

=  l4G„.(iL)/x-  Ai,  a]*'*  .  at  621  =  0 

In  the  general  loading  case  ^22  =  <i  one  can  calculate  (aii)cr  from  (3)), 
providing  that  GAL)  and  G/AL)  are  taken  from  (38)  with  6  =  a. 

Accordingly,  the  strength  of  thin  plies  becomes  related  only  to  onset  of  unstable 
type  L  cracking,  much  as  the  strength  of  a  thick  ply  was  related  to  type  7  crack¬ 
ing. 

If  experimentally  measured  strength  data  are  available  for  thin  plies  of  variable 
thickness,  then  one  can  always  find  G,.(L)  and  Gn.(L)  of  the  composite  system 
in  question. 

At  this  stage  it  may  well  be  necessary  to  consider  the  possibility  of  interface 
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damage  discussed  at  the  end  of  Section  4.  As  a  consequence.  {,  may  increase 
beyond  their  theoretical  estimate.  In  the  absence  of  adequate  experimental  infor¬ 
mation  one  can  tentatively  take 


=  l.  and  hAL)  =  a  (49) 

With  this  adjustment,  one  can  calculate  ply  strength  from  (48)  and  (49)  providing 
that  GiAL)  and/or  G„AL)  are  known.  In  the  general  loading  case  oij  =  q  On  one 
can  calculate  (aj2)tr  from  (31).  if  GaL)  and  G„(Z.)  are  taken  from  (38)  at  6  =  a. 

As  long  as  the  strength  of  plies  of  intermediate  thickness  cannot  be  completely 
analyzed  with  available  data,  it  .suffices  to  consider  only  thick  and  thin  plies.  The 
transition  between  these  two  extremes  can  be  defined  as  the  intersection  of 
strength  curves  calculated  from  (34),  (35),  and  (48),  respectively.  Since  at  either 
side  of  this  transition  point  the  ply  strength  is  controlled  by  either  type  (T)  or 
type  (Z.)  cracking,  there  will  be  a  discontinuous  change  in  the  critical  width  of 
crack  nucleus  from  5aT)  to  6AL).  Of  course,  if  it  were  possible  to  analyze  plies 
of  intermediate  thickness,  this  discontinuity  would  be  bridged. 

9.  COMPARISON  WITH  EXPERIMENTS 

To  illustrate  the  theoretical  derivation,  we  interpret  results  presented  by  Bailey, 
et  al.  (3),  and  by  Crossman.  Wang,  et  al..  (15,16).  Figure  5  shows  the  first  set 


figure  5.  Comparison  of  thaoretical  rasults  with  experimental  data. 
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Table  1.  Material  properties  of  selected  piles. 

E-Glaas/Epoxy 

T300/934 

E,  (GPa) 

42 

163.4 

Et  (GPa) 

14 

11.9 

>'1. 

0.27 

0.3 

(10-‘/»C) 

4.3 

0.35 

ay  (10  V»C) 

14.3 

28.8 

A9  ("C) 

-125 

-125 

of  results.  Several  (0/90,),  coupons  made  of  a  E-glass/epoxy  laminated  plate  were 
loaded  in  tension,  and  the  laminate  stress  or  strain  was  recorded  at  first  failure 
of  the  90°  ply  (c.f..  Figure  13,  p.  616  in  [3]).  To  interpret  the  data,  laminated 
plate  theory  was  used  to  calculate  initial  thermal  stresses  in  the  90°  ply  after 
cooling  from  the  curing  temperature  to  room  temperature,  Ad  =  — 125  °C. 
Then,  the  transverse  normal  stress  au  caustui  in  the  90°  ply  by  loading  of  the 
laminate  was  found,  and  superimposed  with  thermal  stress.  In  this  way,  the 
original  data  points  were  converted  to  ply  stress  022.  and  ti>en  plotted  in  Figure 
5.  Thermoelastic  properties  of  the  ply  were  taken  from  Table  1;  these  were  found 
in  Table  1  in  [3],  and  confirmed  ^  comparison  of  calculated  and  measured 
laminate  properties. 

Next,  ply  toughness  G/^(L)  —  250  J/m‘  was  found  by  inversion  of  (48),  such 
that  the  resulting  ( curve  fitted  the  experimental  points  at  low  ply  thick¬ 
nesses. 

Then,  experimental  points  for  thick  plies  were  used  to  find  the  average  strength 
of  thick  plies  as  (aaa)cr  =  0089  GPa.  This  strength  exceeds  the  value  measured 
on  free  ^°  plies,  which  was  found  as  0.056  GPa,  and  must  be  converted  to  con¬ 
strained  90°  ply  strength  of  0.056*  (1.12  V?)  =  0.0885  GPa,  with  the  help  of 
(36).  The  prediction  agrees  well  with  the  measured  values.  The  transition  be¬ 
tween  thin  and  thick  plies  is  found  at  =0.5  mm  in  this  case. 


Table  2.  Critical  widths  of  type  T  crack  nuclei  In  thick  piles, 
calculated  from  (19). 


E-GlaM/Epoxy 

T300/934 

Gfc(T)  (J/m’) 

ibJJl  (mm) 

2i,.(T)  (mm) 

200 

0.231 

0.372 

400 

0.461 

0.745 

600 

0692 

1.117 

800 

0.922 

1.489 

Thick  ply  strength 

(GPa) 

0.089 

0.064 

AA  ii’‘r  tfi  iT^ 


fv*-- 
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IT300/9341 


O  (t25/90n)s 
tt  (252/-252/902), 
««  (90^), 

6  Onset  of  delomi- 


(0/90„/0) 


,Gic(L)=220  J/m2 


(5=22)  cr 


O-  - O - 


o  Crossmon,  Wong  (1982) 
•  Wang  (1984) 


INNER  90"  PLY  THICKNESS  2a  (mnn) 
f/gure  6.  Comparison  of  theoretical  results  with  experimental  data. 

No  definite  evaluation  of  6u(T)  can  be  made,  because  G,AT)  was  not 
measured.  However,  expected  values  at  specified  Gu  are  shown  in  Table  2.  For 
example,  if  GiAT)  is  taken  as  equal  to  G/.(L)  =  250  J/m^,  then 
25, AT)  =  0.288  mm.  In  any  event.  5, AT)  <  a.  At  the  transition  point  between 
thin  and  thick  plies  2a  =  0.5  mm,  hence  26, .(F)  <  0.5  mm.  At  the  measured 
strength  of  thick  plies  (0.089  GPa),  this  suggests  that  G,AT)  <  433  F/m'.  In 
reality,  both  26, AT)  and  G,AT)  should  be  much  smaller,  so  that  type  T cracking 
may  take  place  even  as  ply  thickness  approaches  the  transition  value  of  0.5  mm. 

Figure  6  shows  similar  results  for  a  T300/934  graphite/epoxy  system,  which 
appeared  in  (15,16).  Data  for  three  different  layups.  (±  25/90,),, 
(25,/ —  25,/9U,),,  and  (0/9u,/0),  were  superimposed  to  find  G,AL)  in  the  90° 
plies.  The  results  are  apparently  not  affected  by  the  differences  in  layup  of  outside 
layers.  Also  shown  is  the  strength  of  a  free  (90,),  ply,  which  again  conforms  with 
(36).  The  G,c(Z-)  values  are  somewhat  lower  than  those  found  hy  Slepeiz  and 
Carlson  [17]  for  glass  and  graphite-epoxy  systems. 

Note  that  the  G„  values  shown  in  Figures  5  and  6  are  actually  values  of  G,Ah 
with  taken  as  equal  to  unity,  because  of  possible  interface  damage.  If  i,  was 
taken  as  equal  to  calculated  values,  e  g.  0.8  in  E-glass/epoxy,  then  one 
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would  obtain  Gu  =  200  J/m^  in  Figure  5,  which  compares  better  with  the 
Gtc  =5  120  ±  30  7/m'  indicated  for  this  material  in  Table  1  of  [3). 

Finally,  we  recall  that  G(L)  does  not  depend  on  crack  length,  hence  type  L 
cracks  nu^h  be  arrested  rather  easily  due  to  local  reduction  in  applied  stress  or 
ply  toughness  Gu{L).  Figure  2  in  [5]  shows  photographs  of  crack  patterns  in  plies 
of  various  thicknesses.  An  examination  of  these  results  reveals  many  arrested 
cracks  in  a  thin  ply  (0.7S  mm)  where  htilure  depends  on  type  L  cracldng.  Only 
very  few  arrested  cracks  can  be  seen  in  thicker  plies  which  start  to  fail  by  type 
T  cracks,  and  where  G{L),  although  still  constant,  should  be  much  larger  than 
Gc(L)  when  the  type  L  crack  propagates  along  a  wide  front  b  =  a  >  b^i.L). 


la  DISCUSSION 

The  theoretical  results  can  predict  strength  of  plies  in  a  laminate  as  a  function 
of  ply  thickness  and  applied  stress  state,  provided  experimental  values  of  ply 
toughnesses  in  (13)  are  known.  Also,  the  correct  form  of  (12),  and  the  ratios  g{T) 
and  g{L)  in  (25)  and  (26)  are  needed  in  case  of  combined  loading.  While  some 
of  these  parameters  may  not  be  readily  available  in  certain  applications,  the 
results  are  of  value  in  estimating  the  benefit  that  enhanced  ply  toughness  may 
have  on  strength. 

Furthermore,  the  results  suggest  that  strength  of  thick  plies  may  be  sig¬ 
nificantly  reduced  by  damage  which  introduces  crack  nuclei  of  sufficient  size  into 
the  ply.  Free  edge  damage  caused  in  the  process  of  delamination  or  machining, 
as  well  as  impact  damage  or  penetrations  can  be  sources  of  such  initial  cracks. 
For  plies  which  have  been  damaged  in  this  manner,  one  should  assume  that  crack 
nuclei  of  width  26  «  2a  are  present  from  the  outset,  and  that  all  ply  cracking 
takes  place  in  th'  fiber  direction  Xj.  The  strength  of  such  damaged  plies, 
regardless  of  the’  thickness,  is  given  by  thin  ply  formulae  (46)  and  their  equiva¬ 
lents,  and  can  f  substantially  lower  than  that  of  undamaged  thick  plies. 

Another  application  of  the  present  results  can  be  made  if  the  growth  rate  of 
the  crack  nucleus  6  in  (14)  is  known  for  specific  loading  histories.  As  long  as  the 
material  properties  in  (13)  are  given  constants,  one  can  calculate  6c(r),  and  from 
that  the  effect  of  the  loading  history  on  ply  strength. 
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ABSTRACT 

A  summary  of  recent  results  Is  presented  on  the 
subject  of  progressive  ply  cracking  In  fibrous  lamin¬ 
ates.  First,  evaluation  of  stiffness  changes  caused  by 
systems  of  aligned  silt  cracks  which  are  parallel  to 
the  fiber  direction  In  a  unidirectional  composite 
lamina  Is  discussed.  Results  obtained  by  Che  self- 
consistent  method  are  presented.  Next,  a  procedure  for 
estimating  Instantaneous  crack  densities  and  stiffness 
changes  In  a  lamina  subjected  Co  prescribed  strain  his¬ 
tory  Is  outlined.  Specific  examples  and  comparisons  of 
analytical  and  experimental  results  are  presented  for 
two  graphite-epoxy  systems. 

INTRODUCTION 

In  many  fibrous  composite  systems  the  failure 
strain  of  the  fiber  far  exceeds  chat  of  the  matrix. 
Under  load,  the  difference  Is  usually  accommodated  by 
matrix  cracking.  This  is  frequently  observed  In  raono- 
tonically  or  cyclically  loaded  laminated  places,  where 
each  layer  may  contain  a  system  of  aligned  slit  cracks 
which  grow  in  the  direction  of  the  fibers  and  across 
the  thickness  of  the  ply.  Such  cracks  are  often  called 
ply  cracks  or  transverse  cracks,  although  It  is  more 
appropriate  to  reserve  the  latter  for  cracks  which  are 
perpendicular  to  the  fiber  axis,  and  refer  to  cracks 
which  grow  parallel  to  the  fiber  direction  as  axial 
cracks. 

In  polymer  matrix  composites  axial  matrix  cracking 
typically  starts  at  low  strain  levels  In  the  weakest 
off-axis  ply.  As  loading  continues,  cracks  appear  In 
other  off-axis  plies;  also,  their  density  Increases 
until  It  reaches  a  certain  saturation  level.  For 
example.  In  statically  loaded  (0/90)^  graphite-epoxy 
laminates  the  minimum  crack  spacing  was  observed  to  be 
equal  to  3. 5-4.0  thicknesses  [1].  In  metal  matrix  com¬ 
posites,  matrix  cracking  appears  to  be  caused  only  by 
cyclic  loads  which  exceed  the  shakedown  limit  of  the 
laminate  (21.  Under  such  circumstances  the  matrix 
experiences  cyclic  plastic  straining  and,  consequently, 
low-cycle  fatigue  failure.  Both  axial  and  transverse 


cracking  Is  present,  the  former  In  off-axis  plies,  the 
latter  In  zero  degree  plies.  The  crack  patterns  and 
densities  are  generally  similar  to  those  found  In  poly¬ 
mer  matrix  systems.  However,  saturation  density 
Increases  with  load  amplitude,  and  It  Is  not  unusual  to 
find  cracks  as  close  as  one  ply  thickness. 

In  polymer  matrix  systems  each  crack  apparently 
propagates  very  rapidly,  as  a  Griffith  crack.  In  metal 
matrix  systems,  under  cyclic  loading  the  cracks  grow 
slowly.  However,  each  single  crack  creates  only  a 
small  Increment  In  average  crack  density  of  the  ply. 
Therefore,  transverse  cracking  Is  progressive  In  the 
sense  that  crack  density  Increases  gradually  with 
applied  load,  regardless  of  the  rate  of  local  crack 
progagatlon. 

In  a  typical  part  of  a  laminated  composite  struc¬ 
ture,  removed  from  concentrated  loads  and  free  edges, 
matrix  cracking  Is  the  Initial,  low-stress  damage  mode 
under  applied  load.  It  is  eventually  followed  by  other 
types  of  damage,  such  as  de lamination  between  layers 
and  fiber  breaks;  but  these  appear  at  relatively  high 
loads  which  may  exceed  allowable  design  magnitudes.  In 
contrast,  matrix  cracks  grow  at  low  loads,  and  they  can 
significantly  Impair  stiffness  and  strength  of  com¬ 
posite  structures,  especially  those  containing  many 
off-axis  plies.  For  example,  fatigue  tests  on  both 
polymer  and  metal  matrix  laminated  plates  Indicate  that 
stiffness  and  residual  static  strength  reductions 
caused  by  cracking  In  plies  may  be  as  high  as  lO-SOX 
after  2  x  10^  cycles  of  loading  [2,3).  It  Is  therefore 
desirable  to  consider  the  effect  of  matrix  cracking  on 
composite  properties  In  design.  Furthermore,  matrix 
cracking  Is  unavoidable,  whereas  certain  other  damage 
modes,  such  as  free  edge  de lamination  can  be  prevented 
If  appropriate  precautions  are  taken  In  design. 

Sufficiently  general  theoretical  models  of  pro¬ 
gressive  cracking  In  composite  laminates  are  apparently 
not  to  be  found  In  the  literature.  Such  results  as  are 
available  for  angle-ply  laminates  have  been  obtained 
from  finite  element  calculations  [1),  while  other 
Studies  have  focused  on  (0/90,,)g  laminates  (4). 
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The  purpose  of  this  research  Is  to  develop  a  pro¬ 
cedure  for  prediction  of  crack  densities  In  Individual 
plies  of  a  laminated  composite  structure  as  a  function 
of  applied  load,  and  to  evaluate  the  effect  of  cracks 
on  stiffness  .of  the  structure. 

The  analysis  has  been  performed  In  the  following 
steps: 

a)  Evaluation  of  overall  thermoelastic  properties 

of  a  fibrous  composite  which  contains  a  cer¬ 
tain  density  of  aligned  slit  cracks. 

b)  Evaluation  of  crack  densities  and  stiffness 

changes  In  a  single  ply  which  Is  strained  In  a 
prescribed  way. 

c)  Evaluation  of  crack  densities  and  stiffness 

changes  In  a  laminated  composite  plate  which 
Is  subjected  to  prescribed  loading.  Specific 
examples  were  solved  for  laminated  plates 
under  In-pl  .ne  loads. 

STIFFNESS  CHANGES  CAUSED  BY  A  SYSTEM  OF  CRACKS  IN  A  PLY 

First  we  are  concerned  with  evaluation  of  overall 
compliance,  thermal  expansion  coefficients  (or  thermal 
strain  and  stress  vectors),  and  strain  and  stress  aver¬ 
ages  In  the  phases  of  a  fibrous  lamina  which  contains 
aligned  silt  cracks,  and  Is  subjected  to  uniform  mech¬ 
anical  loads  and  thermal  changes.  The  approach  to  the 
problem  was  outlined  In  reference  [5],  where  we  sug¬ 
gested  that  the  effect  of  matrix  crack  systems  on  pro¬ 
perties  of  fibrous  composites  can  be  analyzed.  In  prin¬ 
ciple,  by  the  same  techniques  which  are  commonly  used 
In  evaluation  of  elastic  constants  of  composite  mater¬ 
ials  and  fibrous  laminates,  e.g.,  by  the  self- 
consistent  method. 


TWO-PHASE  MODEL  OF  A  CRACKED  LAMINA 


*3 


Figure  1. 

lA.  An  infinite  fibrous  medium  with 
aligned  slit  cracks, 

IB.  A  fiber  lamina  with  parallel 
slit  cracks. 


THREE-PHASE  MODEL  OF  A  CRACKED  LAMINA 


Figure  2. 

2k.  An  infinite  medium  with  aligned 
fibers  and  slit  cracks, 

2B.  A  fiber  monolayer  with  cracks. 

The  essential  approximation  In  the  evaluation  of 
stiffness  and  compliance  changes  of  laminates  consists 
in  the  replacement  of  a  cracked  layer.  Figure  lb,  by  an 
effective  medium  which  contains  many  cracks.  Figure  la. 
The  crack  densities  can  be  exactly  matched  to  provide 
Identical  stiffnesses.  However,  the  cracks  in  the  layer 
are  not  entirely  surrounded  by  the  layer  material. 
Instead  they  interact  with  neighboring  layers  which  have 
different  elastic  properties.  This  Interaction  Is 
limited  to  the  vicinity  of  the  crack  tip,  thus  It  may  be 
Important  in  analysis  of  local  crack  growth  at  the 
Interface,  but  It  has  only  a  minor  effect  on  lamina 
stiffness.  We  note  that  a  similar  approximation  Is  com¬ 
monly  accepted  In  evaluation  of  elastic  moduli  of  lam¬ 
inated  composite  structures  reinforced  by  monolayers  of 
large  diameter  fibers. 

When  the  composite  Is  reinforced  by  monolayers  of 
fibers,  such  as  boron  or  silicon  carbide,  the  cracks  and 
fibers  may  be  of  similar  size.  The  appropriate  model  Is 
shown  In  Figure  2.  It  Is  analogous  to  that  of  Figure  1, 
except  that  It  contains  three  phases  (fiber,  matrix,  and 
cracks),  whereas  the  model  of  Figure  1  can  be  reduced  to 
two  phases  (composite  "matrix",  and  cracks)  Accord¬ 
ingly,  the  models  In  Figures  1  and  2  are  relocred  to  a 
two-phase  and  three-phase  models,  respectively.  From 
the  practical  standpoint,  the  effect  of  model  choice  on 
composite  stiffness  Is  small.  The  simpler  two-phase 
model  Is  thus  sufficient  for  analysis  of  all  fibrous 
systems . 

In  the  three  phase  model  we  designated  the  fiber  as 
phase  1,  matrix  as  phase  2,  and  cracks  as  phase  3.  In 
the  two  phase  model  phases  1  and  2  are  Joined  in  a  homo¬ 
geneous  "matrix"  and  designated  as  phase  2;  the  voids  or 
cracks  remain  phase  3. 
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The  self-consistent  analysis  of  the  cracked  con- 
poslte,  as  outlined  In  [5-7],  starts  with  a  composite 
geometry  In  which  the  cracks  are  Initially  regarded  aa 
ellltlcal  cylindrical  Inclusions*  A  self-consistent 
estimate  of  overall  stiffness  of  this  medium  la 
obtained.  Next,  the  Inclusions  are  evacuated,  l.e., 
replaced  by  voids.  Finally,  the  aspect  ratio  of  the 
voids  Is  adjusted  so  that  In  the  limit  the  voids  change 
to  cracks. 

To  evaluate  crack  density,  the  cracks  In  Figure  1 
are  first  replaced  by  elliptic  cylindrical  voids,  with 
a,  b,  denoting  the  major  and  minor  semiaxes.  If  n  Is 
the  number  of  voids  per  unit  area  In  the  xjx2 -plane, 
then  Che  volume  fraction  of  voids  Is  equal  to  C3  »  nabn, 
and  C2  +  C3  “  1.  Next,  let  the  voids  be  reduced  to 
cracks,  i.e.,  6  •  b/a  ♦  0.  Then 


which  are  expressed  In  terms  of  compliances  M  of  the 
effective  medium  as: 


-  “22“33"’^3  .  >J  »l 

22  rj -  (a  +  o  ) 

M33  1  2 


-  («44  »55^' 


(a]  +  ap 


2 

“  xa  n5 


where  6  ”  ^na^  Is  Che  crack  density  parameter.  In  fact, 
$  Is  equal  Co  the  number  of  cracks  of  fixed  length, 
e.g.,  2a,  In  a  square  of  side  2a.  For  example,  If  the 
cracks  are  located  In  a  ply.  Figure  lb,  then  6  measures 
the  distance  between  regularly  spaced  cracks  In  terms  of 
ply  thickness  2a.  At  6  ”  1,  the  distance  between  cracks 
Is  equal  to  2a,  as  6  decreases  the  distance  Increases 
and  at  8  "0  the  cracks  vanish.  We  recall  that  the 
observed  minimum  distances  between  cracks  In  a  satur¬ 
ation  state  are  equal  to  3. 5-4.0  ply  thicknesses,  l.e., 
7a-8a  [1].  Corresponding  values  of  6  are  0.28-0.25,  but 
values  as  high  as  8  ”  1  were  observed  in  the  B-Al  system 
12).  Therefore  0  £  B  £  1  is  the  appropriate  range  of 

On  the  macroscale,  Che  cracked  unidirectional  com¬ 
posite  of  Figure  1  can  be  regarded  as  an  orthotropic 
homogeneous  solid.  The  elastic  properties  of  the 
"matrix"  are  identical  with  those  of  the  fibrous  com¬ 
posite  and  can  be  easily  evaluated.  When  cracks  are 
Introduced,  the  macroscopic  or  overall  elastic  modujl 
of  the  solid  change.  To  make  Che  concept  of  overall 
moduli  meaningful,  It  Is  necessary  to  consider  overall 
uniform  loading.  Thus,  we  Introduce  uniform  overall 
average  stresses  ?  and  strains  C",  with  components 
arranged  in  (6x1)  column  vectors  and  related  by  consti¬ 
tutive  equations^ 


where  aj  and  02  are  the  roots  of 


(M22M33-M23)a2  -  {M33Mg6  +  2(Mi2M33  -  Mi3M23))a 


+MJ1M33  -  M]3  -  0  (7) 

These  results  Imply  Chat  only  three  compliance 
coefficients  M22 »  M44,  and  M66  are  affected  by  the 
Introduction  of  cracks,  the  remaining  six  terms  In  M 
are  unchanged,  l.e.,  they  remain  equal  to  those  of  the 
uncracked  fiber  composite.  In  particular 
000 
”11  *  »11  •  ”12  ’  ”12  .  M13  •  «13 


000 
M23  -  M23  ,  M33  -  M33  ,  M55  »  M55 


M22  -  M22  6(M22M33  -  <8) 

M44  «  M44  +  B(M44M55)' (9) 


”b6  “  B(M22M33-M23)  (M11M33-M13) 

,  1/2  .  l/2^,u 
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ai.d  from  (7): 


where  L,  and  M  are  the  overall  stiffness,  and  compli¬ 
ance  (6x6)  matrices  of  the  cracked  composite,  respec¬ 
tively.  M  =  when  the  inverse  exists.  Effectl(re 
properties  of  uncracked  fibrous  material  (phase  2)  are 
denoted  by  L°,  M°,  or  by  L2 ,  M2. 

For  a  dilute  concentration  of  aligned  slit  -racks, 
the  matrices  L  and  M  can  be  evaluated  by  the  self- 
consistent  method.  This  was  done  In  Reference  (5);  the 
results  are: 

L  =  LO  -  r  1°  A  h  (3) 
M  =  MO  *•  fT  A  ,  (4) 


2  2 
”  ^”i1”33-Mi3)/(M22M33-M23) 


+  “2  ’  ^'^33^66  2(M]2M33-Mi3M23)1/ 


IM22M33-M23I  . 

The  unknown  shear  compliance  M44  can  be  (  btained 
directly  from  (9): 


”44  "  ”44  I'e  «55  ^  "”44  ”55>‘’J 
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The  matrix  A  has  only  Chr**e  nonzero  components. 


In  (5-7],  (6x6)  oiatrtces  are  denoted  by  capital 
Latin  or  Greek  letters,  e.g.,  L,  M,  A,  A,  B,  P,  Q,  and 
(6x1)  matrices  by  lower  case  bold  face  letters,  e.g., 

Q. 


The  remaining  unknowns  M22  and  are  found  from  (8) 
and  (10). 

These  results  can  be  utilized  In  (4)  for  a  more 
direct  evaluation  of  the  three  nonzero  components  of 
the  matrix  A  In  (6): 

^IJ  -  («1J  -  Mlj)/^  •  (iJ  ■  22.  66)  (12) 

Of  course,  Che  same  result  follows  from  (6)  and 
(7).  In  any  case,  the  may  now  be  substituted  Into 

(3),  and  Che  components  of  Che  overall  stiffness  L  can 
be  found  In  a  closed  form.  The  resulting  expressions 
are  given  In  reference  [7,  eqn.  (A-6)]. 

The  above  procedure  leads  also  to  expressions  for 
thermal  response  of  Che  cracked  lamina.  Stress  and 
strain  averages  In  Che  composite  matrix  (phase  2)  and 
cracks  (phase  3)  were  evaluated  as  well.  All  results 
were  found  In  closed  form,  as  functions  of  A  In  (6)  or 
(12).  Specific  forms  appear  In  (7]. 

It  Is  seen  that  the  results  are  remarkably  simple, 
and  similar  to  Chose  that  are  routinely  used  for  evalua¬ 
tion  of  elastic  moduli  of  an  uncracked  composite  medium. 
This  similarity  Is  particularly  useful  In  applications 
of  the  results  to  laminated  structures.  Existing 
theories  for  reinforced  materials  can  be  adapted  Co 
Incorporate  Che  effect  of  matrix  cracking.  Indeed,  one 
may  regard  cracks  as  negative  reinforcement  which 
reduces  stiffness  of  a  composite,  essentially  In  Che 
same  way  as  fiber  reinforcement  enhances  It. 

PROGRESSIVE  MATRIX  CRACKISG  IN  A  PLY 

This  subject  Is  discussed  extensively  In  (8).  Ac¬ 
cumulation  of  transverse  cracks  In  a  uniformly  strained 
fibrous  lamina  which  Is  a  part  of  a  laminated  plate  Is 
analyzed  there  by  two  methods.  In  the  first  one, 
thickness  averages  of  strains  and  stresses  In  unbroken 
segments  of  the  cracked  lamina  are  evaluated  In  an 
approximate  way  and  combined  with  fracture  criteria  for 
extension  of  Initial  flaws,  and  for  propagation  of  silt 
cracks  In  Che  lamina.  It  Is  found  chat  this  approach 
cannot  be  applied  In  practice  at  this  time  oecause  of 
Insufficient  Information  about  the  residual  strength  of 
the  uncracked  segments.  An  alternative  method  Is  pro¬ 
posed  which  regards  the  cracked  lamina  as  an  effective 
medium.  Volume  averages  of  stresses  and  strains  are 
estimated  by  a  self-consistent  procedure,  and  combined 
with  effective  failure  criteria.  This  approach  can  be 
applied  with  minimum  of  experimental  Information  about 
ply  strength,  such  as  first  ply  failure  data.  The 
effective  failure  criterion  Is  satisfied  by  adjustment 
of  volume  stress  averages  through  Increments  of  crack 
density,  and  the  resulting  ply  stiffness  reductions, 
under  Increasing  strain.  Crack  densities  and  stiffness 
changes  may  be  overestimated  In  this  case,  but  the  re¬ 
sults  remain  on  the  safe  side  In  structural  applica¬ 
tions.  Specific  results  were  found  for  two  Gr/Ep  sys¬ 
tems  . 

The  analysis  of  progressive  cracking  In  an  effec¬ 
tive  medium  Is  made  as  follows:  Suppose  that  a  com¬ 

posite  ply  Is  subjected  to  current  uniform  strain  t, 
and  that  the  corresponding  crack  density  Is  Bq  ^  0.  A 
strain  Increment  dc,  and  a  thermal  change  d0  are  ap¬ 
plied,  and  a  new  magnitude  of  B  Is  sought.  B  Is  regard¬ 
ed  as  a  continuous  function  of  strain,  c.f.,  (7),  Sec¬ 

tion  2.  The  overall  constitutive  equation  of  the  ply 
Is  given  by  (2). 

The  solution  can  be  obtained  In  the  following 


ut  f’  -  c  +  dc  ,  0'  -  e  +  de  (13) 

be  the  new  applied  strain  and  thermal  change.  The 
corresponding  average  stress  In  the  ply  Is  In  analogy 
with  (2 ) 

a*  -  f  ?•  +  6-1*  . 

where  L'  3  L'(B),  t'  “  t'(8)  are  material  properties 

from  (4)  and  (21)  in  [7)  at  as  yet  unknown  value  of 

B  -  80  +  dB  :  dB  >  0  .  (15) 


The  stress  components  of  Interest  In  the  effective 
failure  criteria  below  are,  according  to  equations  (21) 
and  (A-3)  to  (A-6)  In  (7); 


o*  -  L*  (E*  -  a  0’)  +  f  (e'  -  a  0')  +  L’  (£'  -oO’) 

22  12  11  T  -  22  22  T  23  33  A 


o'  -  2L'  e' 
23  44  23 


The  effective  failure  criteria  are  selected  In  the  form 


la*  \ 

22 
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23\ 

l- 

In  1 

Or,  In  analogy  with  the  Tsal-Wu  form: 

__  __2__2 

F  o'  +  F  (o'  )  +  F  (o’  )  -  1  (18) 

2  22  22  22  66  23 

For  B  >0  the  magnitudes  of  overall  strength 

parameters  ^2.  ^22.  ^66  ®Ee  assumed  to  be 

known  functions  of  crack  density  B  and  ply  thickness 
2a.  Their  values  must  be  adjusted  after  each 
Increment  In  B.  If  these  parameters  are  known  only  In 
terms  of  first  ply  failure,  then  the  Lj^j  reduction 
alone  mist  balance  the  stresses  In  (17)  or  (18).  The 
crack  density  and  stiffness  reduction  are 
overestimated  In  such  a  case,  but  model  predictions 
remain  on  safe  side  In  structural  applications. 

To  determine  If  d£  and  dg  actually  cause 

additional  cracking.  It  Is  first  necessary  to  evaluate 
the  stresses  (16)  with  L'  ■  l.(6o)  and  t'”t(6o).  l.e., 
with  the  material  properties  at  the  original  crack 
density  So-  If  these  stresses  make  the  left-hand  side 
of  (17)  or  (18)  larger  than  unity,  then  an  Increment 
dB  Is  required  to  reduce  the  coefficients  L'  ,  and 

IJ 

thus  the  stresses  to  the  level  allowed  by  the  particular 
failure  criterion. 

We  recall  that  both  L  and  t  are  decreasing  func¬ 
tions  of  B,  thus  the  ply  softens  under  Increasing  strain 
but  remains  elastic  during  unloading  and  reloading  to 
the  current  stress  level  (17)  or  (18). 

To  Illustrate  typical  results  provided  by  (17), 
consider  a  graphite-epoxy  ply  with  properties  given  In 
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Table  I.  Theoretical  predictions  of  relative 
changes  in  axial  (El)  and  trans¬ 
verse  (Et)  Young's  moduli  of  lamin¬ 
ated  plates  at  given  values  of  crack 
density  6  in  all  plies. 


table  II  in  I?].  The  ply  Is  deformed  by  prescribed 
strains  C22  and  eij  In  several  radial  directions,  while 
stresses  are  adjusted  so  that  on  -  012  -  013  ’  033  * 

For  greater  clarity  of  presentation,  the  ply  strengths 
are  taken  as  constant: 

Oj  -  o"  -  28. 4A  MPa  ,  Ojj  -  a"  -  37.95  MPa 

These  values  were  measured  for  first  ply  failure  °f 
unconstrained  ply.  hence  they  represent  relatively  low 
strengths. 


The  results  appear  In  Figure  3.  In  each  case,  the 
applied  strain  is  plotted  on  the  horizontal  axis,  the 
average  stresses  and  crack  densities  on  the  vertical. 
The  constant  sLrengths  levels  are  plotted  as  well.  It 
Is  easy  to  see  that  If  strengths  were  increasing  with  g, 
then  the  stress  components  would  be  elevated  accord¬ 
ingly, 

A  remarkable  feature  of  these  results  Is  the 
strong  effect  of  shear  strain  E23  on  crack  density 
increase.  When  only  shear  ^23,  or  only  transverse 
normal  strain  '€22  Is  applied,  one  finds  that  3  increases 
about  twice  as  fast  In  shear  when  2623  is  compared  in 
magnitude  toT22,  At  the  same  time,  the  shear  strength 
exceeds  the  transverse  tension  strength  by  33%.  This 
behavior  can  be  related  to  the  rate  of  change  of  the 
relevant  stlffnessed  1,22  end  1.4/,  with  3,  We  recall  from 
Figure  2  in  (5]  that  1.22  decreases  with  3  much  more 
rapidly  than  L44,  from  comparable  initial  levels.  ^Thus 
023  fs  relatively  large2_  much  larger^  than  022  for  2e23  > 
r22.  Even  at  low  E23/E22  ratios  023  keeps  Increasing 
with  applied  strain  while  022  always  decreases  after 
onset  of  cracking. 

In  addition  to  stress  changes,  the  ply  experiences 
stiffness  change  under  deformation  conditions  of  Figure 
3.  These  are  summarized  in  Figure  A. 

PROGRESSIVE  CRACKING  IM  LAMINATED  PLATES 

The  above  procedure  for  analysis  of  crack  accumul¬ 
ation  in  a  strained  ply  can  be  combined  with  laminated 
plate  theory  and  thus  extended  to  cracking  laminated 
plates.  The  theoretical  background  has  been  fully 
developed,  and  several  examples  have  been  solved.  A  com¬ 
plete  'description  of  the  results  will  appear  in 
reference  19).  Figure  5  illustrates  stress  and  stiff¬ 
ness  changes  In  a  plate  subjected  to  combined  loads. 

As  In  the  case  of  a  lamina,  cracking  in  laminated 
plates  must  be  analyzed  In  an  Incremental  way.  In  many 
applications  one  cannot  predict  the  exact  loading  his¬ 
tory,  or  a  detailed  analysis  may  be  undesirably  complex. 
Under  such  circumstances.  It  is  sufficient  to  estimate 
the  maximum  stiffness  loss  caused  by  the  greatest  pos¬ 
sible  damage  In  individual  plies.  Ply  properties  for 
this  case  were  found  In  [7],  for  large  values  of  3. 
Actual  ply  stiffness  values  In  this  case  are  well 
approximated  Ldien  3  Is  taken  as  equal  to  unity. 

The  effect  of  a  fixed,  uniform  crack  density  In  all 
plies  of  the  laminated  plate,  on  longitudinal  (El)  and 
transverse  (E^.)  in-plane  moduli  of  plates  of  different 
layups  Is  Illustrated  In  Tables  I,  II,  and  III  [10],  It 
Is  seen  that  the  stiffness  loss  caused  by  extensive 
cracking  can  be  very  significant  In  certain  composite 
systems, 

A  comparison  of  the  theoretical  results  with  exper¬ 
imental  observations  Is  presented  In  Figure  6.  Axial 
stiffness  of  a  (0,903)^  E-glass-epoxy  plate  Is  plotted 
against  observed  crack  density  for  a  uniaxial  tension 
test.  The  experimental  points  and  shear  lag  analysts 
were  presented  In  [111*  The  theoretical  (SCM)  curve  was 
calculated  with  the  self-consistent  approximation 
described  earlier.  A  very  good  agreement  of  this  curve 
with  experiments  Is  Indicated.  A  remarkable  feature  of 
this  correlation  Is  that  the  only  material  constants 
used  In  the  analysis  were  the  original  elastic  moduli  of 
Che  uncracked  composite.  With  this  information  the 
stiffness  changes  were  evaluated  for  given  crack  density 
from  (3)  for  the  90°  piles,  and  overall  laminate  stiff¬ 
ness  change  was  found  from  laminated  plate  theory. 
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Figure  5.  Calculated  changes  in  ply  stresses 
3',-j  and  plate  stiffnesses  •  in  a 
graphite-epoxy  plate  loaded  In  simple 
tension  a-g  in  the  0°  direction.  All 
stresses  are  in  local  ply  coordinates, 
the  fiber  direction  coincides  with 
local  xg  axis  of  a  ply. 


Additional  comparisons  were  made  with  experimental 
results  obtained  by  Hwang  and  Hahn  [12]  who  measured 
crack  density  as  a  function  of  applied  stress  or  strain 
In  AS/3501-5A  laminated  plate  coupons  In  tension,  and, 
with  similar  experiments  reported  by  Ryder  and  Crossman 
[13]  for  laminated  T300/5208  plate  coupons. 

Analytical  predictions  of  crack  densities  In  these 
laminated  plate  coupons  were  made  as  follows:  First, 
laminated  plate  theory  was  used  to  calculate  residual 
thermal  stresses  In  each  ply  caused  by  cooldown  from  the 
curing  temperature.  The  temperature  difference  was 
AT»147°C  for  both  systems.  Then  ply  stresses  caused  by 
a  uniaxial  applied  load  were  evaluated  and  superimposed 
on  the  thermal  residual  stresses  In  each  ply.  Next, 
(18)  was  rewritten  for  each  ply  In  the  form  given  by 
Hahn  [14]: 

(1-g)  (^)  +g  (^)  +  &  -1  =  0  (19) 

OJ  3-1  011 

where  the  stresses  022 »  023  sre  taken  In  local  ply  coor¬ 
dinates,  X3  is  parallel  to  local  fiber  direction.  The 
CJi,  Ojj  denote  critical  stresses  at  first  ply  failure. 
The  parameter,  which  accounts  for  failure  mode  inter¬ 
action  In  an  approximate  way  was  taken  as  g-0.1. 

Now,  to  determine  0i  and  On  for  the  laminate  test 
data  shown  In  Figures  7  and  8,  one  calculates  the  local 
stress  Ooy  at  first  ply  failure  In  the  90“  ply.  Since 
023  “  0  In  this  ply,  one  obtains  Oj  from  (19).  Then, 
additional  stress  Is  applied  to  the  laminate,  and  crack 
density  B  In  the  90“  ply  Is  calculated  Incrementally, 
according  to  the  sequence  (13)  to  (19),  The  stiffness  of 
the  90“  ply  gradually  decreases  In  the  process.  When 
the  applied  stress  reaches  the  magnitude  required  for 
first  ply  failure  in  the  +45“  ply,  in  this  ply  Is 

again  found  from  (19). 

These  results  are  listed  In  Table  IV  as  predictions 
(A)  and  (B)  for  the  two  coupons. 

In  a  similar  way,  the  results  In  Figures  9  and  10 
were  interpreted  and  provided  values  of  oj  ando'Tiln  pre" 
dictions  (C),  (D)  and  (E),  Table  IV. 


Figure  6.  Comparison  of  experimental ly  observed 
stiffness  changes  and  crack  densities 
[11]  with  theoretical  predictions. 

The  (SCM)  curve  was  calculated  from 
the  self-consistent  estimate  of  stiff¬ 
ness  in  90“  plies,  Equation  (3). 
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Figure  7.  Comparison  of  predicted  and  observed 
average  crack  spacing  in  a  laminated 
plate  subjected  to  uniaxial  tension. 
Experimental  data  from  Hwang  and 
Hahn  (1982). 
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Figure  8.  Comparison  of  predicted  and  observed 
average  cruck  spacing  in  a  laminated 
plate  subjected  to  uniaxial  tension. 
Experimental  data  from  Hwang  and 
Hahn  (1982). 


Strain  (KT®) 

Figure  9.  Comparison  of  predicted  and  observed 
average  crack  spacing  in  a  laminated 
plate  subjected  to  uniaxial  exten¬ 
sion.  Experimental  data  from  Ryder 
and  Crossman  (1984). 


Table  IV.  Effective  Ply  Strength 
At  First  Ply  Failure 
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Figure  10.  Comparison  of  predicted  and  observed 
average  crack  spacing  in  a  laminated 
plate  subjected  to  uniaxial  exten¬ 
sion.  Experimental  data  from  Ryder 
and  Crossman  (1984). 


The  experloentally  found  values  of  the  critical 
stress  and  at  first  ply  failure  in  (19)  are 
actually  Che  sole  material  parameters  needed  to  fit  the 
data.  The  only  additional  information  entering  the 
analysis  are  original  elastic  moduli  of  Che  undamaged 
laminate,  and  the  empirical  factor  g“0.1  in  (19). 

The  magltudes  of  and  oii  are  fairly  consistent. 
Note  that  predictions  (A)  and  (B)  apply  to  both  0“  and 
45°  plies  in  Figure  ^  and  8.  Differences  between  (A) 
and  (B)  are  caused  by  differences  in  experimental  data. 
Prediction  (C)  is  sufficient  for  all  plies  In  Figure  9, 
except  for  the  -45°  ply  as  explained  below.  Predic¬ 
tions  (D)  and  (E)  in  Figure  (10)  are  different  because 
the  data  points  for  the  two  coupons  are  so  far  apart. 

It  would  be  desirable  to  find  that  one  set  ot  crit¬ 
ical  stresses  fits  all  data  for  a  particular  system,  or 
for  the  two  grafite-epoxy  systems  of  Figures  7  to  10. 
This  is  not  exactly  possible  but  we  note  chat  in  (A) 
(B)  and  (D)  are  similar;  ideally  they  should  coincide 
since  they  refer  to  a  single  system.  Still  they  are 
closer  together  than  those  in  (D)  and  (E)  found  for  not 
only  one  system,  but  also  one  layup.  In  fact  the  two 
coupons  A1  and  A2  were  taken  from  the  same  laminated 
plate. 


Therefore,  one  can  conclude  that  the  differences  in 
effective  ply  strength  at  first  ply  failure  In  Table  IV 
are  caused  by  variations  in  material  properties  that 
affect  the  experimental  data. 

A  surprising  feature  of  the  experimental  results 
shown  in  Figures  7  and  9  is  that  no  cracks  were  observed 
in  the  interior  -45°  plies.  One  would  certainly  expect 
CO  find  cracks  in  these  layers  idilch  experience  the  same 
average  strain  state  as  the  -4-45°  plies  where  numerous 
cracks  were  found.  In  fact,  the  -45°  plies  form  a 
single  ply  of  double  thickness,  and  this  should  make  it 
more  susceptible  to  matrix  cracking  chan  the  +45*  plies 
of  single  thickness.  This  phenomenon  tequires  further 
study. 
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CONCLUSIONS 

!•  An  analytical  technique  has  been  developed  for 
modeling  of  progressive  transverse  matrix  cracking 
In  laminated  composite  plates.  The  analysis  con¬ 
sists  of  three  steps  .  First,  self-consistent 
estimates  of  laminate  stiffness  changes  are  found 
for  a  given  crack  density  from  (3)  and  (A).  Next, 
the  rate  of  crack  density  Increase  under  applied 
load  Is  evaluated  from  (13)  to  (19)  on  the  basis  of 
a  selected  ply  failure  criterion.  Finally,  the  ply 
analysis  Is  made  In  terms  of  ply  stress  state  for  a 
ply  embedded  In  a  laminate  with  several  cracking 
plies.  The  Instantaneous  ply  stresses  follow  from  a 
simple  adaptation  of  laminate  plate  theory. 

2.  A  good  agreement  of  the  theory  with  several  sets  of 
experimental  data  was  found.  The  theoretical  pre¬ 
dictions  require  a  minimum  amount  of  experimentally 
derived  Information,  such  as  first  ply  failure 
stresses  In  the  cracking  ply,  and  elastic  moduli  of 
the  undamaged  laminate.  No  empirical  parameters  are 
required  outside  the  ply  failure  criteria. 

3.  The  theory  makes  it  possible  to  calculate  the 
Instantaneous  crack  density  In  each  ply,  the  Instan¬ 
taneous  stresses  In  each  ply  of  the  laminate,  as 
well  as  the  average  fiber  and  matrix  stresses  In 
each  ply,  after  application  of  each  load  or  strain 
Increment  to  the  laminate.  The  laminate  can  be  sub¬ 
jected  to  combined  three  dlmentlonal  loading,  and 
also  to  varying  uniform  thermal  changes  In  the 
course  of  mechanical  loading. 

A.  The  analysis  can  be  Implemented  through  a  simple 
numerical  routine,  even  on  a  personal  computer.  In 
fact,  the  entire  procedure  for  each  loading  step  is 
rather  similar  to  that  used  In  elasticity  analysis 
of  laminated  plates. 
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ANALYSIS  OF  FIRST  PLY  FAILURE  IN  COMPOSITE 

LAMINATES 

GEORGE  1.  DVORAK 

Department  of  Civil  Engineering,  Rensselaer  Polytechnic  Institute,  Troy.  NY  12180.  U.S.A. 

and 

NORMAN  LAWS 

Department  of  Mechanical  Engineering,  University  of  Pittsburgh.  Pittsburgh,  PA  15260.  U.S.A. 

Abatnct — The  mechanics  of  transver^  cracking  in  an  elastic  fibrous  composite  ply  is  explored  for  the 
case  of  low  crack  density.  Cracks  are  assumed  to  initiate  from  a  nucleus  created  by  localized  fiber 
debonding  and  matrix  cracking.  It  is  found  that  cracks  may  propagate  in  two  directions  on  planes 
which  are  parallel  to  the  fiber  axis  and  perpendicular  to  the  tnidplane  of  the  ply.  In  general,  crack 
propagation  in  the  direction  of  the  fiber  axis  controls  the  strength  of  thin  plies,  while  cracking  in  the 
direction  perpendicular  to  the  fiber  axis  determines  the  straigth  of  thick  plies.  The  theory  relates  ply 
thickness,  crack  geometry  and  ply  toughness  to  ply  strength.  It  predicts  a  significant  increase  in 
strength  with  decreasing  ply  thickness  in  constrained  thin  plies.  The  strength  of  thick  plies  is  found  to 
be  constant,  but  it  may  be  r^uced  by  preexisting  damage.  Results  are  illustrated  by  comparison  with 
experimental  data. 


1.  INTRODUCTION 

In  THIS  paper  we  discuss  some  fracture  mechanics  problems  associated  with  transverse  cracking  in 
composite  laminates.  We  focus  our  attention  at  a  single  ply  in  a  typical  part  of  a  laminated  structure, 
removed  from  free  edges  and  other  stress  concentrations.  The  ply  is  subjected  to  a  certain 
deformation  history  in  the  course  of  incremental  loading  applied  to  the  laminate.  In  particular,  we 
explore  the  initial  stage  of  the  cracking  process  which  consists  of  first  ply  failure,  and  of  few 
subsequent  failures  at  locations  which  are  far  apart  from  each  other,  so  that  the  cracks  do  not  interact. 
Our  goal  is  to  identify  the  failure  modes  which  cause  transverse  cracking,  relate  them  to  appropriate 
fracture  problems,  and  analyse  the  different  stages  of  the  cracking  process. 


2.  CRACK  GEOMETRY 

A  unidirectionally  reinforced  ply  of  thickness  2a  is  embedded  in  a  laminated  plate  or  a  similar 
composite  structure.  The  structure  is  subjected  to  a  certain  incremental  loading  program  which,  at 
each  loading  step,  causes  a  known  instantaneous  strain  state  in  the  ply.  We  assume  that  in  a  crack- 
free  ply  the  local  strains  are  uniform,  or  nearly  uniform,  and  equal  to  their  volume  averages,  or  overall 
strains  i.  The  strain  -tate  i,  the  thermal  change  0,  and  the  resulting  overall  stress  a  in  the  ply  are 
related  by: 


a  =  Lo«-^o.  i  =  Moff-l-0mo  (1) 

where  Lq,  are  stiffness  and  compliance  matrices  of  the  uncracked  ply,  and  Iq.  mo  are  the  thermal 
stress  and  strain  vectors.  We  emphasize  that  a,i  are  ply  averages,  not  laminate  averages. 

In  many  plate  composite  structures,  each  ply  usually  supports  a  plane  state  of  stress.  In  the 
coordinates  of  Fig.  1  this  suggests  that  overall  stress  components  ^22,  ^23  ^tid  <733  are  large  in 
comparison  with  dn,  ffu  and  ^,3.  It  is  easy  to  see  that  stresses  1x22  and  du  favor  cracking  on  x,X3- 
planes.  The  stress  ^33  does  not  affect  these  cracks.  The  remaining  minor  shear  components  d,2  and 
(T,  3  may  exist  in  bent  plates,  but  only  2  may  contribute  to  cracking.  The  normal  stress  d ,  1  is  usually 
small.  Under  such  loading  conditions  the  ply  will  crack  as  indicated  in  Fig.  1 .  The  laminate  layers  on 
either  side  of  the  ply  have  a  certain  thickness  b  and  b'.  Kb  ^  b',  then  we  require  that  b  >  2a,  the  ply 
thickness. 


Fig.  1.  Slit  crack  in  a  uniformly  strained  composite  ply.  Fibers  and  crack  are  aligned  with  Xj  axis. 


3.  CRACK  NUCLEATION 

It  is  commonly  held  that  transverse  cracks  in  polymer  matrix  composites  are  associated  with 
initial  fiber  debonding.  Several  debonds  may  coalesce  into  a  stable  crack  nucleus  which  grows  slowly 
until  it  reaches  a  critical  size,  and  then  propagates  as  a  Griffith  crack.  Since  the  in-plane  contour  of 
each  fiber  debond  resembles  an  ellipse  or  an  elongated  oval,  it  seems  reasonable  to  assume  the  same 
contour  for  the  initial  microcrack.  The  resulting  shape  appears  in  Fig.  2.  Let  28  denote  the  width  of 
the  microcrack.  The  length  in  Xj,  which  may  be  much  larger  than  28,  need  not  be  specified.  The 
thickness  of  the  ply  is  2a,  the  microcrack  is  located  at  mid-thickness  but  this  docs  not  affect  the  results 
that  follow. 

Now,  if  the  ply  is  strained  in  an  incremental  way,  the  microcrack  will  grow  slowly  until  it 
becomes  a  crack  nucleus  of  certain  critical  width  28,.,  and  starts  to  propagate  as  a  Griffith  crack.  The 
slow  growth  process  is  probably  time-dependent  in  most  polymer  matrix  systems,  and  one  should 
relate  the  growth  rate  to  loading  history.  Therefore,  we  write 

8  =  5(i(t)),  for  8  ^ 

The  rapid  crack  propagation  can  be  regarded  as  inviscid,  at  least  in  brittle  resin  matrices. 

4.  STRENGTH  OF  THIN  PLIES 

We  define  thin  plies  by  the  relation 

28  =  2a,  for  8  ^  8„.  (3) 

In  other  words,  the  ply  thickness  is  such  that  the  debonded  region  remains  stable  while  it  grows  in  the 
ply,  and  it  may  become  unstable  only  after  it  has  extended  across  the  entire  ply  thickness.  The 
definition  clearly  depends  on  ply  properties,  such  as  toughness  and  crack  growth  rate  (2).  It  also 
suggests  that  plies  which  contain  through  the  thickness  initial  flaws  or  cracks  induced  by  other  types 
of  damage  must  be  regarded  as  “thin”  plies  regardless  of  their  actual  thickness. 


First  ply  failure 
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Fig.  3.  Extension  by  unit  length  of  a  type  L  slit  crack  in  a  fibrous  composite  ply. 


The  energy  release  rate  of  the  crack  nucleus  can  be  derived  according  to  Fig.  3.  Initially,  the 
nucleus  consists  of  a  prismatic  part  CC,  and  precursor  regions  which  are  approximated  as 
semielliptical  cracks.  The  crack  nucleus  can  extend  only  in  the  fiber  direction  Xj,  and,  as  it  does  so,  a 
unit  extension  CC'  will  increase  the  crack  area  by  2a*  1.  The  prismatic  part  extends  by  unit  length  and 
one  of  the  precursor  regions  translates. 

Let  denote  the  energy  released  by  formation  of  the  prismatic  section  of  unit  length  of  a  ply 
crack  of  width  2a.  This  is  the  interaction  energy  of  the  crack  section.  Then,  the  energy  release  rate  of 
the  crack  can  be  written  as: 

G(L)  =  W;/2a  (4) 

where  (L)  indicates  that  the  crack  extends  in  the  longitudinal  direction  in  the  ply.  As  in  [1],  we  refer  to 
such  cracks  as  longitudinal  or  type  L  cracks. 

It  is  probably  obvious  that  is  equal  to  the  energy  released  by  a  microcrack  of  unit  width  Xj 
which  has  extended  in  the  transverse  direction  x,  to  length  2a.  One  can  write: 

G(L)  =  - iGfTld^,  (5) 

a  Jo 

where  G(T)  is  the  energy  release  rate  for  extension  of  the  microcrack  in  the  transverse  direction  x, . 
Again,  cracks  growing  in  the  x,  direction  will  be  referred  to  as  T-type  cracks. 

If  the  crack  nucleus  was  located  in  a  very  thick  ply,  6  ^  a,  then  W,  could  be  determined  from 
known  solutions  for  cracks  in  orthotropic  solid$[2, 3].  In  a  thin  ply,  the  crack  interacts  with  adjacent 
plies  of  different  elastic  moduli,  and  this  has  an  influence  on  W,.  As  in  [1],  this  can  be  reflected  by 
writing  (4)  in  the  form 

G(L)  =  i7Ca[iJ|A22d22+^ll^*4^23  +  'J|II^66®12]- 

Here  d,j  are  the  overall  stresses  caused  in  the  ply  by  the  applied  strains  iij,  and  by  a  uniform  thermal 
change  i  =  I,  II,  III  are  reduction  coefficients  which  rdfiect  interaction  of  the  crack  in  Fig.  3  with 
the  adjacent  plies,  for  each  of  the  crack  opening  modes  I,  II  or  III.  Note  that  mode  I  is  the  opening 
mode  caused  by  ^22’  mode  II  is  the  longitudinal  shear  mode,  and  mode  III  is  the  transverse  shear 
mode.  The  A,y  are  non  vanishing  components  of  the  crack  tensor  A(° ,  derived  in  [4]  and  [5] ,  and  taken 
for  dilute  concentration  of  cracks: 


where  £l  =  £33,  Ej  =  =  £22  are  the  longitudinal  and  transverse  Young’s  moduli,  Gl  =  L44  = 

Lssin  [43  is  the  longitudinal  shear  modulus,  and  Vl  =  v,3  =  V23  is  the  corresponding  Poisson’s  ratio. 
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Coeflicients  if,  can  be  found  for  many  crack  configurations  in  layered  media,  this  is  discussed  in 
[1].  For  example,  (fi=0.65  in  (0/90)  graphite-epoxy,  and  (f,=0.8  in  (0/90)  glass-epoxy.  However, 
one  should  keep  in  mind  that  the  crack  tips  at  Xj  =  ±a  reside  in  a  soft  matrix  which  may  not  be  able 
to  support  the  high  crack-tip  stresses.  Therefore,  as  the  crack  propagates  in  the  longitudinal 
direction,  local  crack  extension  along  ply  interface,  and  across  this  interface  into  the  adjacent  ply  may 
take  place  to  relieve  the  stress  singularity.  Such  localized  crack  extension  has  been  observed 
experimentally;  it  may  be  present  as  delamination  as  well  as  extensive  fiber  cracks  in  the  next 
ply [6, 7].  This  additional  cracking  will  elevate  W^,&nd  the  reduction  coefficients  if,  in  (6).  Exact  values 
of  are  not  available  for  running  cracks. 

Fracture  strength  of  a  thin  ply  can  be  found  on  the  grounds  that  the  G(L)  reaches  a  certain 
critical  value.  If  dji  ot  ^23  were  acting  alone,  then  one  could  compare  G,(L)  or  G„(L)  to  the 
corresponding  ply  toughness  G|c(L)  or  G„c(L).  However,  under  combined  loading,  the  modes  are 
interrelated.  Specific  criteria  which  describe  this  effect  need  to  be  derived  from  experimental  data. 
Hahn[8]  offers  a  form  based  on  fracture  mechanics  considerations  that  fits  data  obtained  for  certain 
glass-epoxy  and  graphite-epoxy  systems: 


where  g  =  G,,/G„,. 

If  Gc  are  independent  of  loading  mode,  then  g-*  1  and  (8)  assumes  a  familiar  form.  In  actual 
composite  plies  with  epoxy  matrices  p  =  0.1.  Note  that  no  allowance  is  made  for  mode  III  cracking. 
Modes  I  and  II  in  a  composite  ply  always  create  asperities  on  the  fracture  surface,  as  the  crack  seeks 
its  way  between  densely  packed  fibers.  These  asperities  can  be  interlocked  and  thus  impede  mode  III 
cracking.  Also,  in  a  typical  part  of  laminated  plate  or  shell,  the  corresponding  stress  is  usually 
small  or  equal  to  zero.  Therefore,  mode  III  need  not  be  considered  in  most  applications. 

From  (6)  and  (8)  one  can  find  strength  of  thin  plies  as 

(^22)c,  =  [4G,.(L)/(7t^,AS^a)]>'^  at  ^^3  =  0,  (9a) 

(^23)cr  =  [4G||j(L)/(7tiJi|A24a)3*^^,  a.tff22  =  ^-  (9b) 

In  the  general  loading  case  ^33  =  one  can  calculate  {ff22)cT  fro™  (8)- 

Accordingly,  the  strength  of  thin  plies  becomes  related  only  to  onset  of  unstable  longitudinal  (type 
L)  cracking. 

If  experimentally  measured  strength  data  are  available  for  thin  plies  of  variable  thickness,  then 
one  can  find  G|c(L)  and  Gn^fL)  of  the  composite  system  in  question. 


5.  STRENGTH  OF  THICK  PLIES 

In  analogy  with  (3),  we  define  thick  plies  by  the  relationship 

2d  <4  2a,  for  <5  ^4,  (10) 

and  adopt  the  configuration  of  the  crack  nucleus  in  Fig.  2.  There  is  no  interaction  between  the 
microcrack  and  the  adjacent  plies;  hence  ij,  =  1  in  (6),  and  the  crack  width  is  now  equal  to  <5  rather 
than  a.  For  a  nucleus  of  width  6  one  obtains: 

G(L)  =  i7td[A52ff22  +  A24d33  +  AgsCTj2].  ford^a.  (11) 

This  release  rate  is  appropriate  if  the  nucleus  propagates  in  direction,  as  a  type  L  crack.  Of 
course,  propagation  in  the  .x,  direction  is  also  possible.  Such  cracking  is  referred  to  as  type  T 
cracking.  The  corresponding  energy  release  rate  is 


G(T)  =  dWi/dd, 


(12) 


First  ply  failure 
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and  from  (4)  and  (11)  it  follows  that 

G(T)  =  i7t^[Af2di2  +  A24ff23+^66di2],  for  (5  «  a.  (13) 

or 

G(T)  =  2G(L).  (14) 

At  this  point  we  again  appeal  to  the  mixed  mode  failure  criterion  (8).  However,  we  allow  for 
possible  differences  between  toughness  values  in  the  two  modes  by  writing: 


G,.(L)  =  7,G,JT) 

(15) 

^llc(^)  “  rll^IIcfT)- 

Results  reported  in  [9]  and  [10]  indicate  that  y,  ranges  from  0.6  to  0.8  for  certain  graphite-epoxy 
plies.  Data  of  this  kind  are  apparently  very  rare,  therefore  no  general  conclusions  about  y,  and  y,, 
magnitudes  can  be  drawn  at  this  time. 

To  determine  if  a  given  crack  nucleus  of  width  4  will  propagate  as  a  type  T  or  type  L  crack,  one 
can  evaluate  the  critical  values  6^(7)  and  6^(L).  For  example,  for  modes  I  and  II  taken  separately,  one 
can  find  from  ( 14)  and  (S)  that 


^„(L)  =  2y,<5,.(T) 

(16) 

<5„.(L)  =  2y„<J„JT). 

In  a  mixed  mode  case  one  needs  to  invoke  (8).  Since  available  data  indicate  that  2yi  >  1 ,  one  finds  that 
<),5(T)  <  <5,c(L).  If  one  assumes  that  <5j(T)  <  (5j(L)  unoer  mixed  loading  as  well,  then  it  follows  that 
first  ply  failure  of  thick  plies  occurs  as  a  result  of  transverse  (type  T)  cracking,  followed  by  longitudinal 
{type  L)  cracking. 

The  crack  starts  to  propagate  as  the  type  T  crack  in  the  x,  direction.  After  the  crack  has 
propagated  under  constant  stress  to  width  6  =  4(L),  it  also  may  start  to  propagate  in  the  x, 
direction,  as  the  type  L  crack.  This  bidirectional  propagation  continues  until  2S  has  reached  the  ply 
thickness  2<i.  In  any  case,  the  type  L  cracking  continues  in  the  Xj  direction  at  any  6  >  6,j(L)  across  the 
entire  loaded  area  of  the  ply  where  (8)  is  satisfied. 

The  results  derived  so  far  suggest  that  if  many  thick  plies  of  different  thickness  are  tested  in  such  a 
way  that  the  loading  history  s(r)  is  kept  constant  in  each  test,  then  6^  in  (2)  and  ( 16)  are  also  consunt, 
and  the  measured  strength  must  be  constant,  regardless  of  actual  ply  thickness. 

For  ffjj  =  0,  022  O'  ihe  strength  is  given  by 

(^22)c'  =  [2G,.(T)/(irAS2  AclT))]  ( 17) 

For  022  =  0,  ^23  5^  O'  the  strength  is  given  by 

(ff23)c,  =  [2G„,(T)/(2rA24  ^„c(T))]‘'^  ( 18) 

Constant  strength  is  seen  in  experimental  results  discussed  in  the  sequel.  Similar  relations  can  be 
found  from  (8)  for  combined  loading. 

Equations  ( 17)  and  ( 18)  can  be  used  to  predict  strength  only  if  the  (3,(T)  and  G,.(T)  are  known. 
This  is  usually  not  the  case  since  6j(T)  cannot  be  observed.  Therefore,  actual  values  of  thick  ply 
strength  need  to  be  measured  directly  in  experiments. 

One  possible  application  of  ( 17)  and  ( 18)  is  in  experimental  work  aimed  at  evaluation  of  growth 
rates  of  debonded  regions  under  different  loading  conditions.  Such  results  would  provide  an 
experimental  foundation  for  (2). 

In  conclusion  of  the  theoretical  part  we  note  that  one  can  also  examine  failure  of  plies  of 
intermediate  thickness,  where  the  crack  nucleus  interacts  with  adjacent  plies.  This  problem  is 
discussed  in  [1].  As  in  the  case  of  thick  plies,  implementation  of  the  analysis  requires  knowledge  of 
material  properties  which  is  not  available  at  this  time. 

It  is  clear  that  a  transition  exists  between  the  constant  strength  of  thick  plies  ( 17),  (18),  and  the 
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Fig.  4.  Comparison  of  theoretical  results  with  experimental  data. 


variable  strength  of  thin  plies  (9).  This  transition  takes  place  at  a  certain  value  of  ply  thickness  2a  that 
can  be  determined  from  strength  measurements  on  plies  of  different  thickness;  as  shown  in  the  sequel. 

6.  COMPARISON  WITH  EXPERIMENTS 

To  illustrate  the  theoretical  derivation,  we  interpret  results  presented  by  Wang[10],  Batley  et 
a/.[ll],  and  by  Crossman  and  Wang[12].  Figure  4  shows  the  first  set  of  results.  Several  (0/90,), 
coupons  made  of  an  E-glass  epoxy  laminated  plate  were  loaded  in  tension,  and  the  laminate  stress  or 
strain  was  recorded  at  first  failure  of  the  90°  ply.  To  interpret  the  data,  laminated  plate  theory  was 
used  to  calculate  initial  thermal  stresses  in  the  90°  ply  after  cooling  from  the  curing  temperature  to 
room  temperature,  d(f  =  ~  125°C.  Then,  the  transverse  normal  stress  <722  caused  in  the  90°  ply  by 
loading  of  the  laminate  was  found,  and  superimposed  with  thermal  stress.  In  this  way,  the  original 
data  points  were  converted  to  ply  stress  ^22.  and  then  plotted  in  Fig.  4.  Thermoelastic  properties  of 
the  ply  were  taken  from  Table  1. 

Next,  ply  toughness  G,c(L)  =  250  J/m^  was  found  by  inversion  of  (9a)  such  that  the  resulting 
(^22)0,  curve  fitted  the  experimental  points  at  low  ply  thicknesses. 

Then,  experimental  points  for  thick  plies  were  used  to  find  the  average  strength  of  thick  plies  as 
(<722)cr  =  0.089GPa.  The  transition  between  thin  and  thick  plies  was  found  at  2a  =  0.5mm  in  this  case. 

No  definite  evaluation  of  5,t(T)  can  be  made,  because  G,j(T)  was  not  measured.  However, 
expected  values  at  specified  0,^  are  shown  in  Table  2.  For  example,  if  G,,(T)  is  taken  as  equal  to 
G|c(L)  =  250  J/m^,  then  2^,5(T)  =  0.288  mm.  In  any  event,  ^tfT)  <  a.  At  the  transition  point  between 
thin  and  thick  plies  2a  =  0.5  mm,  hence  2d,,(T)  <  0.5  mm.  At  the  measured  strength  of  thick  plies 
(0.089  GPa),  this  suggests  that  G|c(T)  <  433  J/m^.  In  reality,  both  2^|,(T)  and  G|c(T)  should  be  much 
smaller,  so  that  type  T  cracking  may  take  place  even  as  ply  thickness  approaches  the  transition  value 
of  0.5  mm. 

Figure  5  shows  similar  results  for  a  T300/934  graphite/epoxy  system,  which  appeared  in  [10, 12]. 
Data  for  three  different  layups,  ( ±25/90,),,  (252/-252/9O2),,  and  (0/90„/0),  were  superimposed  to 


Table  I.  Material  properties  of  selected  plies 

E-Glass-epoxy  T300/934 
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Table  2.  Critical  'widtHs  <)|c  of  type  T  crack  nuclei  in  thick  plies,  calculated  from 

cq.(17) 


Cic(T)  (J/m**) 

E-Glass-epoxy 

2(515(11  (mm) 

‘  T3<X)/934 

26ic(T)  (mm 

100 

0.115 

0.I8S 

200 

0.231 

0.372 

400 

0.461 

0.746 

600 

0.692 

1.117 

800 

0.922 

1.489 

Thick  ply  strength 

(GPa) 

0.089 

0.064 

INNER  90*  PLY  THICKNESS  2a  (mm I 

Fig.  5.  Comparison  of  theoretical  results  with  experimental  data. 

find  G,,(L)  in  the  90“"  plies.  The  results  are  apparently  not  significantly  alTected  by  the  differences  in 
layup  of  outside  layers.  Note  that  the  0,^  values  shown  in  Figs  4  and  5  are  actually  values  of 
with  taken  as  equal  to  unity  because  of  interface  damage.  If  was  taken  as  equal  to 
calculated  values,  e.g.  =  0.8  in  E-glass-epoxy,  then  one  would  obtain  =  200  J/m^  in  Fig.  5, 
which  compares  better  with  the  Gi^  =  120±30  J/m^  indicated  for  this  material  in  Table  1  of  [11]. 

7.  CONCLUSIONS 

(a)  First  ply  failure  stress  of  plies  of  different  thicknesses  can  be  predicted  theoretically,  providing 
that  ply  toughness  values  and  initial  flaw  sizes  are  known.  In  cases  of  thin  plies  the  ply  thickness 
itself  determines  the  initial  flaw  size  and  ply  strength  can  be  directly  related  to  toughness  via  eqs 
(9).  In  case  of  thick  or  moderately  thick  plies  the  analysis  cannot  be  easily  implemented  for  lack  of 
experimental  information  about  relevant  material  properties. 

(b)  If  a  ply  of  any  thickness  contains  a  through-the-thickness  flaw — such  flaws  can  be  caused  by 
impact,  penetrations,  or  similar  types  of  damage — then  it  must  be  regarded  as  thin  ply  in  strength 
estimates.  Therefore,  it  is  conservative  to  regard  all  plies  as  thin  plies,  and  neglect  the  additional 
strength  which  may  be  found  in  undamaged  thick  plies. 

Acknowledgement — This  work  was  supported  by  a  grant  from  the  Air  Force  OfTice  of  Scientiric  Research.  Dr  M.  Hejaziand  Mr 
C.  J.  Wung  provided  assistance  in  interpretation  of  experimental  data. 

REFERENCES 

[I]  G.  J.  Dvorak  and  N.  Laws.  Analysis  of  progressive  matrix  cracking  in  composite  laminates.  11.  First  ply  failure.  J.  compos. 
Mater,  in  press. 


770 


G.  J.  DVORAK  and  N.  LAWS 


'21  A.  N.  Stroh.  Diilocations  and  cracka  in  aniaotropic  daaticity.  Phil.  Mag.  3,  62S  (1938). 

'3'  N.  Laws.  A  note  on  interaction  energiea  aaaociated  with  cracka  in  aniaotropk  media.  Phil.  Mag.  36,  367  (1977). 

'4'  G.  i.  Dvorak,  N.  Lawa  and  M.  Hejazi,  Analyaia  of  progreaaive  matrix  cracking  in  compoaite  laminaiea.  I.  Thermoelaatk 
propertiea  of  a  unidirectional  compoaite  with  cracka.  J.  compos.  Main.  IS  ( 1983). 

[3]  N.  Lawt.  G.  J.  Dvorak  and  M.  Hejazi,  Stiffheaa  changea  in  unidirectional  compoaitea  cauaed  by  crack  syatema.  Mech. 
Mater.  2,  123  (1983). 

[6]  K.  W.  Garrett  and  J.  E.  Bailey,  Multiple  tranaverae  fracture  in  90°croaa>ply  laminatea  of  a  glaaa  Obre-reinforced  polyester. 
J  Mater.  Sci.  II,  137  (1977). 

[7]  K.  L.  Reifsnider.  W.  W.  Stinchcomb,  E.  G.  Henneke  aitd  J.  C.  Duke,  Fatigue  damage-strength  relationships  m 
composite  laminates.  AFWAL-TR-83-3084,  Vol.  1  (September  1983). 

[8]  H.  T.  Hahn  and  T.  Johnneason,  Fracture  of  unidirectional  composites,  theory  and  applications,  in  Mechanics  of 
Composite  Materials — 1983  (Edited  by  G.3.  Dvorak),  AMD,  Vol.  38.  p.  133.  New  York  (1983). 

[91  W.  D.  Bascom.  Private  communication  ( 1984). 

flO^  A.  S.  D.  Wang,  Fracture  mechanics  of  sublaminate  cracks  in  composite  materials.  Compos.  Technol.  Rev.  6, 43  ( 1984). 
1 1]  J.  E.  Bailey,  P.  T.  Curtis  and  A.  Parvizi,  On  the  transverse  cracking  and  longitudinal  splitting  behavior  of  glass  and 
carbon  fibre  reinforced  epoxy  cross  ply  laminates  and  the  effect  of  Poisson  and  thermally  generated  strain.  Proc.  R.  Soc. 
A366,  399  ( 1979). 

[12]  F.  W.  Crossman  and  A.  S.  D.  Wang,  The  dependence  of  transverse  cracking  and  ddamination  on  ply  thickness  in 
graphite/epoxy  laminates.  ASTM  STP  T7S,  118  (1982). 


/Jir.  1.  SolUb  Sintcims  Vol.  23.  No.  9.  pp.  1269-1283,  1987 
Ptiaud  in  Great  Bnuin. 


0020-7683/87  13.00+  .00 
<0  1917  Ferpmoo  ioiiraalf  Ltd. 


THE  EFFECT  OF  FIBER  BREAKS  AND  ALIGNED 
PENNY-SHAPED  CRACKS  ON  THE  STIFFNESS 
AND  ENERGY  RELEASE  RATES  IN 
UNIDIRECTIONAL  COMPOSITES 

N.  Laws 

Department  of  Mechanical  Engineering,  University  of  Pittsburgh,  Pittsburgh,  PA  1 5261, 

U.S.A. 


G,  J.  Dvorak 

Department  of  Civil  Engineering,  Rensselaer  Polytechnic  Institute,  Troy,  NY  12181,  U.S.A. 

(Received  4  June  1986 ;  in  revised  form  1 1  November  1986) 

Abstract — The  loss  of  stiffness  due  to  penny-shaped  cracks  associated  with  fiber  breaks  in  a  uni¬ 
directional  composite  is  the  main  theme  of  the  paper.  Explicit  results  are  given  for  the  non¬ 
trivial  Hashin-Shtrikman  bound  together  with  the  estimates  obuined  from  the  self-consistent  and 
differential  schemes.  In  addition  the  paper  contains  some  results  on  energy  release  rates  for  two 
different  crack  growth  mechanisms.  It  is  shown  that,  in  theory,  the  differential  scheme  enjoys  a 
distinguished  position. 


1.  INTRODUCTION 

This  paper  is  mainly  concerned  with  the  changes  of  stiffness  and  strength  of  a  unidirectional 
fiber-reinforced  solid  due  to  fiber  breaks  accompanied  by  penny-shaped  cracks  at  the  ends 
of  the  broken  fibers. 

The  literature  on  the  effect  of  crack  distributions  on  the  response  of  solids  falls  into 
two  categories.  First  one  can  consider  periodic  distributions  of  cracks,  see  e.g.  Delameter 
et  u/.[l].  Second,  and  this  is  by  far  the  major  part  of  the  literature,  one  can  consider 
random  distributions  of  cracks.  Within  this  second  category,  the  pioneering  paper  is  due  to 
Budiansky  and  0'Connell[2],  although  a  less  general  approach  was  given  independently  by 
Salganik[3].  The  work  of  Budiansky  and  0’Connell(2]  was  especially  directed  at  predicting 
the  loss  of  stiffness  of  an  isotropic  solid  due  to  a  volume  distribution  of  randomly  oriented 
elliptical  cracks — leading  to  isotropy  of  the  cracked  solid.  Additional  work  on  the  effect  of 
slit  cracks  on  the  stiffness  of  anisotropic  solids  has  been  given  by  Gottesman  et  a/.[4),  and 
by  Laws  and  Dvorak[5, 6]. 

It  is  particularly  relevant  to  the  analysis  of  this  paper  to  call  attention  to  the  work  of 
Hocnig[7]  who  was  the  first  to  consider  anisotropic  distributions  of  cracks.  In  addition,  it 
is  important  to  recognize  the  contributions  of  Mura  and  Taya[8]  and  Taya[9]  which  were 
addressed  to  the  problem  of  determining  the  effect  of  fiber  breaks  on  the  response  of 
unidirectional  fiber-reinforced  materials. 

In  this  paper  we  extend  the  work  of  various  authors[7-9]  in  that  we  obtain  self- 
consistent  estimates  for  the  reduction  in  stiffness  of  unidirectional  composites  containing 
penny-shaped  cracks.  In  addition  we  obtain  the  only  non-trivial  Hashin-Shtrikman  bound 
on  the  moduli  of  the  cracked  solid.  We  also  derive  the  appropriate  differential  scheme  model 
for  the  loss  of  stiffness. 

We  show  how  both  the  self-consistent  and  differential  scheme  results  are  entirely 
consistent  with  the  Hashin-Shtrikman  bound.  Further  we  show  explicitly  that  the  results 
of  Mura  and  Taya[8, 9]  coincide  with  the  Hashin-Shtrikman  bound. 

As  a  further  illustration  of  the  results  presented  herein,  we  pay  brief  attention  to 
isotropic  solids  containing  distributions  of  aligned  penny-shaped  cracks  and  amongst  other 
things,  recover  some  results  first  given  by  Hoenig[7).  The  relevance  of  this  analysis  in  the 
study  of  micro-cracking  in  ceramics  is  discussed  by  Laws  and  Brockcnbrough(lO). 
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Finally  we  investigate  the  influence  of  cracks  on  energy  release  rates.  Here,  we  consider 
two  distinct  cases.  First,  we  consider  situations  in  which  all  cracks  extend  simultaneously. 
Second,  we  consider  only  one  crack  in  the  effective  cracked  solid.  It  turns  out  that  the 
differential  scheme  enjoys  a  distinguished  position  as  far  as  the  evaluation  of  energy  release 
rates  is  concerned.  Whether  or  not  this  is  of  physical  significance  is  an  open  question. 

2.  PHYSICAL  MOTIVATION 

Consider  a  unidirectional  graphite-epoxy  composite  subject  to  loading  parallel  to  the 
fibers.  It  is  known  that  such  composite  specimens  eventually  suffer  fiber  breakage.  The 
density  of  such  fiber  breaks  can  be  quite  large  before  final  fracture  occurs.  At  the  end  of 
each  broken  fiber,  one  usually  finds  penny-shaped  cracks  (whose  normals  are  in  the  fiber 
direction).  Under  continued  loading,  these  cracks  may  extend^see  the  discussion  of  Mura 
and  Taya[8, 9]  and  the  references  contained  therein. 

In  addition  there  is  ample  evidence  that  in  cross-ply  composite  laminates,  one  often 
sees  fiber  breaks  in  the  0°  plies  at  the  end  of  transverse  cracks  in  the  90°  plies.  Furthermore, 
some  recent  work  by  Laws  and  Brockenbrough[10]  indicates  the  relevance  of  considering 
aligned  penny-shaped  micro-cracks  in  ceramics. 

Another  application  arises  in  fibrous  composites  which  can  develop  aligned  matrix 
cracks  on  planes  perpendicular  to  the  fibers.  Such  cracks  are  confined  to  the  matrix  and  do 
not  extend  through  the  fibers.  They  have  been  observed  in  0°  plies  of  cyclically  loaded 
boron-aluminum  laminates  by  Dvorak  and  Johnson[l  1]. 

It  is,  therefore,  desirable  that  the  effects  of  aligned  penny-shaped  cracks  on  the  stiffness 
and  strength  of  both  composites  and  single  phase  materials  be  more  fully  understood. 


3.  ANALYSIS 

In  this  section  we  present  the  mechanics  of  the  effect  of  distributions  of  aligned  penny¬ 
shaped  cracks  on  the  stiffness  and  strength  of  solids.  For  simplicity  we  shall  discuss  cracked 
composites  by  regarding  the  uncracked  fibrous  composite  as  an  effective  homogeneous 
material.  This  is,  of  course,  commonplace  in  the  theory  of  composites,  but  some  further 
explanation  is  required  here  since  the  cracked  composite  will  contain  both  long  and  short 
(broken)  fibers. 

Now  Laws  and  McLaughlin[  1 2],  followed  by  Chou  et  al.  [13],  have  quantified  the  effects 
of  fiber  aspect  ratio  on  the  overall  moduli  of  aligned  short  fiber-reinforced  composites.  The 
general  conclusion  is  that  for  aspect  ratios  greater  than  100,  say,  the  stiffness  of  the 
composite  is  insensitive  to  fiber  length.  Since  the  applications  discussed  in  Section  2  all 
indicate  broken  fibers  whose  aspect  ratios  are  much  larger  than  100,  there  is  no  loss  of 
generality  in  considering  only  long  fibers.  If  we  interpret  correctly.  Taya[9]  arrives  at  the 
same  conclusion,  by  somewhat  different  methods. 

In  addition  for  transverse  cracks  in  composite  laminates  the  use  of  a  model  with  cracks 
in  an  otherwise  homogeneous  solid  is  well  established. 

The  notation  and  basic  ideas  presented  here  are  taken  directly  from  the  work  of  Laws 
and  Dvorak[5,6].  Fourth-order  tensors  are  denoted  by  upper  case  letters,  e.g.  L,  A  and 
symmetric  second-order  tensors  are  denoted  by  bold-face  letters,  e.g.  s,  a.  The  unit  fourth- 
order  tensor  is  denoted  by  /  and  the  inverse  of  a  non-singular  fourth-order  tensor  A  is 
denoted  by  A''. 

Consider  a  linear  elastic  solid  whose  stress  tensor,  a.  and  strain  tensor,  «,  are  related 
through 


a  =  Ls,  e  =  Ma,  LM  =  ML  =  /. 


(1) 


We  use  a  standard  6x6  matrix  notation  for  the  stiffness  tensor  L  and  compliance  tensor 
M.  Since  we  are  here  concerned  with  materials  which  are  at  worst  transversely  isotropic 
with  respect  to  the  coordinate  axis  O.r,,  it  follows  that  eqns  (1)  may  be  written  in  the  form 
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(2) 


(3) 


(4) 


We  use  the  suffix  0  to  refer  to  the  initial  uncracked  solid.  Thus,  fur  example  Lo  is  the 
stiffness  of  the  uncracked  solid,  whereas  L  is  the  stiffness  of  the  cracked  solid.  The  volume 
concentration  of  cavities  is  denoted  by  c. 

It  is  convenient,  but  not  essential,  to  develop  the  theory  by  considering  a  family  of 
aligned  spheroidal  cavities  in  a  transversely  isotropic  matrix.  A  typical  spheroid  is  taken  to 
be 


X 

a 


2 

2 


1. 


(5) 


In  addition  we  suppose  all  cavities  to  be  of  equal  size.  This  assumption  is  not  essential  as 
far  as  theory  is  concerned,  but  may  well  be  essential  in  any  application  to  penny-shaped 
cracks  at  the  ends  of  broken  fibers  in  a  unidirectional  composite.  Ultimately  we  will  obtain 
the  required  results  for  aligned  penny-shaped  cracks  by  allowing  the  aspect  ratio 


£  =  bia 


(6) 


to  approach  zero. 

Many  methods  have  been  proposed  to  predict  the  effect  of  reinforcement  on  the 
effective  moduli  of  composites.  It  is  not  our  purpose  here  lo  give  a  critical  survey  of  the 
various  methods.  Rather,  we  give  the  required  formula  for  the  Hashin-Shtrikman  bounds, 
the  self-consistent  method  and  the  differential  scheme. 

In  order  to  get  the  Hashin-Shtrikman  bounds  we  refer  to  the  work  by  Willisfl4].  We 
assume  that  the  cavities  are  randomly  located  so  that  the  cracked  solid  exhibits  overall 
transverse  isotropy.  Since  this  envisaged  distribution  of  oblate  spheroidal  cavities  will 
conform  to  the  statistics  assumed  by  Willis[141,  it  follows  that  the  overall  compliance  of 
the  cracked  composite  must  satisfy 


M  ^  Mo  +■  - Qo'- 

I  — c 


(7) 


Here  Qo  is  a  tensor  which  depends  on  the  uncracked  compliance  Mo  and  the  aspect  ratio 
£.  The  components  of  Qocan  befound  explicitly,  sec  Laws  and  McLaughlin[12]  or  Laws[I5]. 
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The  right-hand  side  of  inequality  (7)  is  the  Hashin-Shtrikman  lower  bound  for  M ,  denoted 
here  by  Af  ‘ 


M  =  Mo  -f-  - - Qo  '• 

I  — c 


Since  the  damaged  material  contains  cavities  the  Hashin-Shtrikman  upper  bound  for  M  is 
here  infinite. 

Next  we  recall  that  the  self-consistent  estimate  is  given  by  Laws  et  a/.  [5]  as 


M  =  Mo+cQ-'. 


Finally,  it  is  easy  to  read  off  the  differential  scheme  estimate  from  the  work  of  McLaugh- 
lin{16] 


dc  1-c'^ 


M  =  Mq  when  c  =  0. 


In  terms  of  the  stiffness  tensor  L  it  is  easy  to  show  that  the  upper  bound  L*  is  given  by 


Also  the  self-consistent  result  is 


L  =  Z.o(/+cQ-'Lo]‘ 


whereas  the  differential  scheme  gives 


dL  1 


dc  1  — c 


LQ-'L 


L  =  Lo  when  c  =  0. 


We  now  obtain  the  required  results  for  penny-shaped  cracks  by  proceeding  to  the  limit 
as  e  -*  0.  Let  rj  be  the  number  of  cracks  per  unit  volume  then 


c  =  -  ntja  e 


n 

=  -ac 


where  the  crack  density  parameter  a  is  defined  by 


a  =  8>fa\ 


Note  that  the  crack  density  parameter  is  noi  the  same  as  the  crack  density  parameter 
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introduced  by  Budiansky  and  0’Connell[2].  The  choice  of  eqn  (16)  is  motivated  by  the  fact 
that  a  =  1  corresponds  to  an  average  of  one  crack  of  diameter  2a  in  each  cube  of  side  2a. 

In  applications  of  the  theory  to  composites  reinforced  by  aligned  continuous  fibers, 
the  crack  can,  in  the  first  instance,  be  visualized  as  a  fiber  break  which  has  extended  into 
the  surrounding  matrix  and  has  been  arrested  by  adjacent  unbroken  fibers.  This  crack  is 
then  regarded  as  a  crack  in  an  effective  composite  medium.  For  practical  purposes  it  is 
desirable  to  relate  the  radius  of  this  crack  to  the  fiber  radius,  rr,  and  the  fiber  volume 
fraction,  ty.  Estimates  can  be  obtained  in  various  ways.  For  example  if  the  microstructure 
is  such  that  the  fibers  are  located  in  a  close  packed  hexagonal  array  and  that  failure  of  one 
fiber  creates  a  crack  which  extends  through  the  matrix  until  it  reaches  the  surfaces  of 
neighboring  fibers,  then  the  resulting  penny-shaped  crack  has  radius 


In  the  same  spirit  one  obtains  the  estimate 

for  a  square  array.  For  typical  volume  fractions  occurring  in  practice  the  hexagonal  array 
estimate  exceeds  the  square  array  estimate  by  about  13%. 

As  noted  by  Eshelby[17]  in  the  isotropic  case,  and  by  Hoenig[7]  and  Laws[lS]  for 
orthotropic  materials,  the  limit  as  e  approaches  zero  needs  to  be  handled  with  care.  The 
essential  point  here  is  that  whereas  Q  becomes  singular,  the  product  bQ~  '  remains  finite. 
Thus  let 


limrg"'  *  A 

with  an  analogous  definition  for  Aq.  The  components  of  A  are  given  in  the  Appendix. 
In  the  limit  of  aligned  penny-shaped  cracks  it  now  follows  from  eqn  (8)  that 


A/  =  Mo  -h  ^aAo. 
o 


(19) 


Further,  the  self-consistent  estimate  is  obtained  from 


A/=  A/„  -f--aA 
0 


(20) 


whereas  the  differential  scheme  yields 


dM  n 
da  6 


(21) 


The  dual  formulae  for  the  stiffness  tensor  are 


L  —  £.o|^/  + 


(22) 


the  self-consistent  result  being 
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=  Lo 


/  +  7  ol\Lq 
o 


and  the  differential  scheme  giving 


:r  =  ~  Z 

da  6 


Since  the  only  non-zero  components  of  A  (shown  in  the  Appendix)  are  Aj,  and  A 44 
(=  A55)  it  follows  from  eqns  (20)  and  (21)  that  only  and  iV/44  (=  M^)  are  predicted 
to  change.  This  is,  of  course,  only  to  be  expected. 

It  is  instructive  to  consider  the  rather  special  case  when  the  uncracked  body  is  isotropic. 
Thus  let  £0,  vq  be  respectively  the  initial  Young’s  modulus  and  Poisson’s  ratio.  Then,  from 
the  Appendix,  we  see  that 


= 


4(1 -V5) 


8(1 -v^) 

n£o(2-Vo)' 


With  the  help  of  a  standard  notation 


A/)3  —  — ,  A/44  — 


it  now  follows  from  eqn  (19)  that 


£0  i-l-ia(l-v5) 


Co  l+^a(l-Vo)/(2-Vo)' 

We  emphasize  that  eqn  (25)  is  precisely  the  formula  given  by  Taya[9].  Hence  we  see 
that  the  Mori-Tanaka[19]  back-stress  analysis,  which  is  used  by  Taya[9],  yields  the  non¬ 
trivial  Hashi'-Shtrikman  bound.  It  is  noteworthy  that  it  is  possible  to  show  that  the  assump¬ 
tions  of  the  Mori-Tanaka[19]  back-stress  analysis  lead  to  formulae  which  are  coincident 
with  the  Hashin-Shtrikman  bounds  in  other  situations — but  we  do  not  include  details  here. 

We  remark  that  there  is  no  difficulty  in  obtaining  eqns  (19)-(21)  by  direct  methods 
instead  of  using  an  argument  based  on  cavities.  In  fact  Gottesman  et  a/.[4]  have  given  such 
an  argument  to  obtain  the  bounds  and  self-consistent  results  for  a  solid  with  slit  cracks. 
We  note  that  the  fully  general  analysis  is  given  by  Laws  and  Brockenbrough[10]. 

Numerical  results  can  be  found  from  eqns  (19)  to  (21)  or  from  eqns  (22)  to  (24).  In 
practice,  we  have  found  it  easier  to  evaluate  the  compliances  from  eqns  (19)  to  (21)  and 
then  to  determine  the  stiffness,  when  required,  by  matrix  inversion. 

For  the  most  part  we  present  results  for  the  loss  in  stiffness  of  unidirectional  fiber- 
reinforced  materials  containing  a  distribution  of  penny-shaped  cracks  with  common  orien¬ 
tation  perpendicular  to  the  fiber  direction.  In  other  words  we  attempt  to  model  the  loss  of 
stiffness  of  unidirectional  fiber-reinforced  materials  due  to  fiber  breaks.  Data  for  the 
uncracked  composite  is  taken  directly  from  Table  2  of  the  paper  by  Dvorak  et  a/. [6]  which 
is  for  a  T300/5208  graphite-epoxy  system  with  volume  fractions  Cf  =  0.2, 0.4. 0.6.  As  noted 
earlier  the  only  compliances  which  change  are  Mu.  A/44  ( =  A/ss)-  In  Figs  1  and  2  we  give 
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the  self-consistent  and  differential  scheme  predictions  for  these  two  compliances.  There  is 
no  need  to  display  the  Hashiu-Shtriknian  bound  since  it  coincides  with  the  common  tangent 
at  a  =  0.  We  observe  that  for  each  value  of  C(  the  self-consistent  result  is  always  greater 
than  the  differential  scheme  result.  For  small  a  (<0.l  say)  the  two  results  for  Afjs  are 
coincident  and  equal  to  the  Hashin-Shtrikman  bound.  On  the  other  hand,  we  see  from  Fig. 
2  that  there  is  virtually  no  difference  between  the  self-consistent,  the  differential  scheme  and 
the  Hashin-Shtrikman  results  for  jV/44  for  a  <  0.5.  Thus  for  practical  purposes  it  suffices 
to  take 
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Fig.  4.  Longitudinal  shear  modulus  for  various  T300/5208  systems:  (a)  bound  - ; 

(b)  s.c.m. - ;  (c)  d.s. . 

^  =  (l  -f  ^a£2A?,^ 

|^=(l+^aGgAj.)  . 

We  note  that  the  self-consistent  estimate  for  both  El  and  Gl  is  always  lower  than  the 
differential  scheme  result.  Clearly  both  results  are  consistent  with  the  Hashin-Shtrikman 
bound. 

As  for  the  reduction  in  vlt  we  note  that  since  M,,  =  and  M,,  =  M,,,  it  follows 

that 

v2t  El  • 

Hence  the  fractional  reduction  of  vlt  is  equal  to  the  fractional  reduction  of  El-  Also  the 
actual  reduction  in  v^t  can  be  found  from  the  reduction  of  El  merely  by  a  change  in  scale. 
For  clarity  Figs  5  and  6  display  only  the  fractional  reduction  for  El  and  Gl  when  C(  =  0.6. 
We  remark  that  the  appropriate  curves  for  Cf  =  0.2  and  0.4  are  almost  indistinguishable 
from  the  curves  indicated  in  Figs  5  and  6.  By  way  of  comparison  Figs  7  and  8  show  the 
fractional  reduction  of  El  and  Gl  for  an  isotropic  material  with  Poisson's  ratio  vq  =  0.3. 
Two  remarks  are  in  order.  First  the  fractional  reduction  curves  indicated  in  Figs  6  and  7  are 
insensitive  to  the  choice  of  Poisson’s  ratio.  0.2  ^  v<,  $  0.4.  Second  the  fractional  reductions 
indicated  in  Figs  5  and  6  compared  with  those  in  Figs  7  and  8  show  that  there  is  no 
possibility  of  constructing  "master  curves”  for  anisotropic  materials. 

4.  ENERGY  RELEASE  RATES 

In  this  section  we  calculate  energy  release  rates  for  two  different  mechanisms  of  crack 
growth.  First,  we  consider  a  solid  which  contains  a  family  of  cracks  each  of  radius  a  and 
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Fig.  8.  Normalized  longitudinal  shear  modulus  for  isotropic  materials  with  Poisson's  ratio  =  0.3  -. 
(a)  bound - ;  (b)  s.c.m. - ;  (c)  d.s. . 


allow  each  crack  to  extend  by  the  same  amount.  Second,  we  consider  a  single  penny-shaped 
crack  in  an  otherwise  homogeneous  (cracked)  solid.  We  show,  inter  alia,  that  the  differential 
scheme  enjoys  a  peculiar  status  in  that  it  predicts  equality  of  the  two  energy  release  rates. 

Thus  we  first  consider  a  solid  containing  a  family  of  aligned  penny-shaped  cracks  each 
of  radius  a.  The  corresponding  crack  density  is  a.  The  solid  is  subject  to  macroscopically 
uniform  loading  with  applied  stress  d,  which  is  also  the  average  stress.  All  cracks  arc 
assumed  to  be  open.  The  total  energy  of  the  cracked  solid  is 


TudS 


(27) 


where  T,  u  are  respectively  the  tractions  and  displacements  on  the  boundary  S.  It  is  easy 
to  show  that  eqn  (27)  may  be  rewritten  in  the  form 


£(a)  =  -  j  Kff  •  A/(a)ff  when  x  #  0, 

£o  =  —^Vd'  Mod  when  ot  =  0. 


(28) 


Thus  the  total  energy  released  by  the  introduction  of  cracks  is  (Eg  -  E(oi)).  Accordingly  the 
energy  released  by  each  crack,  is  given  by 


£o-£(a) 

r,y 


=  ^d’[M{%)-Mo]d. 


(29) 


Suppose  now  that  each  crack  extends  from  radius  a  to  radius  {a  +  5a)  while  the  total 
number  of  cracks  remains  constant.  This  corresponds  to  an  increase  in  crack  density  from 
a  to  (a  +  <5a).  The  crack  extension  force,  or  energy  release  rate  per  unit  length,  C*,  of  each 
crack  is  given  by 


G  -  -L  —  — 

^  2na  d'x  da 


(30) 


Here  the  subscript  A  is  used  to  signify  that  Ga  is  the  energy  release  rate  of  each  crack  when 
all  cracks  extend  simultaneously.  From  eqn  (16) 
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da 

d<2 


=  24^a^ 


so  from  eqns  (29)  and  (30) 


Ga  = 


(31) 


The  lower  bound  A/  in  eqn  (19)  gives 


—  QG  *  ^0^ 

while  the  self-consistent  estimate  of  M  in  eqn  (20)  yields 


and  the  differential  scheme  result,  eqn  (21),  gives 


=  ao’  \a. 


Next  we  focus  on  a  single  penny-shaped  crack  in  the  cracked  solid.  We  now  regard  the 
cracked  solid  as  a  homogeneous  medium  with  effective  compliance  A/(a).  Only  one  crack 
is  now  present  in  this  effective  medium  and  we  wish  to  determine  the  energy  release  rate, 
Gs,  for  extension  of  this  single  crack. 

In  these  circumstances  we  can  use  the  results  of  Laws[l  5]  to  obtain  the  required  energy 
released  by  extension  of  this  single  crack  of  radius  a 

(a)  =  \a. 


Hence 


Gs 


1 

2ita  da 


ad  •  \d. 


(32) 


We  emphasize  that  different  models  give  rise  to  different  values  of  M.  Since  the  components 
of  A  are  given  in  terms  of  the  components  of  M  by  the  formulae  of  the  Appendix,  it  follows 
that  each  of  the  three  models  gives  rise  to  a  different  value  of  A  and  hence  of  Gs. 

Whether  we  consider  all  cracks  simultaneously  growing  in  self-similar  fashion  or  a 
single  growing  crack,  it  is  clear  that  both  G^  and  Gs  are  average  crack  extension  forces.  Of 
course  if  we  are  considering  crack  extension  in  a  single  phase  brittle  solid  there  is  no 
conceptual  problem.  However,  in  the  case  when  the  cracks  are  located  in  an  effective 
composite  matrix  some  further  comments  are  in  order.  Thus  consider  an  arbitrary  growing 
penny-shaped  crack  in  a  unidirectional  fiber-reinforced  composite — with  the  plane  of  the 
crack  perpendicular  to  the  fiber  direction.  In  the  first  place  this  crack  will  arise  from  a  fiber 
break  followed  by  extension  through  the  matrix  until  the  circumferential  crack  tip  reaches 
neighboring  fibers.  If  there  is  further  growth  then  part  of  the  crack  tip  must  lie  in  the  matrix 
material  whereas  the  remainder  must  lie  in  the  fibers.  Since  here  we  consider  cracks  in  an 
effective  solid  it  is  clear  that  the  energy  release  rates  Ga  and  Gs  do  not  apply  when  the  crack 
radius  is  smaller  than  the  value  given  in  eqn  (17)  or  eqn  (18).  However,  for  subsequent 
growth  an  average  crack  extension  force  is  precisely  the  quantity  which  is  demanded  by  the 
physics  of  the  problem. 

Taya(9]  has  reported  some  work  on  the  second  case  (Gs)  but  his  work  appears  to  have 
little  in  common  with  the  work  described  here. 
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Fig.  9.  Ratio  of  energy  release  rates  for  various  T300/S208  systems  according  to  s.c.m. 


From  eqns  (21),  (29)  and  (30)  we  obtain  the  surprising  result  that  the  differential 
scheme  predicts  that  Ga  =  Gs.  Thus  the  differential  scheme  can  provide  a  useful  border 
between  those  models  which  predict  Ga  to  be  greater  or  smaller  than  Gs-  However,  one 
ought  to  bear  in  mind  that  there  is  no  a  priori  reason  why  a  particular  model  should  not 
give  Ga  >  Gs  for  small  a  and  yet  Ga  <  Gs  for  large  a. 

Since  a  general  analysis  of  eqns  (29)  and  (30)  is  complicated  by  the  fact  that  both  refer 
to  mixed  mode  loading,  we  concentrate  on  mode  I  loading  in  which  the  only  non-zero 
applied  stress  is  From  eqns  (31)  and  (32)  we  have 

Ga  6  dMn/da 

Gs  ^  A33 

It  is  perhaps  useful  to  emphasize  here  that  both  Ga  and  Gs  depend  upon  crack  density.  This 
may  be  contrasted  with  the  usual  situation  in  fracture  mechanics  wherein  a  single  crack 
extends  in  a  material  of  fixed  properties.  Nevertheless  we  see  from  eqns  (31)  and  (32)  that 
if  Ga/Gs  >  1  extension  of  all  cracks  is  indicated  rather  than  extension  of  a  single  crack. 
When  Ga/Gs  <  1  the  opposite  conclusion  applies. 

Numerical  results  can  be  obtained  for  the  T300/5208  graphite-epoxy  systems  con¬ 
sidered  earlier.  In  particular  Fig.  9  shows  the  self-consistent  and  (trivially)  the  differential 
scheme  estimates  for  Ga/Gs.  Since  Ga/Gs  ^  1  both  models  suggest  extension  of  all  cracks. 
By  way  of  comparison  Fig.  10  shows  the  results  for  initially  isotropic  solids  in  which 
case  the  value  of  Ga/Gs  is  insensitive  to  the  choice  of  initial  Poisson's  ratio  in  the  range 
0.2  <  Vo  <  0.4. 

A  further  general  conclusion  on  energy  release  rates  for  a  single  crack  can  be  obtained 
from  eqn  (30)  and  the  results  of  Section  3.  For  simplicity  consider  the  mixed  mode  loading 
d)}  54  0,  d2)  #  0  so  that  eqn  (30)  reduces  to 

Gs  =  aA33((f,3)--(-aA44(ff23)^- 

Thus  A33  and  A44  may  be  interpreted  as  energy  release  rate  factors.  Indeed  Figs  1 1  and  12 
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contain  the  self-consistent  and  differential  scheme  results  for  A33  and  A44  for  the  T300/5208 
systems  considered  earlier.  It  is  clear  that  both  models  predict  that  at  fixed  crack  density, 
A33  and  A44  decrease  as  the  volume  fraction  of  fiber  increases.  Hence  Gs  decreases  as  the 
volume  fraction  of  fiber  increases  in  agreement  with  some  results  of  Taya[9].  Further  it  is 
evident  that,  for  fixed  Cf,  Gg  increases  with  a. 

Finally,  we  note  that  in  the  fibrous  composite  medium  the  crack  density  will  increase 
gradually  under  incremental  load,  as  a  result  of  fiber  breaks  at  randomly  distributed 
locations.  Each  new  crack  will  require  for  its  formation  at  least  the  same  amount  of  crack 
energy  as  the  first  crack,  i.e.  lV{a.)  ^  Wo-  Equation  (29)  indicates  that  this  will  be  the  case 
for  all  models  considered  herein.  Therefore,  progressive  cracking  is  governed  by  spatial 
variation  of  fiber  strength ;  sufficient  crack  energy  is  available  for  each  fiber  break  to  produce 
a  new  crack. 
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APPENDIX 

For  completeness  we  give  here  the  formulae  for  the  non-zero  components  of  the  A  tensor,  defined  in  eqn 
(18).  for  a  penny-shaped  crack  in  a  transversely  isotropic  material.  The  crack  lies  in  a  plane  normal  to  the  axis  of 
transverse  isotropy.  As  is  shown  by  Laws[15] 

,  2y.V!(Vi+y,)  Mii-iVfj 

M— 

_ 4(y,+y^)(M],-M]^)(2M„)'^ _ 

“  n{Mu(2M,,)'^  +  (y,+y2)(M„  +  Mn)iM„-M,,)''^] 

where  yj  and  yj  are  the  roots  of 

(M\,-M]r)x^-[MuM„  +  lM„(M,,-M,^)\x+MuM„-M\,=0. 

When  the  material  is  isotropic  with  Y oung's  modulus  E  and  Poisson's  ratio  v,  it  is  easy  to  see  that  y ,  =  y  j  =  1 . 
Hence 


)t£(2  — v) 


These  results  are  in  complete  agreement  with  those  of  Eshelby(I7J. 
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ABSTRACT 

In  their  past  work,  Dvorak  and  Johnson  had  shown  that  fatigue  daaage  in 
metal  composite  la^nates  takes  place  when  the  applied  loads  cause  cyclic 
plastic  straining  in  the  matrix.  The  present  work  describes  an  incremen¬ 
tal  analysis  of  the  plastic  deformation  and  damage  accunuletlon  processes 
in  elastic-plastic  composite  laminates.  It  is  baaed  on  a  strain  space 
foraulation  of  a  plasticity  theory  of  fibrous  composites  and  on  a  self- 
consistent  analysis  of  stiffness  changes  and  local  fields  In  cracked 
fibrous  layers  which  has  been  extended  here  to  elastic-plastic  laminates. 


1.  INTRODUCTION 

In  most  fiber  composite  systems,  the  strength  of  the  fiber  far  exceeds 
chat  of  Che  matrix,  hence  damage  caused  by  matrix  cracking  Is  comnonly  ob¬ 
served  in  fibrous  composites  and  laminates*  This  Is  certainly  an  unwel¬ 
come  disadvantage  of  composite  materials,  and  of  many  ocher  heterogeneous 
media  such  as  concrete  and  rocks,  which  La  not  encountered  to  a  similar 
extent  In  metals  and  Indeed  la  biosc  macroscoplcally  homogeneous  materials. 
The  difficulty  with  composites  Is  that  damage  may  start  at  relatively  low 
overall  loads  which  do  not  utilize  the  superior  strength  afforded  by  the 
reinforcement.  Typical  consequences  of  damage  are  reduction  of  stiffness, 
redistribution  of  Internal  stresses  which  may  impair  strength,  and  expo¬ 
sure  of  the  mlcrostruccure  to  possible  environmental  attack. 

In  some  systems  and  In  certain  applications,  damage  does  not  cause 
significant  problems.  For  exa^)le,  the  loss  of  stiffness  and  strength 
chat  may  be  caused  by  cracking  of  a  polymer  matrix  is  often  quite  small, 
because  the  matrix  itself  makes  only  a  limited  contribution.  On  the  other 
hand,  mechanical  properties  of  metal  and  ceramic  matrix  cooiposftes  are 
more  sensitive  to  the  effects  of  matrix  damage  because  the  matrix  may 
carry  a  major  part  of  the  applied  load  in  these  systems. 

In  the  last  fifteen  years,  damage  modeling  and  analysis  have  become  a 
major  activity  In  mechanics  of  composites.  Recent  surveys  by  Hashln 
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(1983)  and  Vang  (1984)  describe  some  of  the  earlier  work  that  pertains  to 
brittle,  mostly  polymer  systems.  Since  the  geometry  and  evolution  of 
damage  In  composites  are  often  rather  complex,  predictive  modeling  of  the 
process  and  of  Its  consequences  Is  very  difficult.  These  problems  seem  to 
be  magnified  In  metal  matrix  systems  where  damage  Is  seldom  caused  by 
monotonlc  or  short  variable  loads,  but  frequently  appears  under  sustained 
cyclic  plastic  straining  which  leads  to  fatigue  cracking  In  the  matrix. 
This  particular  problem  Is  the  subject  of  the  present  paper.  While  cer¬ 
tainly  not  denying  the  difficulties  Involved,  we  outline  a  relatively 
simple  micromechanical  approach  to  fatigue  damage  analysis  In  metal  matrix 
conq>o8lte  laminates.  The  next  section  presents  some  essential  results 
which  evaluate  the  effect  of  transverse  matrix  cracks  on  thermomechanical 
properties  and  stress  distribution  In  an  elastic  fibrous  ply.  Section  4 
reviews  some  basic  concepts  of  plasticity  of  fibrous  composites.  Section 
5  Is  a  discussion  of  certain  physical  aspects  of  damage  development  In 
mecal  matrix  laminates  which  Identify  damage  as  a  shakedotm  mechanism,  and 
pave  the  way  to  quantitative  modeling  of  damage  development.  The  details 
of  the  analysis  are  omitted  here,  but  selected  results  are  used  to  Illu¬ 
strate  the  key  results. 


2.  ELASTIC  PROPERTIES  OF  A  CRACKED  PLY 

The  first  step  in  the  analysis  Is  the  evaluation  of  overall  compliance  and 
stiffness  tensors  of  a  fibrous  ply  which  Is  embedded  in  a  laminated  plate 
and  contains  a  certain  density  of  transverse  cracks.  Fig.  1.  For  modeling 
purposes,  the  cracks  are  represented  by  aligned  slit  cracks  which  extend 
In  the  fiber  direction  chat  coincides  with  Che  X3  axis  in  Che  local  co¬ 
ordinate  system  of  each  ply.  The  remaining  parts  of  the  laminace  are  re¬ 
presented  here  by  layers  of  thickness  b  and  b*.  In  reality,  the  cracks 
are  not  uniformly  spaced  but  their  average  density  In  a  large  representa¬ 
tive  volume  V  of  Che  ply  can  be  described  by  a  crack  density  parameter 
8  ■  4a2n/V,  where  n  is  the  number  of  cracks  in  a  unit  area  of  the  trans¬ 
verse  xxX2  plane,  and  2a  Is  the  crack  width  equal  to  the  ply  thickness. 


Fig.  1 


Schematic  representation 
of  transverse  cracks  In 
a  fibrous  conposlte  ply. 
Fibers  are  aligned  In  the 

X  direction. 
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Suppose  Chat  the  stiffness  and  coo^llance  tensors  and  of  the 
ply  are  known  in  the  undamaged  state.  These  quantities" can  be" estimated 
In  terms  of  fiber  and  matrix  moduli  and  fiber  volume  fractions  by  approxi¬ 
mate  averaging  techniques,  or  bracketed  by  the  Hashln-Rosen  bounds.  Exact 
evaluation  of  the  corresponding  censors  for  a  damaged  ply  would  have  to 
account  not  only  for  the  interaction  between  the  cracks,  but  also  for  the 
Interaction  of  each  crack  with  the  adjacent  plies  of  arbitrary  orienta¬ 
tion.  To  circumvent  the  considerable  difficulty  posed  by  these  Inter¬ 
actions,  and  particularly  by  the  latter,  we  assume  that  the  properties  of 
Che  cracked  ply  can  be  approximated  by  those  of  a  homogeneous  medium  which 
contains  a  certain  density  of  aligned  silt  cracks  of  width  2a.  Then,  the 
overall  stiffness  L(B)  and  compliance  M(6)  at  certain  crack  density  8  can 
be  determined  by  the  following  self-consistent  approximation,  originally 
suggested  by  Budlansky  and  O’Connell  (1976). 

In  the  absence  of  cracks,  the  potential  energy  of  a  homogeneous  body 
of  volume  V,  under  uniform  overall  stress  o  can  be  written  In  the  form 

Eq  -  -4  V  o.M  0  .  (1) 

t  “  "O” 

The  energy  released  by  formation  of  a  single  crack  In  the  homogeneous  body 
can  be  written  In  terms  of  a  certain  crack  tensor  Aq  as: 

W  -  4  ir  a2  o.A  0  .  (2) 

I  "  “O" 

Now,  If  many  aligned  slit  cracks  are  introduced  In  such  a  way  that  the 
body  remains  statistically  homogeneous,  Chen  the  effect  of  Interaction 
between  the  cracks  on  Che  energy  released  by  a  single  crack  can  be  ex¬ 
pressed  In  terms  of 

W(6)  -  Y  w  a2  o.A(0)o,  (3) 

which  Is  equal  to  the  energy  released  by  a  single  crack  In  a  homogeneous 
medium  of  compliance  M(8),  or  stiffness  L(S).  Of  course,  the  number  of 
cracks  in  the  volume  V  Is  equal  to  n  •  0V/4a2.  Therefore,  the  potential 
energy  of  the  medium  with  open  slit  cracks  of  density  8  can  be  written  as 

-  Y  V  o.M(8)o  -  -  Y  V  o.MqO  -  y  h  w  a2  o.A(8)o  .  (4) 

This  must  hold  for  any  V  and  o,  hence  the  effective  compliance  of  the 
cracked  medium  is 

M(8)  -  Mo  8  A(8).  (5) 

Next,  use  the  Identity  ^  ■  I,  where  1  Is  the  identity  Censor,  to 

find  the  effective  stiffness 

"  to  ®  to  A(8)  L(8)  . 


(6) 


Evaluation  of  the  crack  censor  A(B)  In  a  transversely  Isotropic  medium 
with  aligned  slit  cracks  was  described  by  Laws,  Dvorak,  and  Hejazl  (1983), 
and  Dvorak,  Laws,  and  Hejazl  (198S)> 

The  overall  Chermoelastlc  constitutive  equations  of  the  cracked  ply 
can  be  written  in  the  form 


L  e  -  0  1 


M  0  +  0  m  , 

•m  ••  « 


(7) 


where  m  is  thermal  strain  vector  of  expansion  coefficients  and  I  - 
is  che'thermal  stress  vector.  One  can  establish  from  Che  above^analysls 
Chat  Q  ■  of  Che  undamaged  state,  and  chat  t  •  [I  ~  8L(8)  A(8)]tg. 
Also,  if  Che  overall  strain  of  the  cracked  medium  is  written^as  the  sum 


e  +  e 
“o  "c 


s  • 
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of  Che  average  strain  tg  in  the  uncracked  ligaments,  and  the  average 
strain  eg  accommodated  by  opening  of  all  the  cracks,  then  Dvorak  ec  al. 
(1985)  show  that 


-  A  €  -  0  a 


A  e  -  0  a 
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where 
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(10) 
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and  A  -  A(S),  L  -  L(8). 

Similar  results  can  be  developed  for  cracks  in  zero-degree  plies. 
Such  plies  typically  contal-n  matrix  cracks  on  planes  perpendicular  to  the 
fiber.  In  well-made  systems  such  cracks  bypass  Che  fibers,  and  the  cor¬ 
responding  crack  geometry  must  be  analyzed,  c.f.,  Wung  (1987). 


3.  RELAXATION  SURFACES 

Plastic  deforMCion  of  heterogeneous  media  is  a  very  complex  process.  The 
Initial  stages  are  dominated  by  local  yielding  at  inhomogeneities,  but 
such  contained  plastic  flow  does  not  significantly  affect  overall  stiff¬ 
ness.  Indeed  the  onset  of  overall  yielding  is  typically  associated  with 
extensive  plastic  straining  in  the  aggregate.  One  may  then  assume  that 
Che  overall  yield  condition  of  a  fibrous  ply  can  be  expressed  in  terms  of 
Che  matrix  yield  condition  and  the  average  matrix  stress.  Recent  experi¬ 
mental  observations  of  plastic  behavior  of  fibrous  composites  (Dvorak 
et  al.  1988)  suggest  that  in  certain  but  not  all  cases  Che  experimentally 
detected  overall  yield  surface  of  the  composite  medium  can  be  approximated 
in  this  way.  Several  composite  plasticity  models  have  been  developed 
along  these  lines;  in  the  sequel  it  will  be  convenient  to  use  Che  strain 
space  formulation  of  Wung  and  Dvorak  (1985). 


Let  Che  matrix  relaxation  condition  be  expressed  In  the  Mlses  form 
2  T  2 

g  (e  -  B  )  =  2G  (e  -  g  )  C(e  -  8  )  -  k  -  0  ,  (12) 

IB  IQ  ^  flS  m  IQ  D  OQ  Ill 

where  Gq  is  the  matrix  shear  modulus,  C  Is  a  constant  matrix,  Cq  Is  Che 
matrix  strain,  Is  Che  position  of  Che  center  of  the  relaxation  surface 
In  eg,  space,  and  k  Is  the  matrix  yield  stress  In  simple  shear.  If  the 
matrix  strain  concentration  factors  are  written  In  Che  form  used  In  (9), 
Chen  the  composite  relaxation  surface  In  Che  overall  strain  space  e  Is: 

2  T  T  2 

g(e  -  8)  -  2G  (e  -  6)  A  C  A  (e  -  8)  -  k  -  0  ,  (13) 

"  “  m  *  '  "me  "me 

where  the  concentration  factors  ^mc  '>>ust  be  derived  from  a  chosen  model  of 
Che  composite  mlcrostructure.  This,  and  the  development  of  hardening  and 
flow  laws  which  can  be  applied  to  fibrous  plies  In  a  metal  matrix  laminate, 
was  described  In  Che  above  reference. 

(Then  a  system  of  cracks  of  density  8  exists  In  a  ply  which  has  under** 
gone  some  plastic  strain  Cg  followed  by  unloading,  the  total  strain  In  the 
layer  with  fully  open  cracks  can  be  written  as 

e-Mo+e  ,  o«Le+o  (14) 

P  R 

where  M,  L  are  given  by  (5)  and  (6),  Cp  is  the  plastic  strain,  and  Is 
Che  relaxation  stress  which  Is  given  by  »  L  Cp.  Of  course,  Che  un* 
broken  ligaments  In  the  ply  are  now  regarded  as'an  effective  homogeneous 
material  which  has  the  properties  of  the  composite.  At  Che  onset  of 
further  plastic  loading  Che  ligament  material  carries,  on  average,  the 
part  of  overall  ply  strain  given  by  ^  In  (10).  Again,  this  strain  Is  not 
uniformly  distributed  In  the  ligament.  However,  If  one  extends  the  argu¬ 
ment  leading  to  derivation  of  the  relaxation  surface  of  the  undamaged 
composite  to  the  present  case,  the  relaxation  surface  of  the  damaged  ply 
with  open  cracks  is  obtained  In  the  form 

2  T  T  2 

g  (e  -  8)  -  2G  (e  -  8)  (A  A  )  C(A  A  )(e  -  8)  -  k  -  0.  (15) 

8  "  "  m  "  "  "me  "o  *  ~mc  "o  " 


When  the  cracks  are  closed,  the  relaxation  surface  of  the  ply  Is  again 
given  by  (13),  but  adjustments  must  be  made  If  the  cracks  close  with  a 
relative  shear  displacement  which  cannot  be  adjusted  by  relative  sliding. 
The  distinction  between  open  and  closed  cracks  Is  determined  by  Che  In¬ 
equality  chat  assures  crack  opening: 


c  0 

e  -  e  -  €  >  0, 

22  22  22 


(16) 


In  Che  coordinate  system  of  Fig.  1 
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Fig.  2  shows  an  illuscratlon  of  relaxation  surfaces  (IS)  of  a  (0/90)g 
laminated  plate  subjected  to  ln~plane  uniaxial  load  from  0  to  500  MPa  and 
to  load  cycles  from  500  MPa  to  SO  MPa.  Note  that  a  steady  state  response 
la  obtained  after  the  second  cycle.  In  this  example,  the  matrix  follows 
the  Ziegler  hardening  rule,  details  of  the  calculation  were  described  by 
Wung  (1987).  The  surfaces  and  the  loading  path  are  plotted  In  the 
In-plane  strain  coordinates  622,  F33.  The  top  bars  Indicate  laminate 
strains,  which  are  equal  In  all  piles;  the  X3  direction  coincides  with  the 
fiber  orientation  In  the  0”  ply,  and  X2  Is  the  fiber  direction  in  the  90* 
ply.  In  this  Illustration,  the  relaxation  surface  equations  (14,15)  have 
been  transformed  Into  the  laminate  strain  coordinates.  The  ctacks  have 
been  added  after  completion  of  the  load  cycle,  and  the  relaxation  surfaces 
(15)  are  therefore  plotted  from  the  end  position  at  50  MPa.  The  dashed 
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Motion  of  ply  relaxation  surfaces  during  first  and  second  loading 
cycle,  and  expansion  of  the  surfaces  at  different  values  of  the 
damage  parameter  S. 
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line  indicates  the  crack  opening  condition  (16).  As  one  would  expect,  a 
larger  overall  strain  is  needed  to  cause  the  initial  yield  strain  in  the 
matrix  of  a  cracked  ply.  Therefore,  the  open  crack  branches  of  the  relax¬ 
ation  surfaces  expand  with  Increasing  crack  density.  Again,  the  scalar  6 
denotes  the  crack  density,  while  the  vector  8  defines  the  position  of  the 
center  of  the  relaxation  surface. 


4.  THE  DAMAGE  PROCESS 

Damage  development  in  metal  matrix  Couipoaltes  is  often  observed  only  under 
sustained  cyclic  loading  in  the  platlc  range.  The  damage  process  in  these 
systems  has  many  unusual  aspects  which  deserve  attention  since  they  affect 
the  subsequent  modeling  and  analysis. 

We  focus  on  a  laminated  place  which  consists  of  many  fibrous  layers 
which  at  all  times  remain  perfectly  bonded.  This  constrains  the  plies  to 
identical  in-plane  strains,  and  a  fairly  straightforward  procedure  based 
on  laminated  plate  theory  can  be  used  to  derive  overall  plate  properties 
from  Che  ply  properties.  In  what  follows  we  assume  that  all  plies  are 
made  of  the  same  composite  material,  the  only  distinction  between  Che 
plies  being  Che  In-plane  orientation  of  the  fiber  which  also  determines 
Che  orientation  of  cracks.  Of  course,  laminates  niade  of  dissimilar  plies 
can  be  considered  as  well.  The  plate  is  subjected  only  Co  in-plane 
mechanical  stresses  which  are  assumed  to  be  uniform,  and  also  to  a 
spatially  uniform  change  in  temperature. 

Each  ply  consists  of  an  elastic-plastic  fibrous  composite  medium;  the 
macroscopic  elastic  response  of  a  ply  is  described  by  (7),  with  Che  over¬ 
all  properties,  ^  in  the  undamaged  state,  and  L(B),  H(8),  given  by 
(S)  and  (6),  in  Che  presence  of  damage.  The  elastic  state  terminates  when 
Che  stress  state  in  any  ply  satisfies  Che  yield  condition  (13)  in  the  un¬ 
damaged  state,  and  also  when  (16)  is  not  satisfied  in  a  damaged  ply. 
Similarly,  (IS)  and  (16)  describe  the  yield  condition  of  a  damaged  ply. 
Normality  of  the  plastic  strain  increment  to  the  matrix  yield  surface 
guarantees  norntalicy  at  the  ply  and  also  at  the  laminate  level.  This 
opens  Che  way  to  formulation  of  constitutive  equations  in  the  plastic 
range. 

Under  cyclic  loading  of  a  metal  matrix  laminate  in  the  plastic  range, 
the  plastically  strained  plies  eventually  develop  a  certain  fatigue  damage 
state  which  is  dominated  by  transverse  cracks  of  the  type  shown  in  Fig.  I 
(Dvorak  and  Johnson  1980).  However,  other  types  of  damage  such  as  local¬ 
ized  delamination  at  the  fiber-matrix  interface  and  between  plies,  as  well 
as  fiber  splitting  by  the  transverse  cracks  are  frequently  observed.  Of 
course,  when  two  adjacent  plies  have  transverse  crack  systems,  Che  cracks 
will  Intersect  at  ply  boundaries.  In  fact,  if  the  laminate  is  made  of 
momolayers  reinforced  by  large  diameter  fibers,  such  as  boron  or  silicone 
carbide,  the  ply  boundaries  are  not  well  defined  and  cracks  in  one  ply 
will  extend  up  to  the  layer  of  Inclined  fibers  in  Che  next  ply.  Modeling 
of  the  various  types  of  damage  would  be  very  difficult.  However,  one  may 
account  for  the  additional  damage  by  Increasing  the  effective  density  of 
transverse  cracks  in  each  ply.  This  can  be  understood  as  addition  of  more 
cracks  of  width  2a  to  chose  which  are  physically  present,  or  as  a  magnifi¬ 
cation  of  the  width  of  the  existing  cracks.  The  second  alternative 
appears  to  be  more  plausible  as  it  may  well  represent  the  actual  reason 


for  an  increase  In  W(8)  in  (3),  which  is  the  extension  of  the  existing 
cracks  by  delamination  and  ply  boundary  crossing.  Then,  the  scalar  0  may 
no  longer  be  regarded  as  a  measure  of  the  actual  crack  density,  but  as  a 
damage  parameter  to  be  used  in  the  results  of  Section  2,  which  are  still 
supposed  to  predict  the  effect  of  damage  on  stiffness  and  on  averages  of 
local  fields. 

The  experimental  results  of  Dvorak  and  Johnson  (1980)  on  B/Al  lamin¬ 
ates,  as  well  as  subsequent  studies  by  Johnson  on  SlC/Tl  plates  suggest 
that  cyclic  plastic  straining  of  the  matrix  is  Che  principal  cause  of 
fatigue  damage  growth.  No  damage  is  typically  observed  if  the  laminate  is 
loaded  by  an  elastic  load  cycle,  either  in  the  initial  or  shakedown  state. 
This  argument  can  be  extended  to  damaged  laminates.  In  particular,  one 
can  assume  chat  all  damage  growth  will  terminate  if  the  laminate  reaches 
an  elastic  deformation  state.  As  Illustrated  by  the  example  in  Fig.  2, 
elastic  deformation  can  be  restored  under  an  initially  inelastic  cycle  of 
loading,  if  the  amount  of  damage  in  Che  plies,  and  Che  applied  strain 
cycle  cause,  respectively,  expansion  and  translation  of  Che  ply  relaxation 
surfaces  such  that  the  prescribed  load  cycle  or  the  corresponding  strain 
cycle  can  be  accommodated  within  the  new  elastic  range.  The  damage 
evolution  process  can  then  be  regarded  as  a  mechanism  Chat  Che  composite 
laminate  employs  to  reach  an  elastic  state.  In  this  new  state,  the 
originally  inelastic  part  of  the  total  strain  in  each  ply  is  accommodated, 
in  part,  by  the  strain  caused  by  opening  of  the  cracks. 

Viewed  from  a  different  perspective,  the  damage  process  can  be  thought 
of  as  a  part  of  a  shakedown  mechanism  in  the  composite  laminate.  Accord¬ 
ing  to  Che  static  or  Melan  shakedown  theorem,  an  elstlc-pla  .Ic  solid  or 
structure  will  shake  down  if  any  admissible  residual  stress  field  can  be 
found  such  that  its  superposition  with  the  stresses  caused  by  the  applied 
loading  will  not  violate  the  yield  condition  anywhere  in  Che  solid.  In 
ocher  words,  the  laminate  will  shake  down  if  a  subsequent  yield  surface, 
or  its  relaxation  surface  counterpart  in  the  strain  space,  can  be  found 
which  contains  the  applied  load  or  strain  cycle.  Of  course,  shakedown  can 
take  place  only  if  the  structure  is  loaded  within  its  failure  envelope,  and 
if  early  collapse  by  incremental  plastic  straining  can  be  prevented.  That 
is  usually  the  case  in  laminated  places  where  the  elastic  fibers  support  a 
major  part  of  the  load  so  Chat  the  total  strains  are  small,  yet  substan¬ 
tially  larger  chan  the  initial  yield  strains  of  a  ply. 

It  is  useful  to  point  out  chat  cyclic  plastic  loading  of  Che  laminate 
creates  cyclic  plastic  strains  in  Individual  plies  which,  as  illustrated 
again  by  Fig.  2,  tend  to  reach  a  steady  state  after  relatively  few  cycles. 
On  the  other  hand,  the  plastic  deformation  cycle  also  promotes  low  cycle 
fatigue  damage  growth  which,  in  comparison,  proceeds  very  slowly.  Typi¬ 
cally,  several  thousand  cycles  may  be  needed  Co  cause  a  significant  change. 
One  may  then  expect  Che  relaxation  surfaces  Co  translate  much  more  rapidly 
Chan  expand.  The  direction  of  translation  should  be  such  as  Co  minimize 
Che  magnitude  of  plastic  work  per  cycle.  Under  such  circumstances,  Che 
relaxation  surfaces  will  tend  to  translate  into  such  most  favorable  posi¬ 
tions  which  will  assure  Chat,  the  amount  of  expansion  -  which  is  to  say 
extent  of  ply  damage  -  will  reach  only  the  minimum  amount  necessary  for  an 
elastic  accommodation  of  the  loading  program. 

Another  consequence  of  the  large  disparity  between  deformation  and 
damage  races  is  chat  in  each  damage  state,  or  at  a  particular  magnitude  of 
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6  under  a  prescribed  load  or  strain  cycle,  the  effect  of  past  deformation 
history  will  fade  very  quickly.  That  la  to  say  that  the  deformation  field 
In  the  laminate  at  a  particular  state  of  ply  damage  will  be  very  similar 
to  that  which  one  would  reach  If  this  amount  of  damage  was  Introduced  Into 
an  elastic  laminate  prior  to  the  application  of  the  corresponding  load 
cycle.  Of  course,  that  does  not  suggest  that  the  final  shakedown  state 
may  be  reached  by  purely  elastic  deformation  of  an  Initially  damaged 
medium.  We  recall  that  damage  causes  only  expansion  of  the  relaxation 
surfaces.  But  If  translation  Is  also  needed  due  to  the  minimum  damage 
requirement,  the  shakedown  state  under  the  prescribed  load  or  strain  cycle 
may  be  reached  only  after  several  cycles  of  plastic  straining  which  may  be 
needed  to  bring  the  respective  ply  relaxation  surfaces  to  their  final 
location.  However,  the  implication  la  that  the  final  deformation  state  In 
a  damaged  laminate  subjected  to  a  constant  load  cycle  Is  Independent  of 
previous  loading  and  damage  history.  Therefore,  In  arriving  at  a  final 
damage  state,  one  may  follow  any  convenient  path.  For  example,  the  damage 
analysis  of  the  typical  case  of  constant  amplitude  loading,  which  causes 
large  excursions  Into  the  plastic  range  during  isany  Initial  deformation 
cycles,  can  be  replaced  by  analysis  of  damage  caused  by  a  load  cycle  which 
expands  at  a  rate  comparable  to  that  of  damage  growth.  In  this  particular 
case  each  Increment  In  amplitude  Is  followed  after  few  cycles  by  a 
saturalon  damage  Increment  which  restores  elastic  straining  within  a  new 
shakedom  state. 

Of  course,  the  Independence  of  the  final  damage  state  on  the  loading 
path  would  be  expected  to  hold  only  In  elastic-plastic  systems  which  have 
the  ability  to  adjust  their  Internal  permanent  strain  state  In  response 
to  current  cyclic  loading  conditions.  Apart  from  the  above  arguments, 
available  results  obtained  In  tension-tension  fatigue  tests  of  B/Ai  metal 
matrix  coiqioslte  laminates  support  the  path  Independence  concept.  For 
example,  Johnson  (1982)  found  that  a  saturation  damage  state  at  a  certain 
maximum  load  amplitude  could  be  reached  by  dissimilar  loading  sequences. 
Another  support  for  this  concept  was  established  by  experiments  which 
showed  chat  the  amount  of  damage  In  a  laminate  was  determined  prlmsrily  by 
Che  load  amplitude  and  not  by  variations  of  the  mean  stress. 


5.  INCREMENTAL  SHAKEDOWN-DAMAGE  ANALYSIS 

We  now  present  some  results  which  Illustrate  certain  aspects  of  damage  de¬ 
velopment  In  B/At  laminates  subjected  to  cyclic  tension  loading.  Our  ob¬ 
jective  Is  Co  find,  for  several  different  load  amplitudes,  the  amount  of 
damage  In  each  ply  that  is  needed  Co  reach  a  shakedown  state  In  the  lamin¬ 
ate.  Of  course,  the  stiffness  loss  and  the  Internal  stress  distribution, 
particularly  the  fiber  stress,  are  also  of  Interest.  One  specific  lamin¬ 
ate  under  consideration,  of  0/90  layup,  was  already  discussed  In  connection 
with  Fig.  2.  In  addition,  a  similar  analysis  was  performed  for  a  0”  plate. 
The  final  load  cycle  we  wish  Co  reach  Is  from  >  50  MPa  to  «  500 
MPa  In  both  laminates.  The  actual  path  we  follow  starts  with  cycling  of 
the  laminate  to  a  steady  state,  as  In  Fig.  2.  Then,  while  Is  kept 

constant,  8^^,^  Is  reduced  to  bring  the  laminate  Into  an  elastic  state. 

Next,  Is  Increased  in  small  Increments.  After  each  Increment,  cracks 

are  Introduced  In  Che  plastically  deforming  plies  to  the  extent  needed  to 
acommodate  the  deformation  path  within  the  expanded  relaxation  surfaces. 
This  Involves  both  expansion  and  translation  of  the  surfaces,  which  cause 


a  new  plastic  strain  state  and  a  stress  redistribution  through  the  laminate 
Details  of  the  procedure  have  been  described  by  VTung  (1987).  As  an  example 
Fig*  3  shows  the  current  relaxation  surfaces  at  two  levels  of  Note 
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Fig.  3  Translated  and  expanded  relaxation  surfaces  of  a  damaged  laminate 

at  two  different  levels  of  S  . 

max 


Fig.  4  The  effect  of  sustained 
cyclic  loading  on  re¬ 
duction  of  axial  elastic 
modulus  of  two  B/At 
laminates.  Comparison 
of  prediction  with  ex¬ 
perimental  results 
obtained  under  constant 
load  an^litudes. 


ICINT  OF  INITIAL  ELASTIC  MODULUS 
AFTIN  a«lO*  CYCLES 


that  high  values  of  8  are  required  to  accoo^llsh  the  acconuDodatlon,  as 
anticipated  In  Section  4.  The  Increoental  expansion  of  the  loading  range 
continues  until  one  reaches  the  desired  final  magnitude  of  ®max* 

Fig.  4  shows  the  change  In  the  axial  elastic  modulus  caused  by  satu¬ 
ration  damage  In  the  two  laminates  as  a  function  of  the  applied  tension 
stress  range.  The  computed  results  are  plotted  together  with  experimental 
data  of  Dvorak  and  Johnson  (1980).  In  the  experiments,  the  saturation 
damage  state  was  defined  as  the  damage  state  after  2x10^  cycles  at  constant 
stress  aiq>lltude,  as  noted  In  Fig.  4,  but  actual  measurements  of  stiffness 
loss  Indicated  that  damage  usually  stabilized  after  5x10^  cycles.  Finally, 
Fig.  5  shows  the  computed  magnitudes  of  the  axial  stress  in  the  0*  layer 
fibers.  In  the  saturation  state  at  different  levels  of  This  stress 

change  has  been  plotted  up  to  «<iual  to  the  experimentally  observed 

endurance  limit.  Note  that  while  the  endurance  quite  different 

In  the  two  laminates,  the  terminal  fiber  stresses  are  neerly  identical. 

The  Implication  Is  that  fatigue  failure  occurs  In  these  composite  systems 
by  overloading  of  the  0*  fibers.  Of  course,  the  maximum  stress  Is  not  seen 
by  the  fibers  until  the  laminate  reaches  the  saturation  damage  state  at  the 
endurance  ®max*  In  the  Initial  stages  of  damage  development,  a  part  of  the 
load  la  carried  by  the  undamaged  off-axis  plies,  but  as  damage  grows  more 
stress  Is  absorbed  by  the  0*  fibers. 
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Fig.  5  Fiber  stresses  In 
zero-degree  plies 
after  damage-induced 
s  hakedown. 
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6.  CONCLUSION 

The  good  agreement  between  experimentally  measured  stiffness  changes  and 
the  predictions  derived  by  the  present  method.  Fig.  4,  has  also  been  repro¬ 
duced  In  an  application  of  the  method  to  fatigue  analysis  of  a  (0/i4S/90)s 
laminate.  We  recall  that  the  predictions  were  found  for  an  incrementally 
expanding  load  cycle,  while  the  experiments  were  performed  under  constant 
load  aoqilltudes.  This  seems  to  Indicate  that  the  approximations  Introduced 
In  the  oiodel  strike  a  reasonable  balance  between  simplicity  and  accuracy, 
and  that  the  damage  process  Is  actually  path  Independent.  Of  course,  a  more 
realistic  model  of  the  elastic-plastic  behavior  of  the  cracked  ply  would  be 
desirable.  However,  the  significant  results  are  the  identification  of  the 
saturation  damage  state  as  a  damage-induced  shakedown  mechanism,  and  of  the 
path-independence  of  the  damage  process.  These  aspects  of  the  problem 
could  be  exploited  in  design  of  a  more  efficient  approach. 
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